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Abstract 


The  results  of  investigations  of  the  diffusion  of  oxygen  in 
silicate  melts,  the  structure  of  melts,  and  the  redox 
equilibria  in  silicate  melts  are  presented  in  this  thesis. 
Five  separate  investigations  were  undertaken,  all  of  which 
were  directed  toward  gaining  an  understanding  of  the 
chemistry  of  silicate  melts,  and  in  particular  the  role 
played  by  oxygen  in  those  melts. 

The  diffusion  of  oxygen  was  investigated  in  six 
different  melt  compositions.  Three  of  the  compositions  were 
in  the  synthetic  system  diopside-anorthi te .  Oxygen  diffusion 
was  measured  in  those  melts  at  one  bar  pressure  by  means  of 
isotopic  exchange.  The  dif f usivit ies  determined  are  similar 
in  magnitude  to  divalent  cation  dif fusivities  and  obey  the 
compensation  law  for  the  diffusion  of  divalent  cations  in 
silicate  melts.  Those  observations  suggest  that  there  is  a 
significant  contribution  to  the  diffusivity  of  oxygen  from  a 
"cation-like"  diffusion  mechanism.  Eyring  model 
calculations,  using  the  oxygen  diffusion  data  and  melt 
viscosity  data  from  the  literature,  suggest  that  the  size  of 
the  average  diffusing  species  is  similar  to  that  of  the 
Si044"  anion. 

The  diffusion  of  oxygen  was  also  measured  in  three 
basaltic  melts  (an  olivine  nephelinite,  an  alkali  basalt, 
and  a  tholeiitic  basalt)  at  various  temperatures  and 
pressures  up  to  21  kilobars.  The  dif fusivities  were 
determined  by  monitoring  the  rate  of  reduction  of  Fe3+  to 
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Fe2  +  in  the  melts.  The  oxygen  di f f usivit ies  measured  are 
approximately  the  same  as,  or  slightly  greater  than, 
divalent  cation  di f f usivit ies  in  basaltic  liquids.  The 
diffusivity  of  oxygen  shows  an  abrupt  decrease  in  all  three 
melts  at  approximately  the  same  pressure  as  the  change  in 
the  liquidus  phase  from  olivine  to  pyroxene.  The  decreases 
in  oxygen  diffusivity  are  interpreted  as  being  related  to 
decreases  in  the  proportion  of  02'  anion  in  the  melts  during 
reactions  which  decrease  the  proportion  of  olivine  building 
units  and  increase  the  proportion  of  pyroxene  building  units 
in  the  melt.  The  results  suggest  that  oxygen  diffuses 
principally  as  the  02'  anion  in  basaltic  melts. 

The  structure  of  lead  orthosilicate  melts  was 
investigated  as  a  function  of  the  thermal  history  of  the 
melt  by  means  of  infrared  spectroscopy.  The  melts  were 
observed  to  become  increasingly  polymerized  with  decreasing 
rate  of  cooling  and  with  increasing  duration  of  isothermal 
soaking.  The  presence  of  silicate  anions  larger  than  Si044- 
in  the  melts  was  demonstrated,  which  requires  that  02‘ 
anions  also  be  present.  The  proportion  of  02"  anions  in 
Pb2Si04  melts  depends  on  the  thermal  history.  Therefore,  the 
ratio  of  non-bridging  oxygens  to  tetrahedral  cations  (NBO/T) 
also  depends  on  thermal  history. 

The  study  of  the  structure  of  Pb2Si04  melts  provides 
direct  evidence  of  the  presence  of  02'  anions  in  an 
orthosilicate  melt.  That  result  when  combined  with  the 
implications  of  the  diffusion  studies  suggests  that  02‘ 
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anions  are  present  in  much  more  polymerized  melts  and  that 
the  proportion  of  those  anions  in  a  given  melt  depends  on 
the  thermal  history  of  the  melt.  That  suggestion  requires 
that  melt  models  which  do  not  incorporate  02"  anions  be 
reassessed . 

The  two  investigations  of  iron  redox  equilibria  in 
basaltic  melts  provide  an  improved  predictive  equation  for 
the  oxidation  of  iron 
which  may  account  for 
(graphite  and  Fe°)  in 
too  oxidized  to  allow 


in  natural  melts  and  a  model  process 
the  presence  of  reduced  phases 
quenched  rocks  which  are  apparently 
those  phases. 
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I  .  Introduction 


A.  Introductory  remarks 

Geologists  and  geochemists  have  been  studying  silicate 
melts  for  most  of  the  twentieth  century.  Those  studies  have 
progressed  from  determinations  of  phase  equilibria  and 
physical  properties  in  simple  synthetic  systems  (e.g.  Bowen, 
1913  and  1915;  Kozu  and  Kani,  1935)  to  investigations  of 
melt  structure  (e.g.  Taylor  and  Brown,  1979;  Mysen  et  al . , 
1980d;  and  de  Jong  and  Brown,  1980)  and  the  thermochemistry 
of  silicate  melts  (e.g.  Carmichael  et  al . ,  1977;  Richet  and 
Bottinga,  1980;  and  Weill  et  al . ,  1980).  A  present  emphasis 
in  the  study  of  silicate  melts  is  the  study  of  kinetic 
phenomena  (e.g.  Watson,  1979;  Kirkpatrick  et  al . ,  1981;  and 
Karsten  et  al . ,  1982).  Such  studies  have  provided  a  great 
deal  of  information  about  the  crystallization,  structure  and 
chemistry  of  silicate  melts.  However,  one  aspect  of  the 
study  of  silicate  melts  that  has  received  comparatively 
little  attention  is  the  role  played  by  oxygen  in  silicate 
melts . 

This  thesis  is  an  attempt  to  fill  some  of  the  gaps  in 
our  knowledge  of  the  behavior  of  oxygen  in  silicate  melts. 

To  that  end,  five  more  or  less  independent  investigations 
were  undertaken.  The  results  of  these  researches  are 
presented  in  chapters  2  through  6  of  this  thesis.  Each  of 
chapters  2  through  6  is  written  as  an  independent  paper  and 
has  its  own  introduction  and  conclusions.  Therefore,  while 
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this  chapter  briefly  describes  the  rationale  for  the  study 
as  a  whole  and  to  some  extent  for  the  individual  chapters, 
the  detailed  rationale  for  each  study  is  left  to  the 
beginning  of  that  chapter.  Similarly,  the  concluding  chapter 
(chapter  7)  is  principally  a  summary  of  the  individual 
conclusions  of  chapters  2  through  6. 

B.  Rationale  behind  the  research 

Oxygen,  together  with  silicon  and  aluminum,  constitutes 
the  structural  framework  of  all  silicate  melts.  Oxygen 
accounts  for  more  than  90%  of  the  volume  in  all  geological 
silicate  melts.  Furthermore,  the  coordination  polyhedra  of 
all  the  cations  in  a  silicate  melt  are  formed  by  oxygen. 
Highly  charged  cations  (principally  Si 4  +  and  Al3  +  )  are  bound 
to  oxygen  "anions"  to  form  the  silicate  (or  aluminosilicate) 
anionic  units  which  define  the  structure  of  the  melt.  Other, 
less  highly  charged  cations  (e.g.  Mg2+,  Ca 2 + ,  Fe2+,  Na+,  and 
K + )  occupy  the  interstitial  "sites"  between  the  anionic 
units.  There  may  also  be  oxygen  anions  (02“)  occupying  the 
interstitial  "sites"  as  Bockris  et  al .  (1948)  suggested  on 
the  basis  of  their  electrical  conductivity  data.  However, 
the  models  presented  by  some  workers  (Mysen  et  al . ,  1980a, 
1980b,  1980c,  1980d,  and  1981)  do  not  allow  for  oxygen 
anions  in  melts  with  NBO/T  (ratio  of  non-bridging  oxygens  to 
tetrahedral  cations)  of  less  than  3.  Conversely,  the  models 
of  melt  structure  presented  by  other  workers  (e.g.  Masson  et 
al . ,  1970;  Hess,  1  97  1  ,  1  975,  and  1977;  and  Smart  and 
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Glasser,  1978)  require  the  presence  of  02"  anions  in  all 
melts.  That  difference  in  interpretation  is  significant  and 
must  be  resolved  if  the  structure  and  chemistry  of  silicate 
melts  are  to  be  understood. 

Oxygen,  because  it  accounts  for  most  of  the  volume  in 
silicate  melts,  is  in  a  unique  position  relative  to  the 
other  constituent  elements  of  silicate  melts.  The  chemical 
and  kinetic  behavior  of  oxygen  certainly  affects  the 
behavior  of  all  other  components  in  the  melt.  Consequently, 
investigations  of  the  chemistry  and  kinetic  behavior 
(diffusion)  of  oxygen  in  melts  can  provide  a  great  deal  of 
information  about  the  overall  chemistry  and  structure  of 
those  melts. 

Oxygen  diffusion  in  silicate  melts 

Comparison  of  the  diffusivity  of  oxygen  with  the 
dif f usivit ies  of  cations  can  provide  an  indirect  assessment 
of  the  structural  role  of  oxygen  in  melts.  If  oxygen 
diffusion  is  of  the  same  magnitude  as  cation  diffusivity  it 
may  be  concluded  that  the  mechanism  of  oxygen  diffusion  is 
similar  to  that  for  cations  (i.e.  oxygen  diffuses  as  a  free 
anion  -  02").  Similarly,  if  the  magnitude  of  oxygen 
diffusion  is  close  to  that  of  the  tetrahedral  cations  (Si4+ 
and  Al 3  + ) ,  oxygen  diffusion  may  be  linked  to  tetrahedral 
cation  diffusion.  In  other  words,  oxygen  diffuses  by  the 
movement  of  silicate  (or  aluminosilicate)  anionic  units. 

That  implies  that  the  oxygen  in  the  melt  is  dominantly  bound 
to  the  tetrahedral  cations. 
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Oishi  et  al .  (1975)  and  Yinnon  and  Cooper  (1980)  have 
observed  that  oxygen  diffusion  in  silicate  melts  is 
inversely  proportional  to  the  viscosity  of  the  melt.  If  that 
correlation  is  valid,  models,  such  as  that  due  to  Eyring 
(Glasstone  et  a/.,  1941),  may  be  used  to  model  the  diffusing 
species  from  the  combined  oxygen  diffusion  and  viscosity 
data . 

Chapter  2  presents  the  results  of  a  study  of  the 
diffusivity  of  oxygen  in  three  melts  along  the  join  diopside 
(CaMgS i  20 6  )  -  anorthite  (CaAl2Si  208  )  at  one  atmosphere  total 
pressure.  That  system  was  selected  for  the  investigation 
because  the  melts  become  increasingly  polymerized  as  the 
anorthite  composition  is  approached.  In  addition,  there  are 
good  viscosity  data  for  the  system  (Scarfe  et  al . ,  1980). 
Thus,  the  system  diopside-anorthi te  is  ideal  for  an 
investigation  of  the  relationships  among  melt  structure, 
viscosity,  and  oxygen  diffusivity. 

Pressure  can  have  a  considerable  effect  on  the 
structure  of  silicate  melts,  as  has  been  illustrated  by  many 
workers  (e.g.  Waff,  1975;  Kushiro  et  a/.,  1976;  Fujii  and 
Kushiro,  1977;  Scarfe  et  a/.,  1979;  Mysen  et  al . ,  1980d;  and 
Tyburczy  and  Waff,  1983).  Because  of  the  special  position  of 
oxygen  in  silicate  melts  (referred  to  above),  the  variation, 
if  any,  of  the  diffusivity  of  oxygen  with  increasing 
pressure  should  provide  considerable  information  about  the 
variation  of  the  structure  of  silicate  melts  with  pressure. 
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Chapter  3  is  a  study  of  the  effect  of  pressure  on 
oxygen  diffusivity  for  three  basaltic  melts.  The  melts 
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investigated  are  natural  rock  melts  and  were  chosen  because 
viscosity  and  phase  equilibrium  data  were  available  for 
those  or  very  similar  compositions.  Consequently,  the  study 
allows  a  further  investigation  of  the  apparent  correlation 
between  oxygen  diffusivity  and  viscous  flow.  The  study  also 
investigates  the  possibility  that  oxygen  diffusivity  may 
reflect  structural  transformations  in  the  melt  as  a  function 
of  pressure,  which  are  suggested  either  theoretically  (Waff, 
1975),  or  by  changes  in  the  liquidus  mineralogy  observed  for 
basaltic  liquids  (e.g.  Takahashi,  1980). 

Effect  of  thermal  history  on  melt  structure 

Gotz  et  al .  (1980)  and  Kirkpatrick  et  al .  (1981)  have 
shown  that  different  phases  can  be  made  to  crystallize  from 
a  given  melt  composition  by  varying  the  thermal  history  of 
the  melt.  Such  changes  imply  that  the  structure  of  the  melt 
depends  on  the  thermal  history.  Therefore,  an  investigation 
of  the  effect  of  thermal  history  on  the  structure  of  a 
suitably  chosen  melt  could  provide  some  very  useful 
information  about  the  structure  of  the  melt  and  could  also 
yield  some  valuable  insights  into  the  role  of  oxygen. 

Orthosilicates  define  an  end  member  as  far  as  the 
structures  of  geologically  significant  silicate  melts  are 
concerned.  The  other  end  member  is  Si02  melt.  All  geological 
melts  fall  between  those  two  extremes.  Indeed,  the  models 
put  forth  by  Mysen  and  co-workers  (Mysen  et  al . , 
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1980b,  1980c,  1980d,  1981  and  1982b)  utilize  the 
orthosilicates  as  the  low  silica  endmember.  Mysen  et  al . 
(1982b)  propose  that  melts  with  NBO/T  greater  than  3  contain 
only  Si044'  monomers,  Si2076"  dimers,  and  02~  anions.  If 
NBO/T  is  fixed  for  a  given  melt  composition  (as  Mysen  et 
al . ,  1982b  argue)  then  the  relative  proportions  of  those 
anions  will  also  be  fixed.  The  result  is  that  the  structure 
of  the  melt  should  be  invariant.  On  the  other  hand,  if  NBO/T 
varies  with  variations  in  the  thermal  history  of  the  melt 
then  the  relative  proportions  of  the  various  anionic  species 
may  also  change  as  a  function  of  thermal  history.  Such  a 
melt  could  contain  a  diverse  assembly  of  silicate  anions  in 
contrast  to  the  relatively  simple  melt  envisioned  by  Mysen 
et  al .  (1982b)  and  would  certainly  have  different  physical 
properties.  In  any  event,  an  orthosilicate  composition  melt 
is  not  likely  to  have  a  large  fraction  of  highly  polymerized 
silicate  anions  and  it  should  be  possible  to  identify  the 
contributions  of  the  various  anions  to  the  infrared  spectrum 
of  the  melt.  Therefore,  an  orthosilicate  melt  appears  to  be 
a  good  choice  for  an  investigation  of  the  effect  of  themal 
history  on  silicate  melt  structure.  Chapter  4  reports  the 
results  of  an  investigation  of  the  effect  of  thermal  history 
on  the  structure  of  lead  orthosilicate  melts. 

Iron  oxidation  state  and  oxygen  activity  in  melts 

Chapters  5  and  6  are  principally  concerned  with 
applications  of  the  oxygen  diffusion  data  to  petrologic 
problems.  Both  chapters  deal  with  the  oxidation  state  of 
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iron  in  silicate  melts.  Chapter  5  is  a  study  of  the 
relationship  of  the  oxidation  state  of  iron  in  melts  to 
temperature,  oxygen  fugacity,  and  bulk  composition.  That 
study  developed  as  an  outgrowth  of  chapter  3,  where  it  was 
necessary  to  determine  the  equilibrium  oxidation  state  of 
iron  in  three  melts  at  a  given  temperature  and  oxygen 
fugacity.  Oxygen  diffusion  data  were  used  to  determine  the 
required  time  to  equilibrate  the  melts.  Those  calculations 
pointed  out  the  likelihood  that  much  of  the  iron  oxidation 
state  data  in  the  literature  were  not  equilibrium  data. 
Consequently,  the  study  reported  in  chapter  5  developed  as 
an  attempt  to  expand  the  data  base  for  the  oxidation  state 
of  iron  in  melts  with  new  equilibrium  data.  Chapter  5 
presents  new  iron  oxidation  state  data  for  basaltic  liquids 
and  develops  and  discusses  a  predictive  model. 

Chapter  6  presents  the  results  of  a  study  of  the 
interrelationship  of  redox  reactions  involving  iron  in  a 
basaltic  melt  and  the  diffusion  kinetics  of  oxygen.  Loh  et 
al .  (1981)  have  shown  that  the  oxygen  activity  in  melts 

containing  polyvalent  cations  varies  greatly  with  changes  of 
temperature.  They  ascribe  the  changes  in  oxygen  activity  to 
rapid  redox  reactions  involving  the  polyvalent  cations.  The 
effects  they  observed  are  principally  kinetic  and  are  due  to 
the  difference  in  the  rates  of  cationic  redox  reactions 
(fast)  and  oxygen  diffusion  in  the  melt  (slow).  Chapter  6 
presents  the  results  of  a  study  of  the  feasibility  of  using 
an  oxygen  activity  probe  (Arculus  and  Delano,  1981)  to 
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directly  monitor  the  activity  of  oxygen  in  a  basaltic  melt 
during  heating  and  cooling. 


II.  Oxygen  Diffusion  in  Three  Silicate  Melts  along  the  Join 

Diops i de-An orthi te 


A.  Introduction 

Oxygen  self-diffusion  in  silicate  melts  and  glasses  has 
recently  been  a  very  active  area  of  research.  Table  1 
summarizes  the  bulk  of  the  recent  work  for  silicate  melts. 
The  reasons  for  this  activity  are  quite  numerous,  but  are 
centered  around  the  fact  that  oxygen  is  the  major 
constituent  of  silicates.  Consequently,  studies  of  the 
kinetic  behavior  of  oxygen  may  yield  a  great  deal  of 
information  about  the  structure  and  chemistry  of  silicate 
melts . 

In  an  attempt  to  systematically  investigate  the 
relationships  between  oxygen  self-diffusion  and  melt 
composition,  the  diffusivity  of  oxygen  has  been  measured  in 
three  melts  along  the  join  diops ide-anorthi te .  Those 
compositions  (Table  2)  were  chosen  because  samples  were 
available  upon  which  Scarfe  et  a/.  (1980)  have  made 
viscosity  measurements.  This  allowed  an  investigation  of  the 
validity  of  the  Eyring  model  relating  the  activation 
energies  of  oxygen  diffusion  and  viscous  flow.  Such  a  study 
is  suggested  by  the  close  correspondence  of  the  experimental 
data  of  Yinnon  and  Cooper  (1980)  and  the  Eyring  equation 
model  values  for  oxygen  diffusion  coefficients  in  a 
soda-lime-silica  melt. 
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Table  1.  Oxygen  diffusion  data  for  silicate  melts. 


Composition  (wt%)  T°C  D0  E  Ref. 


(cm2/s) 

( kcal ) 

Na20Ca0Si02 

1000-1300 

4.2x10' 

53 

1 

25K2075Si02 

850-1000 

4x10' 

62 

2a 

36K2064Si02 

700-1000 

2.4x10- 1 

46 

2b 

Basalt 

1280-1540 

2 . Ox  1 0  4 

90 

3 

Basalt 

1  160-1360 

5.9x10' 

54 

1  0a 

Andesite 

1260-1360 

6.0x10' 

60 

10b 

40CaO20Al 2O340SiO2 

1275-1435 

4.7x10s 

85 

4 

40CaO20Al 2O340SiO2 

1320-1540 

1.8x10' 

54 

5a 

40Ca020Al20340Si02 

1000-1450 

4.5x10° 

33 

6a 

45Ca01 5Al2O340SiO2 

1000-1450 

4.5x10° 

33 

6b 

25CaO 10Al2O365SiO2 

1000-1450 

4.5x10° 

34 

6c 

1 6Na  20 1 2Ca072Si 0  2 

800-1470 

4.5x10" 2 

38 

5b 

20Na  2O80SiO2 

1061-1395 

7.9x10- 2 

44 

5c 

37K2063Si02 

820-  902 

1.7x10' 

60 

7 

7K2068Pb025Si02 

578-  678 

3.2x10° 

72 

8 

97Pb03Si02 

1000-1400 

8.3x10° 

24 

6d 

94Pb06Si02 

1000-1400 

8.3x10° 

25 

6e 

90PbO1 0SiO2 

1000-1400 

8.3x10° 

26 

6f 

Na  20Ca0Mg0Si02 

700-1083 

6.3x10- 2 

50 

9a 

20K2020Sr060Si02 

680-1121 

8.5x10" 2 

52 

9b 

References:  1.  Doremus  (1960),  2.  May  et  a/.  (1974),  3. 
Muehlenbachs  and  Kushiro  (1974),  4.  Koros  and  King 
(1962),  5.  Oishi  et  a/.  (1975),  6.  Sasabe  and  Goto 
(1974),  7.  deBerg  and  Lauder  (1980),  8.  deBerg  and 
Lauder  (1978),  9.  Yinnon  and  Cooper  (1980),  10. 
Wendlandt  (1980) 


B.  Experimental  methods 

The  method  used  in  this  study  to  measure  oxygen 
diffusivity  is  that  of  Muehlenbachs  and  Kushiro  (1974).  The 
method  involves  exchange  of  oxygen  between  a  spherical 
sample  of  known  oxygen  isotopic  composition  with  an 
"infinite"  reservoir  of  gas,  also  of  known  oxygen  isotopic 


* 

. 


composition.  After  the  exchange  experiment  the  entire  sample 
is  analyzed  for  its  oxygen  isotopic  composition.  This 
configuration  and  analytical  method  allow  the  solution  of 
the  diffusion  equation  (Jost,  1960): 

U'  8 Of-<f’  6Og)/(tf1  80i-  80g)=6*-  2£n"  2exp( -n  2^2Dt/r 2 )  ( 1  ) 

n 

where  S'  80f  is  the  isotopic  composition  of  the  sample  after 
the  exchange  experiment.  S'  8Oi  is  the  isotopic  composition 
of  the  starting  material,  S'  80g  is  the  isotopic  composition 
of  the  sample  upon  equilibration  with  the  gas,  t  is  the 
duration  of  the  experiment,  r  is  the  radius  of  the  sample, 
and  D  is  the  diffusion  coefficient.  The  application  of 
equation  1  to  calculate  diffusion  coefficients  is  based  on 
several  assumptions  including:  1)  That  volume  diffusion  of 
oxygen  is  the  rate  limiting  process;  and  2)  That  convection 
does  not  occur  within  the  sample. 

The  results  of  this  study  confirm  the  observations  of 
Oishi  et  al .  (1975)  that  volume  diffusion  of  oxygen  is  the 
rate-limiting  process  in  long  duration  experiments  but  may 
not  be  in  short  duration  runs,  as  discussed  below. 

The  assumption  of  non-convection  within  the  sample  was 
tested  theoretically.  Rayleigh  numbers  were  calculated  for 
the  worst  case  experiments  (highest  temperature)  for  each 
composition  and  found  to  be  much  less  than  one.  As  a 
Rayleigh  number  of  greater  than  657.5  is  required  for 
convection  to  occur  (Chandrasekhar,  1961),  it  is  concluded 
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that  convection  did  not  occur  in  the  melts  used. 

Equation  1  is  a  series  expansion.  Consequently,  it  is 
important  to  know  how  many  terms  are  required  for 
convergence.  If  enough  terms  are  not  used  too  small  a 
diffusion  coefficient  will  result.  The  critical  parameter  in 
determining  the  required  number  of  terms  is  the  product  Dt 
(diffusion  coefficient  times  the  experimental  duration).  As 
Dt  increases  the  number  of  terms  required  for  convergence 
decreases.  A  small  exchange  fraction  requires  many  more 
terms  than  does  a  large  exchange  fraction.  To  alleviate  this 
problem  a  computer  program  was  used  to  fit  equation  1  to  the 
data  for  each  experiment.  That  program  adds  terms  until  the 
diffusion  coefficient  converges  within  0.1%.  The  minimum 
number  of  terms  used  was  arbitrarily  taken  as  100. 

The  starting  materials  were  aliquots  of  the  same 
glasses  used  by  Scarfe  et  al .  (1980).  In  order  to  increase 
the  accuracy  of  determination  of  the  exchange  fractions,  the 
starting  glasses  were  enriched  in  1 60  by  hydrothermal 
exchange  with  1 60  enriched  water.  The  glasses  were  finely 
crushed,  sealed  into  Pt  capsules  with  an  aliquot  of  1 60 
enriched  water,  and  exchanged  at  800°C  and  1  kilobar 
pressure.  The  exchanges  were  carried  out  for  at  least  seven 
days.  The  resulting  hydrated  glasses  were  finely  ground  and 
fused  for  1-2  hours  at  1450°C  in  an  Ar  atmosphere.  The 
glasses  were  then  reground  and  fused  for  an  additional  hour 
at  1450°C  in  argon.  The  resultant  glasses  were  coarsely 
crushed  and  stored  at  25°C. 
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Table  2.  Melt  compositions  used  (theoretical  values,  wt .  %) 
and  other  parameters. 


Cl 

C2 

C3 

Si02 

55.49 

50.33 

48.11 

A1203 

— 

15.39 

21.99 

CaO 

25.90 

23.49 

22.46 

MgO 

18.61 

10.79 

7.44 

E1 ( kcal ) 

39.2 

52.3 

51.9 

T°C ( 1 iq) 

1391 

1274 

1388 

1 80(i) (%o) 

-25.8 

-62.3 

-73.2 

1 80(g) (%o) 

+  23.5 

E  is  the  activation  energy  for  viscous  flow.  Both  E  and 
T(liq)  are  taken  from  Scarfe  et  al .  (1980). 


All  experiments  were  done  on  spheres  of  glass  (10-70 
mg)  suspended  on  weighed  loops  of  250  micrometer  Pt  wire. 

The  glass  spheres  were  made  by  placing  chips  of  the  starting 
glass  on  the  wire  loops  and  thoroughly  fusing  them  with  an 
oxy-acetylene  torch  to  yield  a  bubble-free  sphere.  The  torch 
fusion  time  was  held  to  a  minimum  (always  less  than  30 
seconds)  at  the  lowest  temperature  required  to  melt  the 
sample,  that  temperature  was  always  below  1760°C  (the 
melting  point  of  Pt ) .  A  reducing  flame  was  used  for  the 
fusion  to  minimize  oxygen  exchange  between  the  sample  and 
air.  For  each  of  the  starting  compositions,  a  sphere 
prepared  in  this  manner  was  analyzed  for  its  oxygen  isotopic 
composition.  Two  spheres  were  analyzed  for  composition  3  and 
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found  to  have  oxygen  isotopic  compositions  within  0.3%o  of 
one  another.  It  was  concluded  that  the  sample  preparation 
did  not  appreciably  alter  the  isotopic  composition  of  the 
starting  material.  That  conclusion  is  further  supported  by 
the  similarity  of  the  results  of  isothermal  experiments  in 
which  spheres  of  different  radius  were  used  (Table  3  and 
fig.  1).  Thus,  those  compositions  were  taken  as  the  initial 
isotopic  compositions  (/180i)  and  are  given  in  Table  2. 
Several  samples  of  the  exchange  gas  were  collected  both 
before  and  after  passing  through  the  furnace.  Isotopic 
analyses  of  the  gas  samples  were  within  0.2%o  of  each  other. 
The  mean  of  those  values  ( +23 . 5%o )  was  taken  as  cf1  80g  (Table 
2).  The  sphere-loop  combinations  were  weighed  and  the 
resultant  sample  mass  was  used  in  conjunction  with  melt 
densities  (calculated  using  the  method  of  Bottinga  et  al . , 
1982)  to  compute  the  sample  radii.  Radii  are  estimated  to  be 
accurate  within  ±2%.  The  samples  were  also  weighed  after  the 
experiments.  A  change  in  sample  weight  greater  than  2%  was 
taken  as  cause  to  reject  the  experiment. 

The  experiments  were  done  in  an  electrically  heated 
vertical  tube  furnace  equipped  with  a  gas-tight  alumina 
muffle  tube.  Before  each  experiment,  the  furnace  was  set  to 
the  run  temperature  and  allowed  to  equilibrate.  The  sample 
was  then  introduced  into  the  cold  zone  of  the  furnace 
(temperature  less  than  220°C)  and  the  furnace  was  sealed. 

The  flow  of  exchange  gas  (oxygen)  was  started  and  the  system 
was  left  to  allow  the  furnace  atmosphere  to  be  purged  of 
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air.  After  30  minutes  the  sample  was  lowered  into  the  hot 
zone  of  the  furnace,  thus  initiating  the  experiment.  At  the 
end  of  the  experiment  the  sample  was  removed  from  the 
furnace  and  placed  in  a  stream  of  cool  air.  The  samples 
cooled  to  ambient  temperature  (25°C)  within  20-30  seconds. 
The  samples  were  stored  at  110°C  until  oxygen  isotope 
analysis.  As  it  is  necessary  to  show  that  the  diffusion 
coefficients  are  independent  of  experimental  duration,  the 
lengths  of  replicate  experiments  at  each  temperature  were 
varied  by  a  factor  of  at  least  1.5  and  up  to  17. 

Oxygen  extraction  was  done  by  the  BrFs  technique 
(Clayton  and  Mayeda,  1963).  The  extracted  oxygen  was  reacted 
with  carbon  to  form  C02,  and  analyzed  on  an  isotope  ratio 
mass  spectrometer  (micromass  602).  The  isotopic  composition 
of  the  sample  is  given  as: 

S'  E0=1000( { ( 1 "O/1 60)smp/( 1 »0/1 60)std}-1 )  %o  (2) 

and  is  reported  relative  to  a  laboratory  standard.  The 
isotopic  analyses  are  considered  to  be  accurate  within  0.1%o 
except  as  noted  below.  As  the  bulk  compositions  of  the 

samples  were  known  the  theoretical  yield  of  oxygen  was 

\ 

calculated  and  used  as  an  indicator  of  the  quality  of  the 
oxygen  extractions  (Table  3).  The  oxygen  yields  were 
determined  by  transferring  the  C02  (produced  by  reaction  of 
the  oxygen  with  carbon)  into  a  known  volume  and  measuring 
the  pressure.  The  yields  ranged  from  24  to  105%.  Because  of 
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the  greater  than  100%  yields  observed  yields  as  low  as  93.4% 
were  accepted.  The  isotopic  compositions  of  the  samples  are 
summarized  in  Table  3. 

The  isotopic  analyses  and  the  radius  calculations  are 
the  primary  identifiable  sources  of  experimental  error.  A 
possible  source  of  error  in  the  isotopic  analysis  of  samples 
which  are  significantly  different  from  the  standard  used  is 
leakage  of  gas  between  the  sample  and  standard  portions  of 
the  mass  spectrometer.  Gas  leakage  causes  the  sample  to 
appear  to  be  closer  to  the  standard  than  is  actually  the 
case.  Furthermore,  the  error  increases  as  the  difference 
between  the  sample  and  standard  increases.  In  the  worst  case 
encountered  in  this  study  a  gas  leak  of  1%  would  result  in  a 
0 . 7%o  error.  Therefore,  it  is  possible  that  there  may  be 
errors  as  large  as  1.0%o  for  the  more  depleted  samples 
analyzed  in  this  study.  Taking  an  error  of  1%o  in  the 
isotopic  analyses  and  propogating  it  through  eqn .  1  yields  a 
worst  case  error  in  logD  of  approximately  0.5.  The  maximum 
scatter  in  replicate  analyses  is  approximately  0.5  logD 
units  which  supports  the  error  analysis.  However,  the  leak 
rate  for  the  mass  spectrometer  used  is  generally  assumed  to 
so  small  as  to  be  negligible  (k.  Muehlenbachs ,  pens.  comm.). 

C.  Results 

The  experimental  results  are  summarized  in  Table  3  and 
shown  graphically  in  figure  1.  Most  rate  processes  in 
silicate  melts  are  exponentially  related  to  the  inverse  of 
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Table  3.  Experimental  conditions  and  results. 


No. 

*-3 

o 

O 

Time 
( sec ) 

1  80 
(%o) 

r  (cm) 

-logD 

(cm 2/s ) 

yield 

(%) 

1  .  1 

1396 

5400 

-11.6 

0.172 

7.32 

96.7 

1.2 

1397 

12060 

3.1 

0.153 

7.05 

94.8 

1  .32 

1398 

7200 

-5.6 

0.  172 

7.10 

80.9 

1.4 

1431 

12000 

16.5 

0.110 

6.83 

99.4 

1.5 

1437 

9860 

14.8 

0.122 

6.83 

95.3 

1.6 

1438 

12960 

13.5 

0.132 

6.82 

101.6 

1.7 

1440 

7200 

14.6 

0.137 

6.49 

96.6 

1.8 

1473 

3600 

-4.2 

0.156 

6.82 

96.5 

1.9 

1473 

5400 

2.  1 

0.  149 

6.77 

96.8 

1.10 

1473 

8100 

15.3 

0.097 

6.82 

98.4 

2.  1  1 

1282 

3600 

-56.3 

0.131 

8.64 

97.7 

2.2 

1297 

61200 

12.8 

0.109 

7.50 

104.3 

2.3 

1298 

33180 

4.6 

0.127 

7.29 

n  .d. 

2 . 4  1 

1300 

3600 

-58.2 

0.192 

8.64 

93.8 

2.5 

1357 

29640 

9.2 

0.114 

7.24 

104.2 

2 . 6  1 

1361 

3600 

-54.5 

0.192 

8.64 

98.7 

2.7 

1419 

14400 

18.8 

0.115 

6.65 

103.0 

2 . 8  1 

1420 

3660 

-47.1 

0.157 

7.68 

98.4 

2.9 

1473 

8400 

8.6 

0.120 

6.66 

101.1 

2.10 

1475 

3600 

-29.6 

0.158 

6.95 

99.0 

3.  1  2 

1386 

10800 

-29.3 

0.109 

7.61 

26.3 

3.2 

1388 

7500 

-35.4 

0.090 

7.73 

101.6 

3.3 

1388 

7500 

-43.2 

0.107 

7.80 

95.7 

3.4 

1395 

12780 

-30.  1 

0.112 

7.64 

101.0 

3.5 

1395 

14640 

-27.6 

0.116 

7.61 

95.8 

3 . 6  2 

1430 

15780 

-26.2 

0.159 

7.34 

86.7 

3 . 7  2 

1430 

90360 

20.5 

0.132 

7.24 

73.0 

3.8 

1473 

3600 

-43.3 

0.116 

7.42 

105.2 

3.9 

1474 

5400 

-40.8 

0.127 

7.44 

99.4 

3  .  1  0  2 

1474 

7620 

-18.7 

0.115 

7.  14 

84.3 

-  Rejected  due  to  surface  effects  (see  text). 

-  Rejected  due  to  poor  oxygen  yield. 
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Figure  1.  Arrhenius  plots  for  the  composition  1  (Cl),  2 

(C2),  and  3  ( C 3 )  oxygen  diffusion  data.  The  equations  of 
the  regression  lines  are  given  in  the  text  as  equations 
4,  5,  and  6,  respectively. 
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the  absolute  temperature.  Such  a  relationship  is  termed  an 
Arrhenius  equation  and  is  given  for  the  case  of  diffusion  by 
equation  3: 


D  =  D  o  exp ( -E/RT )  (3) 

where  D  is  the  diffusion  coefficient,  D0  is  the  Arrhenius 
frequency  factor,  E  is  the  Arrhenius  activation  energy,  R  is 
the  gas  constant,  and  T  is  the  absolute  temperature.  The 
data  for  the  three  compositions  can  be  fit  to  Arrhenius 
equations  to  give  Arrhenius  activation  energies  for  oxygen 
diffusion.  The  resultant  equations  (shown  as  solid  lines  in 
f igure  1 )  are : 

D= ( 1 . 64  x  1 0 1 ) exp( ( -63000 ) /RT)  ( cm 2 / sec )  r  2  =  0 . 39  (4) 

D= ( 1 .35  x10' 1 )exp( (-47000) /RT)  (cm2/sec)  r2=0.75  (5) 

D=(1.29  x 1 0 ~ 2 ) exp ( ( -44000 ) /RT )  (cm2/sec)  r2=0.86  (6) 

for  compositions  1,  2,  and  3,  respectively.  The  activation 
energies  are  63±20  kcal  for  composition  1,  47±9  kcal  for 
composition  2,  and  44±6  kcal  for  composition  3.  The  errors 
quoted  are  calculated  from  the  regression  analyses  using  the 
method  of  Birge  (1932).  The  activation  energies  for 
compositions  1  and  3  are  based  on  data  from  all  of  the 
experiments  not  rejected  due  to  low  oxygen  yield,  but  the 
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activation  energy  for  composition  2  is  not.  The  composition 
2  experiments  form  two  trends  which  intersect  only  at  1782°C 
(fig  2.).  These  two  data  sets  may  be  delineated  in  terms  of 
experimental  duration,  the  value  of  Dt ,  or  exchange 
fraction.  The  experiments  having  a  longer  duration,  a  higher 
value  of  Dt ,  or  a  higher  exchange  fraction  fall  into  the 
upper  set  while  those  experiments  of  short  duration,  small 
exchange  fraction,  or  lower  Dt  value  fall  into  the  lower 
set.  It  is  apparent,  where  only  a  small  fraction  of  oxygen 
was  exchanged,  that  some  process  other  than  volume  diffusion 
is  the  rate-limiting  process.  Hofmann  and  Magaritz  (1977) 
noted  that  a  critical  value  of  Dt  had  to  be  exceeded  to 
decrease  errors  from  competing  transport  processes.  The 
composition  2  data  indicate  that  the  critical  value  of  Dt  is 
approximately  1  x  10” 4  cm2  for  the  experimental 
configuration  used  in  this  study.  Samples  having  smaller  Dt 
values  were  rejected  (Table  3). 

That  experiments  having  higher  exchange  fractions  give 
similar  D  values  regardless  of  the  experimental  duration 
indicates  that  the  competing  process  is  itself  of  short 
duration.  The  competing  process  is  most  likely  the 
establishment  of  a  surface  equilibrium  between  the  sample 
and  the  furnace  atmosphere. 

The  normal  test  in  diffusion  studies  for  surface 
effects  is  to  plot  the  product  of  the  diffusion  coefficient 
and  the  experimental  duration  (Dt)  versus  the  experimental 
duration  (t)  for  each  set  of  isothermal  experiments.  If 
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Figure  2.  Arrhenius  plot  for  the  composition  2  oxygen 

diffusion  results.  Experiments  having  Dt  values  less 
than  1x10“  4  are  shown  as  circles  and  those  having  larger 
Dt  values  are  shown  as  hexagons. 
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there  are  no  surface  or  other  competing  effects,  such  a  plot 
will  be  linear  and  will  pass  through  the  origin.  In  addition 
D  was  plotted  against  experimental  duration  (t)  for  each 
isothermal  series  of  experiments.  The  plots  for  composition 
2  show  that  volume  diffusion  is  not  the  rate  limiting 
process  for  the  low  Dt  experiments.  Therefore,  those 
experiments  were  rejected.  While  the  plots  for  compositions 
1  and  3  show  some  variation  from  the  ideal,  the  absolute 
variation  in  the  D  values  is  not  systematic  and  is  always 
small  (less  than  0.3  log10  units).  Consequently,  all  of  the 
data  were  used  for  compositions  1  and  3  as  noted  above. 

D.  Discussion 

Table  1  summarizes  the  bulk  of  the  previous  oxygen 
diffusion  measurements  for  silicate  melts.  The  results  of 
this  study  may  be  compared  to  those  data  either  from  the 
point  of  view  of  activation  energy,  or  on  the  basis  of  the 
absolute  magnitude  of  the  measured  diffusivity.  Before 
making  a  comparison  on  the  basis  of  activation  energies,  it 
is  desirable  to  consider  the  physical  meaning  of  the 
Arrhenius  activation  energy. 

Mechanisms  of  oxygen  diffusion  in  silicate  melts  have 
been  and  continue  to  be  enigmatic.  It  is,  however,  essential 
that  some  understanding  of  the  mechanism  of  oxygen  diffusion 
in  melts  be  achieved.  Without  such  an  understanding  the 
physical  meaning  of  the  Arrhenius  activation  energy  for 
diffusion  remains  abstruse. 
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Among  the  mechanisms  which  have  been  proposed  are: 
diffusion  of  free  oxygen  anions,  diffusion  by  combined  Si-0 
bond  rupture  and  free  anion  movement,  diffusion  by  oxygen 
transfer  during  reactions  among  large  silicate  anions  (with 
or  without  free  oxygen  anions),  and  diffusion  via  rotation 
and  translation  of  large  silicate  anions  without  bond 
rupture.  Mechanisms  involving  bond  rupture  have  been 
suggested  in  studies  where  large  activation  energies  have 
been  determined  (Muehlenbachs  and  Kushiro,  1974;  Koros  and 
King,  1962),  thereby  relating  the  activation  energy  to  the 
strength  of  the  Si-0  bond.  Oxygen  diffusion  via  free  oxygen 
anion  movement  has  been  proposed  by  those  workers  whose 
studies  yielded  low  activation  energies  (Sasabe  and  Goto, 
1974  and  Wendlandt,  1980).  In  those  cases,  the  activation 
energy  is  associated  with  the  energy  required  to  move  a  free 
oxygen  anion  from  one  "site"  to  another.  Studies  in  which 
intermediate  activation  energies  were  determined  tend  to 
support  diffusion  via  combined  mechanisms  (c.f.  deBerg  and 
Lauder,  1978  &  1980). 

The  Arrhenius  activation  energies  for  oxygen  diffusion 
determined  in  this  study  (C-1:  63  kcal/mole,  C-2:  47 
kcal/mole,  and  C-3:  44  kcal/mole)  are  of  the  same  order  of 
magnitude  as  those  determined  by  Doremus  (1960),  May  et  al . 
(1974),  Oishi  et  al .  (1975),  Wendlandt  (1980),  deBerg  and 
Lauder  (1978  &  1980),  and  Yinnon  and  Cooper  (1980)  for 
oxygen  diffusion  in  other  silicate  melts.  However,  the 
activation  energies  measured  in  this  study  are  considerably 
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smaller  than  those  determined  by  Muehlenbachs  and  Kushiro 
(1974)  and  Koros  and  King  (1962)  for  oxygen  diffusion  in 
natural  basalt  and  CaO-Al 203-Si02  melts,  respectively,  and 
are  greater  than  those  determined  by  Sasabe  and  Goto  (1974) 
for  a  variety  of  silicate  melts. 

The  activation  energies  determined  by  Muehlenbachs  and 
Kushiro  (1974)  and  Wendlandt  (1980)  for  basalt  melts  are 
very  different  (90  and  54  kcal/mole,  respectively).  The  90 
kcal/mole  reported  by  Muehlenbachs  and  Kushiro  (1974)  is  the 
largest  Arrhenius  activation  energy  for  oxygen  diffusion  in 
a  silicate  melt  yet  determined.  The  only  similar  value  is 
that  reported  by  Koros  and  King  (1962)  for  a  CaO-Al 203 -Si02 
melt  (85  kcal/mole).  However,  the  result  of  Koros  and  King 
(1962)  is  questionable  as  Oishi  et  al  .  (1975)  have  reported 
an  activation  energy  of  54  kcal/mole  for  oxygen  diffusion  in 
a  melt  of  the  same  composition.  Sasabe  and  Goto  (1974)  have 
also  studied  the  same  melt  composition  as  Oishi  et  al . 

(1975)  and  Koros  and  King  (1962)  and  report  an  activation 
energy  of  33  kcal/mole.  However,  Sasabe  and  Goto  (1974)  did 
their  measurements  via  oxygen  ion  electrode  and  as  Tran  and 
Brungs  (1980a&b)  have  shown,  such  methods  are  fraught  with 
difficulties.  In  addition  the  results  of  Sasabe  and  Goto 
(1974)  fall  well  off  the  trend  defined  by  all  other  oxygen 
diffusion  measurements  on  a  compensation  law  plot  (fig.  4) 
which  is  discussed  below. 

The  other  basis  for  comparison  is  the  magnitude  of  the 
diffusion  coefficients.  Figure  3  is  an  Arrhenius  plot  of  the 
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results  of  this  study  and  the  data  summarized  in  Table  1. 

The  data  define  a  range  of  oxygen  di f f usivi t ies  spanning 
approximately  2.5  orders  of  magnitude  at  a  given 
temperature.  The  only  data  which  fall  outside  that  range  are 
those  of  Sasabe  and  Goto  (1974)  which  are  from  two  to  three 
orders  of  magnitude  greater  than  the  rest  of  the  oxygen  data 
and  do  not  even  plot  within  the  range  of  fig.  3.  The  scatter 
in  the  data  combines  with  the  paucity  of  the  results  to  make 
it  impossible  to  derive  any  correlation  between  melt 
composition  and  oxygen  diffusivity. 

Notwithstanding  this  lack  of  a  correlation,  fig.  3  is 
useful  in  showing  that  Arrhenius  activation  energies  are  not 
neccessarily  a  good  basis  for  comparison  of  diffusion  data. 
The  results  of  Oishi  et  a/.  (1975)  and  Koros  and  King  (1962) 
for  CaO-Al 203~Si02  melts  illustrate  this  well.  Figure  3 
shows  that  the  results  of  those  two  studies  are  similar  when 
the  actual  magnitudes  of  the  diffusion  coefficients  are 
compared,  even  though  their  activation  energies  are  very 
different . 

The  actual  magnitudes  of  oxygen  diffusion  coefficients 
(Tables  1  &  3)  are  similar  to  those  of  various  divalent 
cations  (see  fig.  3).  In  some  melts  oxygen  has  been  reported 
to  be  the  fastest  moving  species  in  the  melt  (Oishi  et  al . , 
1975).  However,  the  di f f usivities  measured  in  this  study  are 
less  than  those  of  divalent  cations  by  up  to  an  order  of 
magnitude.  Physical  mobility  must  be  related  to  the 
mechanism  of  diffusion.  If  oxygen  moves  as  a  free  anion 
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Figure  3.  Arrhenius  plot  showing  the  results  of  this  study 
(Cl,  C2,  &  C3),  previous  oxygen  diffusion  measurements 
(numbers  refer  to  reference  numbers  as  given  in  Table 
1),  the  field  of  divalent  cation  di f f us i vi t ies  as  given 
by  Hofmann  (1980)  and  Watson  (1981)  (dotted  box),  and 
the  silicon  diffusion  data  of  Towers  and  Chipman  (1957) 
(dashed  line )  . 


(X 


1400  1200 


1000 


800 


ioyr®K 


0*65 


27 


(02~)  its  rate  of  movement  will  be  comparatively  rapid  like 
that  of  a  cation  of  similar  size  and  absolute  charge. 
Similarly,  if  oxygen  diffuses  via  bulk  translation  and 
rotation  of  large  silicate  anions,  its  mobility  is  expected 
to  be  small  compared  to  divalent  cations.  As  oxygen  mobility 
is  generally  similar  to,  or  slightly  less  than  divalent 
cation  mobility,  oxygen  may  diffuse  by  a  combination  of  the 
mechanisms  discussed  above.  The  results  of  this  study 
bracket  the  silicon  diffusion  data  of  Towers  and  Chipman 
(1957).  This  may  be  taken  as  contradictory  to  the  conclusion 
about  oxygen  diffusion  mechanisms.  However,  the  results  of 
Oishi  et  al .  (1975),  Koros  and  King  (1962),  and  Towers  and 
Chipman  (1957)  when  taken  in  conjunction  with  the  results  of 
Yinnon  and  Cooper  (1980),  show  that  comparisons  between  the 
relative  mobilities  of  oxygen  and  cations  for 
composi t ionally  diverse  melts  are  not  necessarily  valid. 

Thus,  the  measured  oxygen  diffusivity  may  be  thought  of 
as  representing  a  weighted  average  of  the  contributions  from 
the  various  mechanisms.  This  hypothesis  may  be  extended  to 
provide  an  explanation  for  the  slope  changes  observed  in 
plots  of  InD  versus  1 /T  (Yinnon  and  Cooper,  1980;  deBerg  and 
Lauder,  1978  &  1980;  and  Oishi  et  al . ,  1975),  by  postulating 
that  such  slope  changes  represent  zones  of  transition  from 
the  dominance  of  one  mechanism  to  another.  This  hypothesis 
is  also  compatible  with  the  fact  that  the  slope  changes  in 
InD  vs  1/T  plots  are  generally  at  or  near  the  glass 
transition  temperature,  where  the  translational  degrees  of 


■ 


28 


freedom  in  the  system  are  lost  with  decreasing  temperature. 

Compensation  law 

Winchell  (1969)  and  Winchell  and  Norman  (1969)  have 
proposed  that  all  ionic  dif f usivities  in  silicate  melts  and 
glasses  become  similar  at  approximately  1497°C.  This 
conclusion  is  based  on  the  observation  of  a  linear 
relationship  between  the  log  of  the  Arrhenius  frequency 
factor  and  the  Arrhenius  activation  energy  (log10D0  and  E, 
respectively:  see  equation  3).  Hofmann  (1980)  has  expanded 
the  data  base  used  by  Winchell  (1969)  and  Winchell  and 
Norman  (1969)  and  reports  the  relationship: 


(7) 


log10D0  =  1.08  x  10"4E  -  5.50 


where  D0  is  in  cm2/sec  and  E  is  in  cal/mole-deg .  The 
frequency  factors  and  activation  energies  determined  in  this 
study  are  within  experimental  error  of  the  line  represented 
by  equation  7,  when  the  errors  inherent  in  equation  7  are 
considered.  The  oxygen  diffusion  results  of  this  study  along 
with  those  in  Table  1  (excluding  Sasabe  and  Goto,  1974) 
define  a  line  essentially  identical  to  Hofmann's  result 
(figure  4)  with  the  equation: 


log , 0D0  =  1.27  x  1  0 '  4E  -  6.52  r2=0.76 


(8) 


From  equation  8  a  critical  temperature  of  1448°C  can  be 
calculated.  That  temperature  is  in  reasonable  agreement  with 
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Table  4.  Calculated  A,  and  (A2x3)  1/2  values  for  the  diffusion 
of  02 "  and  Si04  4  “ . 


02-  SiO 4  4 " 

Ai  7.13x10"9cm  2.42x10_8cm 

U2a  3  )1/2  3 . 59x  1  0  ‘  8cm  6.63x10-8cm 

the  critical  temperatures  calculated  by  Winchell  (1969): 
1234°C,  Winchell  and  Norman  (1969):  1504°C,  and  Hofmann 
(1980):  1 770°C.  Equations  7  and  8  and  fig.  4  show  that  the 
trend  of  the  oxygen  diffusion  data  are  compatible  with  the 
compensation  law  as  proposed  by  Winchell  (1969). 

The  work  of  Sasabe  and  Goto  (1974)  has  been  excluded 
from  the  compensation  law  calculation  as  their  data  fall 
well  off  the  trend  defined  by  all  other  oxygen  diffusion 
data  (fig.  4).  This  observation  in  conjunction  with  other 
considerations  discussed  above  leads  to  the  conclusion  that 
the  results  of  Sasabe  and  Goto  (1974)  may  be  in  error. 
Comparison  with  the  Eyring  model 

Yinnon  and  Cooper  (1980)  and  Oishi  et  al . ,  (1975) 
report  good  agreement  between  measured  oxygen  diffusion 
coefficients  and  those  calculated  using  the  Eyring  equation 
( equation  9 ) . 


D  =kT/?* 


(9) 


3C 


igure  4.  Compensation  law  plot  for  czyger.  diffusion.  The 
results  of  this  study  hexagons  ar.d  all  previous  o x y 
diffusion  data  ’circles  except  those  cf  Sasabe  ar.d  G 
('974)  (triangles)  are  used  to  deterr.ire  the 
compensation  law  equation  fcr  oxygen  dif fusion  solid 
line).  The  compensation  lav  for  cations  as  given  by 
Hofmann  (1580)  is  shown  as  a  dashed  line.  The  equations 
of  the  two  lines  are  given  in  the  text  as  equations  5 
and  7,  respectively . 
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The  parameters  in  equation  9  are  the  diffusion  coefficient 
(D) ,  the  Boltzmann  constant  (A),  the  melt  viscosity  (^ ) ,  the 
absolute  temperature  (T) ,  and  the  characteristic  particle 
translation  distance  (A).  D  values  were  calculated  from  the 
viscosity  data  of  Scarfe  et  al  .  (  1980)  and  an  assumed  value 
of  A  of  2.6  x  10‘8cm  (the  diameter  of  an  oxygen  anion).  The 
calculated  D  values  are  larger  than  those  measured  by 
approximately  an  order  of  magnitude.  This  difference  is 
considerably  larger  than  that  reported  by  either  Yinnon  and 
Cooper  (1980)  or  Oishi  et  a/.  (1975). 

The  use  of  the  Eyring  model  (equation  9)  to  model 
diffusion  requires  the  assumption  that  the  diffusing  species 
is  also  the  viscous  flow  unit.  The  failure  of  the  Eyring 
model  to  adequately  model  the  measured  oxygen  di f f usivi t ies 
of  this  study,  when  oxygen  anion  is  taken  as  the  viscous 
flow  unit  (A  =  2 . 6x 1 0 ~  8 cm) ,  implies  either  that  02-  is  not 
the  diffusing  species  or  that  02~  is  not  the  viscous  flow 
unit.  While  compositions  1  and  2  have  essentially  identical 
oxygen  di f f usivit ies ,  the  decrease  in  the  magnitude  of 
oxygen  diffusivity  between  compositions  2  and  3  is  as  the 
Eyring  model  predicts.  Also,  if  the  value  of  A  in  equation  9 
is  increased,  Eyring  model  di f f usivit ies  can  be  brought  into 
quite  close  agreement  with  the  dif fusivities  measured  in 
this  study. 

The  form  of  the  Eyring  equation  given  in  equation  9  has 
been  simplified  by  an  assumption  about  A.  Without  that 
assumption  the  equation  becomes  (Glasstone  et  al . ,  1941): 
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D  =  A ,  kT/j  *  2  2  3  (10) 

where  23  is  the  interatomic  distance  in  the  direction  of 
transport  and  and  A2  are  the  interatomic  distances 
perpendicular  to  the  transport  direction.  The  assumption 
that  the  product  (  A  ,  A  2  2  3)  gives  the  effective  volume  of 
the  diffusing  species  may  be  used  in  conjunction  with  the 
measured  viscosity  and  diffusivity  to  compute  values  for  A, 
and  ( A  2  X  3 )  1/2 .  The  results  for  composition  3  at  1400°C, 
assuming  that  the  diffusing  species  is  02“  ( V=9 . 2x 1 0 " 2 4cm3 ) 
and  Si04  tetrahedral  anion  ( V= 1 x 1 0 ~ 2 2cm3 ) ,  are  shown  in 
Table  4.  The  value  of  A,  determined  for  the  02~  case  is 
clearly  too  small  to  be  realistic  as  it  is  smaller  than  the 
radius  of  02‘.  Whereas,  the  values  determined  for  the  SiO* 
case  are  more  reasonable.  This  may  be  taken  as  implying  a 
larger  diffusion  unit  than  O2-  anions.  However,  complete 
cation  diffusion  data  for  this  melt  composition  are  required 
to  substantiate  this  supposition. 

E.  Conclusion 

The  results  of  this  study  are  compatible  with  previous 
oxygen  diffusion  measurements  for  a  diverse  selection  of 
silicate  melts.  However,  there  is  a  great  deal  of  scatter  in 
the  data  (fig.  3),  which  precludes  the  identification  of  any 
systematic  trends  or  relationships  between  oxygen 
diffusivity  and  silicate  melt  composition.  Some  studies 
(Oishi  et  al . ,  1975;  Wendlandt,  1980,  and  this  study)  do 
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show  a  rough  inverse  relationship  between  oxygen  diffusivity 
and  the  silica  content  of  the  melt.  However,  that 
correlation  is  far  from  conclusive. 

The  self-diffusion  coefficients  do  not  agree  with  the 
Eyring  equation  unless  mean  ionic  jump  distances  (A  ) 
considerably  larger  than  the  diameter  of  the  oxygen  anion 
are  assumed.  However,  the  sense  of  variation  of  the  actual 
dif f usivities  is  as  the  Eyring  equation  predicts.  Fitting 
the  data  to  the  Eyring  equation  by  varying  the  parameters 
allows  an  estimate  of  the  size  of  the  diffusing  species  to 
be  made.  For  the  compositions  studied,  the  average  size  of 
the  diffusing  species  is  too  large  to  be  an  oxygen  anion  but 
it  is  smaller  than  any  of  the  silicate  anions  present  in 
melts  except  the  monomer  ( Si04  4  " ) . 

Consideration  of  the  results  of  this  study  and  the  bulk 
of  previous  work  shows  that  oxygen  conforms  to  the 
compensation  law  for  cationic  diffusion  in  silicate  melts 
and  glasses.  The  range  of  oxygen  di ff usivities  was  also 
found  to  encompass  the  range  of  divalent  cation 
di f f usivities  in  silicate  melts. 

These  results  lead  to  the  rather  speculative  conclusion 
that  the  diffusion  of  oxygen  in  silicate  melts  may  involve  a 
contribution  from  a  cation-like  diffusion  mechanism  as  well 
as  contributions  from  the  diffusion  of  larger  structural 


units  in  the  melt. 


III.  Oxygen  Diffusion  in  Basaltic  Liquids  at  Elevated 


Pressures 


A.  Introduction 

Measurements  of  component  di f f usivit ies  have  proven 
very  useful  in  increasing  our  understanding  of  a  variety  of 
processes  in  silicate  melts.  Oxygen,  however,  is  unique 
because  it  is  the  major  component  of  all  silicates.  As  such, 
the  kinetic  behavior  of  oxygen  in  silicate  liquids  strongly 
affects  and  may  indeed  control  the  kinetic  and  chemical 
behavior  of  the  other  components  in  the  melt.  In  this 
respect,  Dunn  (1982)  and  Oishi  et  al , .  (1982)  have  shown 
that  oxygen  diffusion  data  can  be  very  useful  in 
understanding  the  structure  of  silicate  melts. 

The  data  set  for  oxygen  diffusion  in  silicate  melts  is 
quite  restricted  (Dunn,  1982).  The  bulk  of  the  data  is  for 
simple  systems  which  are  of  limited  applicability  to 
geological  problems.  With  the  exception  of  Muehlenbachs  and 
Kushiro  (1974),  who  report  oxygen  diffusion  results  for  a 
basalt  liquid,  and  Wendlandt  (1980),  who  presents  data  for 
tholeiitic  basalt  and  andesite  melts,  there  are  no  oxygen 
diffusion  data  for  rock  melts.  Unfortunately,  their  results 
for  basalt  liquids  are  s igni f icanlty  different.  While  that 
difference  probably  reflects  differences  in  experimental 
methods  and  melt  composition,  it  is  important  that  the 
difference  be  resolved. 
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All  existing  oxygen  diffusion  data  have  been  determined 
at  one  atmosphere  pressure.  Oishi  et  al  .  (1975)  have  shown 
that  oxygen  diffusion  in  silicate  melts  is  inversely  related 
to  the  viscosity  of  those  melts.  Recent  work  on  the  pressure 
dependence  of  the  viscosity  of  basaltic  liquids  has  revealed 
that  with  increasing  pressure  the  viscosities  of  those  melts 
either  remain  essentially  constant  (olivine  nephelinite  - 
Scarfe,  1981),  decrease  slightly  (1921  Kilauea  tholeiite  - 
Kushiro  et  al . ,  1976),  or  decrease  significantly  (alkali 
basalt  -  Scarfe,  1981).  Those  results  imply  that  oxygen 
diffusivity  may  increase  with  increasing  pressure  in 
contrast  to  the  behavior  of  Ca  (Watson,  1979)  and  C02 
(Watson  et  a/.,  1982)  (the  only  other  species  whose 
diffusion  has  been  studied  as  a  function  of  pressure). 
Therefore,  the  diffusion  of  oxygen  in  a  series  of  basaltic 
liquids  was  investigated  as  a  function  of  pressure. 

The  compositions  chosen  for  this  study  (Table  5)  are  an 
olivine  nephelinite  from  Hamada ,  Japan,  an  alkali  olivine 
basalt  from  Kettle  River  B.C.,  and  1921  Kilauea  tholeiite. 
The  three  compositions  were  chosen  because  they  increasingly 
approach  silica  saturation  (nephelinite  to  tholeiite)  and 
there  are  both  viscosity  and  phase  equilibrium  data  for 
these  or  similar  compositions  as  a  function  of  pressure.  The 
three  compositions  allow  an  investigation  of  the  effects  of 
temperature,  pressure,  and  bulk  composition  on  oxygen 
diffusion  in  basaltic  liquids. 
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Table  5.  Composition  of  starting  materials. 


C-  1  1 

KR-  1  3 

1921 

Si02 

37.77 

49.63 

48.08 

TiO  2 

2.75 

2.43 

2.97 

A  1  2  0  3 

12.00 

14.63 

12.11 

MgO 

9.29 

5.91 

11.30 

CaO 

12.59 

9.38 

9.44 

FeO1 

2.08 

2.26 

2.44 

Fe  202 

13.46 

10.10 

10.05 

MnO 

0.22 

0.21 

0.18 

Na  20 

4.29 

4.06 

2.25 

k2o 

2.17 

1  .25 

0.51 

P  2  0  5 

2.37 

0.50 

0.28 

Total 

98.99 

100.36 

99.68 

1  -  Volumetric  determination. 

2  -  Determined  by  difference  between  probe  FeO  and  FeO1 . 


B.  Experimental 

The  experiments  were  done  by  reducing  previously 
oxidized  spheres  of  basalt  glass  in  a  solid  media 
piston-cylinder  apparatus  (Boyd  and  England,  1960).  The 
oxidation  state  of  the  samples  was  determined  before  and 
after  the  experiment  by  micro-titration  (Wilson,  1960).  The 
spheres  of  glass  were  packed  with  graphite  powder  into 
graphite  capsules.  In  this  experimental  configuration, 
oxygen  diffuses  out  of  the  spherical  samples  into  an 
effectively  infinite  reservoir  of  graphite.  Analysis  of  the 
entire  sample  after  the  experiment  allows  the  solution  of 
the  diffusion  equation  given  by  Jost  (1960)  to  apply: 
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{FeO(f )  -FeO(e)  }/{FeO(  /)-  FeO(e)  }  =  6^'  2ln  -  2EXP  ( -n  2Tr2Dt/r  2 )  (11) 

n 

where  FeO (f)  is  the  wt .  %  FeO  in  the  sample  after  the 
experiment,  FeO (/)  is  the  initial  wt.  %  FeO  in  the  starting 
glass,  FeO (e)  is  the  equilibrium  wt .  %  FeO  in  the  sample,  D 
is  the  diffusion  coefficient,  t  is  the  experimental 
duration,  and  r  is  the  radius  of  the  sample. 

Equation  11  is  a  series  expansion  and  it  is  important 
to  know  how  many  terms  are  required  for  convergence.  If 
enough  terms  are  not  used,  too  small  a  diffusion  coefficient 
will  result.  The  critical  parameter  in  determining  the 
required  number  of  terms  is  the  product  of  the  diffusion 
coefficient  and  the  experimental  duration  (Dt).  As  D t 
increases  the  number  of  terms  required  for  convergence 
decreases.  To  alleviate  this  problem  a  computer  program 
using  a  large  number  of  terms  (>999)  was  used  to  fit  eqn .  11 
to  the  data  for  each  experiment  using  a  convergence  limit  of 
0.1%. 

The  approach  used  here  measures  the  chemical 
diffusivity  of  oxygen  rather  than  oxygen  self-diffusion. 

That  is  because  of  the  oxygen  chemical  potential  gradient 
imposed  by  placing  a  highly  oxidized  sample  into  extremely 
reducing  conditions  (discussed  below).  The  determination  of 
oxygen  diffusivity  by  monitoring  the  rate  of  reduction  of 
iron  in  the  sample  requires  that  oxygen  diffusion  be  the 
rate  limiting  process.  That  the  reduction  or  oxidation  of 
transition  metal  cations  in  silicate  melts  is  rapid  compared 
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to  oxygen  diffusion  has  been  demonstrated  by  several 
workers.  Semkow  et  al .  (1982)  have  directly  measured  the 
rates  of  oxidation  and  reduction  of  Ni,  Co,  and  Zn  in 
diopside  melts  and  found  them  to  be  fast  compared  to 
diffusion  (The  first  order  rate  constants  are  on  the  order 
of  10"2  cm/sec  at  1500°C).  Loh  et  a/.  (1981)  investigated 
the  oxidation-reduction  behavior  of  arsenic  in 
soda-lime-silica  melts  using  immersed  oxygen  fugacity 
probes.  They  found  that  immediate  oxidation  or  reduction 
occurred  on  cooling  and  heating,  respectively.  The  melts 
then  equilibrated  slowly  by  diffusion  of  oxygen  into  or  out 
of  the  bulk  melt.  Similarly,  Tran  and  Brungs  (1980b)  found 
that  the  rate  of  oxidative  equilibration  was  greatly 
enhanced  by  stirring  the  melt,  thus  decreasing  diffusive 
path  lengths. 

In  addition  to  the  arguments  above,  it  is  necessary  to 
show  that  the  diffusion  coefficients  are  independent  of 
either  sample  size,  or  experimental  duration.  Both  tests 
were  done  in  this  study  (Table  7),  but  the  test  of 
independence  from  sample  size  predominated  due  to  the  fact 
that  two  samples  were  generally  run  together  in  each 
experiment.  The  results  (Table  7)  show  that  the  diffusion 
coefficients  measured  are,  in  most  cases,  independent  of 
both  time  and  sample  size. 

Dunn  (1982)  and  Oishi  et  a/.  (1975)  noted  that  the  time 
required  to  establish  surface  equilibrium  caused  erroneous 
results  in  short  duration  oxygen  diffusion  experiments.  Even 
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though  the  experiments  done  in  this  study  were  short  (<  1 
hr.),  the  establishment  of  surface  equilibrium  is  not 
thought  to  be  a  problem  because  of  the  experimental  design 
and  procedure  used  in  this  study.  After  loading  into  the 
graphite  capsules  with  graphite  powder,  and  before  heating, 
the  samples  were  taken  up  to  a  pressure  3  to  5  kilobars 
above  the  run  pressure.  This  overpressuring  compresses  the 
sample  assembly,  seats  the  thermocouple,  and  forces  the 
sample  into  intimate  contact  with  the  graphite  powder, 
thereby,  assuring  rapid  attainment  of  surface  equilibrium 
upon  heating.  After  that  step,  the  experiment  was  initiated 
by  rapidly  heating  the  sample  to  the  run  temperature  (the 
time  from  900°C  to  run  T  was  generally  less  than  15 
seconds) .  The  experiments  were  terminated  by  shutting  off 
the  power.  The  water  cooled  apparatus  quenched  the  charge  at 
approximately  250°C  per  second.  Several  samples  examined 
optically  were  completely  glassy. 

The  starting  glasses  were  prepared  from  rock  powders  by 
multiple  fusion  in  air  at  1350°C  (nephelinite  and  tholeiite) 
or  1400°C  (alkali  basalt)  in  covered  Pt  crucibles.  The 
samples  were  crushed  between  fusions.  The  homogeneity  of  the 
starting  glasses  was  determined  for  the  major  elements  by 
electron  microprobe  and  for  FeO  by  at  least  four  replicate 
FeO  analyses  (Table  5). 

The  microprobe  analyses  were  done  on  the  University  of 
Alberta  ARL-SEMQ  automated  microprobe  by  wavelength 
dispersive  methods.  Data  reduction  was  done  by  the  method  of 
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Bence  and  Albee  (1968)  using  the  alpha  factors  of  Albee  and 
Ray  (1970).  FeO  analyses  were  done  by  the  method  of  Wilson 
(1960).  Samples  for  FeO  analysis  were  weighed  on  either  a 
Mettler  M5  microbalance  or  a  Cahn  Gram  electrobalance  to  a 
precision  of  ±5  micrograms.  The  titrations  were  done  with  a 
10  microliter  syringe.  These  procedures,  when  used  with 
reagent  solutions  as  described  by  Wilson  (1960),  allow  a 
precision  of  ±0.01  wt.  %  FeO  for  a  10  milligram  sample.  To 
test  the  accuracy  of  the  method  25  aliquots  of  BCR- 1  (8.80 
wt .  %  FeO;  Flanagan,  1973)  were  analyzed  and  gave  a  value 
for  FeO  of  8 . 66±0 . 2 1  wt.  %.  Several  batches  of  the  reagent 
solutions  were  used  over  the  course  of  this  study.  Those 
solutions  were  periodically  standardized  against  each  other 
as  well  as  against  the  initial  batch  of  ammonium  vanadate 
solution.  The  FeO  value  reported  here  for  BCR- 1  indicates 
that  the  split  had  become  slightly  oxidized.  However,  as  the 
FeO  contents  of  the  samples  in  this  study  need  only  be  known 
relative  to  one  another  rather  than  absolutely,  this  does 
not  cause  a  problem. 

Spheres  of  the  starting  glasses  were  prepared  either  by 
directly  grinding  chips  using  a  sphere  grinder  similar  to 
that  described  by  Bond  (1951),  or  by  gluing  chips  of  glass 
onto  loops  of  Pt  wire  and  fusing  the  chips  into  roughly 
spherical  shapes  at  either  1350°C  or  1400°C  in  air  (as  noted 
above)  and  then  grinding  the  spheroids  into  a  final 
spherical  shape.  Several  spheres  prepared  in  the  latter 
manner  were  analyzed  for  FeO  and  found  to  be 
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indistinguishable  from  the  starting  glasses.  The  latter 
method  of  making  spheres  was  adopted  due  to  the  long  times 
(10-15  days)  required  to  grind  irregular  chips  into  good 
spheres.  The  sample  spheres  were  ultrasonically  cleaned, 
weighed,  and  measured  with  a  micrometer  (±0.01  millimeters) 
before  use. 

The  experiments  were  done  in  a  solid  media 
piston-cylinder  apparatus  (Boyd  and  England,  1960)  using 
both  1/2"  (above  8  kbar)  and  3/4"  (8  kbar  and  below) 
assemblies.  In  order  to  ascertain  that  the  results  achieved 
with  the  1/2"  and  3/4"  assemblies  were  the  same,  two  runs 
were  made  for  the  nephelinite  at  16  kbar.  The  first  run 
(#67a  and  #67b  Table  7)  was  done  with  a  3/4"  assembly  and 
the  other  (#111  Table  7)  was  done  with  a  1/2"  assembly.  The 
results  for  samples  67a  and  111  are  indistinguishable  and 
the  errors  associated  with  all  three  results  overlap. 
Temperatures  were  monitored  with  Pt-Pt87Rh13  thermocouples 
without  any  pressure  correction.  Temperatures  are  estimated 
to  have  an  uncertainty  of  better  than  ±10°C.  The  pressure 
calibration  was  done  for  both  types  of  assemblies  by 
monitoring  the  melting  point  of  NaCl  by  DTA .  Pressures 
reported  here  have  an  uncertainty  of  ±0.5  kbar. 

Generally,  two  samples  with  different  radii  were  packed 
in  graphite  powder  in  a  graphite  capsule  (fig.  5).  The 
capsules  were  then  heated  to  red  heat  for  1-2  minutes  to 
drive  off  adsorbed  water  and  final  assembly  was  done  before 
the  capsule  cooled.  Samples  run  in  this  manner  retain  their 
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Figure  5.  Schematic  drawing  of  the  1/2"  sample  assembly 

showing  the  location  of  the  samples.  The  3/4"  assembly 
is  the  same  except  that  the  graphite  heater  and  the 
AlSiMag  end  plugs  are  tapered  to  reduce  the  temperature 
gradient  in  the  10mm  long  capsule  (Kushiro,  1976). 
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shape  (no  distortion  was  observed  for  any  of  the  samples  run 
in  this  study).  The  only  difference  in  the  samples  after  the 
experiments  is  that  they  are  coated  with  graphite.  The 
possibility  that  the  graphite  might  affect  the  FeO  analyses 
was  tested  by  mixing  equal  weights  of  BCR- 1  and  graphite 
powder  together  and  then  analyzing  for  FeO.  No  interference 
from  the  graphite  was  observed.  In  all  cases,  the  entire 
sample  was  analyzed  for  FeO  without  any  pre-analy t ical 
treatment . 

The  experimental  method  described  above  assures  that 
the  system  is  on  the  C-CO-CO 2-0 2  buffer  curve  (French  and 
Eugster,  1965).  That  buffer  curve  is  approximately  one  log 
unit  above  Fe-FeO  (Eugster  and  Wones,  1962)  for  the 
experimental  conditions  employed  in  this  study.  Furthermore, 
that  curve  is  well  below  FeO-Fe  30 4  (Ernst  1960  and  Eugster 
and  Wones,  1962).  The  relative  positions  of  the  three  buffer 
curves  are  insensitive  to  changes  in  pressure  and 
temperature  (fig.  6).  Thus,  the  equilibrium  state  of  the 
iron  in  the  samples  is  expected  to  be  essentially  100% 
ferrous  throughout  the  range  of  temperatures  and  pressures 
investigated.  This  was  verified  for  each  composition  by 
doing  a  long  duration  run  (minimum  4  hours)  and  analyzing 
the  sample  for  both  FeO  and  total  iron  (by  atomic 
absorption).  In  all  cases,  total  iron  (as  FeO)  was 
indistinguishable  from  FeO.  The  total  iron  content  (as  FeO) 
of  the  samples  was  taken  as  FeO(e)  for  use  in  the  diffusion 
calculation . 
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Table  6.  Volatile  contents  of  the  alkali  basalt  starting 
material  and  alkali  basalt  run  for  two  hours  at  1350°C 
and  12  kilobars. 


i  u 


H20 1 (wt .  %) 
CO 2 (wt.  %) 


#61 

<0.06 

0.018±0.002 


#74 

<0.06 

0 . 24±0 . 02 


1  -  wt .  %  H20  includes  the  total  pressure  of  non-condensible 
gas  (taken  as  H2  and  recalculated  as  H20)  in  addition  to 
gas  detected  as  water. 

Analyst  D.  Harris  using  the  method  of  Harris  (1981). 


Figure  6.  Log  of  oxygen  fugacity  versus  pressure  showing  the 

FeO-Fe  304  ( - )  ,  Fe-FeO  ( - ■)  ,  and  C-C0-C02-02 

( - )  oxygen  fugacity  buffer  curves  at  1280°C  (light 

lines)  and  at  1450°C  (heavy  lines). 
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Considerable  care  was  taken  to  assure  that  the  samples 
were  anhydrous.  This  was  verified  in  two  ways.  In  the  first 
case  a  sample  (#62  -  Table  7)  was  packed  into  a  Pt  capsule 
with  graphite  powder.  The  capsule  was  then  heated  to  red 
heat  and  welded  shut.  The  sealed  Pt  capsule  was  loaded  into 
a  1/2"  assembly,  surrounded  by  Fe  203  ,  and  run.  After  the 
experiment,  the  iron  oxide  packing  was  examined  both 
optically  and  by  x-ray  diffraction  and  found  to  consist  only 
of  Fe203,  implying  that  hydrogen  did  not  diffuse  into  the 
sample  from  the  talc  pressure  medium.  Furthermore,  the 
diffusion  coefficient  determined  from  that  experiment  agrees 
with  that  determined  in  another  experiment  done  in  the 
normal  fashion  (#63  -  Table  7).  In  the  second  case,  the 
volatile  contents  of  a  chip  of  the  alkali  basalt  starting 
material  and  a  chip  of  alkali  basalt  run  in  the  normal 
manner  at  1350°C  and  12  kbar  for  2  hours  were  determined  by 
the  method  of  Harris  (1981).  The  results  of  those 
measurements  (Table  6)  show  that  the  starting  material  is 
essentially  anhydrous  and  that  it  does  not  gain  an 
appreciable  amount  of  water  in  a  two  hour  run  (twice  the 
duration  of  the  longest  diffusion  experiment).  However, 
table  7  shows  that  the  C02  content  of  the  glass  increased 
from  0.018  wt.  %  to  0.24  wt.  %  during  the  2  hour  experiment. 
While  Eggler  and  Rosenhauer  (1978)  and  Eggler  et  al .  (1979) 
have  shown  that  dissolved  C02  can  affect  melt  structure, 
their  results  also  show  that  considerably  more  than  0.2  wt . 

%  are  required  to  produce  a  measureable  effect.  Thus,  the 
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effect  of  dissolved  C02  has  been  neglected  in  this  study. 

The  errors  reported  for  the  oxygen  diffusion 
coefficients  in  Table  7  were  calculated  by  propagating  all 
of  the  errors  due  to  the  determination  of  FeO,  sample  mass, 
experimental  duration,  and  sample  radius  through  eqn .  11. 
Those  errors  are  reported  at  the  one  standard  deviation 
level . 


C.  Results  and  Discussion 

The  results  of  this  study  are  summarized  in  Tables  7, 

8,  and  9  and  in  figures  7  through  11.  The  temperature 
dependence  of  oxygen  diffusivity  in  silicate  liquids  is 
adequately  described  by  an  Arrhenius  relationship  (eqn.  12), 
where  D  is  the  diffusion  coefficient,  D0  is  the  Arrhenius 
frequency  factor,  E  is  the  activation  energy,  R  is  the  gas 
constant,  and  T  is  the  absolute  temperature. 

D  =  D0EXP(-E/RT)  (12) 

Arrhenius  lines  were  determined  isobarically  for  the 
nephelinite  and  the  tholeiite  at  12  kilobars  (fig.  7)  and 
for  the  alkali  basalt  at  4,  12,  and  20  kilobars  (fig.  8). 

The  parameters  for  the  regression  lines  shown  in  figures  7 
and  8  are  given  in  Table  8.  The  activation  energies 
determined  in  this  study  span  the  range  of  previous 
determinations  of  activation  energy  for  oxygen  diffusion  in 
natural  rock  melts  (54  kcal  (Wendlandt,  1980)  to  90  kcal 
(Muehlenbachs  and  Kushiro,  1974)). 
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Table  7.  Summary  of  experimental  conditions  and  results. 


No. 

Alkali 

T°C 

basalt 

P ( kbar ) 

(KR-13) 

t (min ) 

FeO ( wt  %) 1 

r  (cm) 

-logD 
(cm2/sec ) 

72a 

1280 

4 

60 

9 . 66±0 . 03 

0.090 

6.5710. 10 

72b 

1280 

4 

60 

1 0 . 22±0 .  13 

0.072 

6.6310. 15 

71a 

1280 

12 

50 

9 . 24±0 . 05 

0.081 

6.6610.  1  1 

71b 

1280 

12 

50 

10.1 5±0 . 07 

0.070 

6.5910. 13 

70a 

1350 

4 

30 

10. 1 1 ±0 . 03 

0.092 

6.1410. 12 

70b 

1350 

4 

30 

9.  1  7±0 . 06 

0.074 

6.5310. 13 

50 

1350 

12 

30 

1  0 . 76±0 . 04 

0.077 

6.2110.07 

58 

1350 

12 

60 

1  0 . 94±0 . 03 

0.087 

6.3610.04 

57 

1350 

12 

240 

1 1 . 74±0 . 05 

0.076 

<6.70 

59 

1350 

20 

45 

8 . 78±0 . 04 

0.126 

6.3510.03 

73a 

1400 

4 

15 

9 . 2  3±0 . 06 

0.114 

5.8410. 13 

73b 

1400 

4 

15 

1 0 . 35±0 . 04 

0.080 

5.9010. 15 

51 

1400 

12 

120 

>11.60 

0.  144 

<5.84 

53 

1400 

12 

30 

1 1 .06±0.04 

0.125 

5.7410.06 

55 

1400 

12 

20 

1 0 . 09±0 . 04 

0.132 

5.7210. 10 

86a 

1400 

14 

20 

9 . 39±0 .  12 

0.073 

6.32  +  0.  14 

86b 

1400 

14 

20 

9 . 54±0 . 04 

0.065 

6.39  +  0.  12 

8  1  a 

1400 

16 

30 

9 . 58±0 . 05 

0.078 

6.4010.12 

81b 

1400 

16 

30 

9 . 29±0 . 07 

0.086 

6.38  +  0.  12 

79a 

1400 

20 

17 

9 . 53±0 . 08 

0.091 

6.0310. 12 

79b 

1400 

20 

17 

9 . 87±0 . 08 

0.066 

6.24  +  0.  12 

56 

1450 

12 

10 

8 . 80±0 . 04 

0.163 

5.48+0.07 

77 

1450 

20 

15 

9.2110.03 

0.122 

5.79  +  0.  1  1 

Nephelinite  ( C— 1 1 ) 


85 

1300 

12 

45 

1 1 .810.07 

0.087 

6.49+0 .09 

82a 

1350 

4 

30 

12.6410.09 

0.079 

6.28+0.09 

82b 

1350 

4 

30 

13. 6010. 05 

0.065 

6.22+0. 13 

65a 

1350 

8 

30 

13.0010. 10 

0.126 

5.78+0. 1 1 

65b 

1350 

8 

30 

12.93+0. 1 1 

0.138 

5.71+0.11 

62 

1350 

12 

30 

12.7310. 10 

0.095 

6.0810. 10 

63 

1350 

12 

31 

1 1 .9210.05 

0.115 

6.07+0.08 

67a2 

1350 

16 

30 

12.8510. 14 

0.101 

6.0110.12 

67b2 

1350 

16 

30 

12.88+0.06 

0.143 

5.7010. 17 

1113 

1350 

16 

18 

13. 1810.09 

0.071 

6.0310.20 

68a 

1350 

21 

30 

12.0310. 19 

0.078 

6.37+0. 1 1 

83 

1450 

12 

10 

12. 1 110.05 

0.096 

5.7010. 10 
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Table  7.  Summary  of  experimental  conditions  and  results 


1921  Kilauea  tholeiite 


No . 

CJ 

o 

E-* 

P( kbar ) 

t (min ) 

FeO(wt%) 1 

r  ( cm) 

-logD 
(cm2/sec ) 

88a 

1300 

12 

120 

10.6310.05 

0.096 

<6.43 

88b 

1300 

12 

120 

10.9410.05 

0.076 

<6.40 

89a 

1300 

12 

30 

10. 0310. 03 

0.086 

6. 1410.06 

89b 

1300 

12 

30 

10.25+0.05 

0.078 

6 . 1610.05 

1  12a 

1350 

4 

15 

10.03+0.03 

0.112 

5.6110.05 

1  12b 

1350 

4 

15 

10.  1410.01 

0.108 

5.6110.04 

1  13a 

1350 

6 

13 

10.2110.02 

0.101 

5.5910.05 

113b 

1350 

6 

13 

10.47+0.02 

0.088 

5.6110.06 

1  14a 

1350 

8 

15 

9.8510.05 

0.097 

5.7810.05 

1  14b 

1350 

8 

15 

10. 2910. 04 

0.078 

5.84+0.06 

90a 

1350 

12 

20 

9.9710. 10 

0.095 

5.8910.06 

90b 

1350 

12 

20 

10. 8010. 03 

0.072 

5.80  +  0.  10 

109a 

1350 

16 

30 

10. 0810. 03 

0.110 

5.9110.04 

109b 

1350 

16 

30 

10.3510.02 

0.100 

5.9110.04 

1  08a 

1350 

20 

30 

10. 1410.02 

0.101 

5.9710.04 

108b 

1350 

20 

30 

10.4210.02 

0.092 

5.9610.05 

105a 

1400 

12 

20 

10.2410.01 

0.110 

5.69+0.04 

105b 

1400 

12 

20 

10. 2410. 02 

0.090 

5.8610.05 

87 

1450 

12 

15 

10.2710.01 

0.114 

5.5210.05 

1  10 

1450 

12 

10 

9.92+0.03 

0.114 

5 . 4510 . 05 

1  -  FeO 

i  s  the 

wt .  %  1 

FeO  in  the 

sample  after  the 

experiment . 

2  -  Experiment  done  in  3/4"  assembly  at  16  kilobars  for 

comparison  with  experiment  #111  (see  note  3). 

3  -  Experiment  done  in  1/2"  assembly  at  16  kilobars  for 

comparison  with  experiments  67a  and  67b  (see  note  2). 
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Figure  7.  Arrhenius  plots  of  the  12  kbar  data  for  the 
olivine  nephelinite  (a)  and  the  tholeiite  (b).  The 
activation  energies  are  62±7  kcal  for  the  olivine 
nephelinite  and  51±4  kcal  for  the  tholeiite. 
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Figure  8.  Arrhenius  plots  for  the  alkali  basalt  data  at  4 
kbar'  (a),  12  kbar  (b),  and  20  kbar  (c).  The  activation 
energies  are  70±7  kcal  (4  kbar),  86±6  kcal  (12  kbar), 
and  7 1 ± 1 4  kcal  (20  kbar). 


v 


51 


Table  8.  Parameters  for  Arrhenius  equations  for  the  isobaric 
oxygen  diffusion  data. 


Composition 

P  ( kb ) 

log, oD0 
(cm2/sec ) 

E ( kcal ) 

r  2 

Nephelinite 

12 

2 . 23±0 . 97 

62±7 

0.94 

Alkali  basalt 

4 

3 . 2 1 ±2 . 5 1 

70±7 

0.81 

Alkali  basalt 

12 

5 . 4  5±0 . 8  1 

86±6 

0.95 

Alkali  basalt 

20 

3 . 23± 1 .87 

7  1  ±  1  4 

0.85 

Tholei ite 

12 

0 . 99±0 . 59 

5  1  ±4 

0.91 

The  parameters  above  were  determined  in  each  case  by  least 
squares  analyses  of  the  entire  data  set.  Errors  were 
calculated  by  the  method  of  Birge  (1932)  and  are 
reported  at  the  one  standard  deviation  level,  r2  is  the 
correlation  coefficient  for  the  regression  analysis. 


The  effect  of  pressure  on  oxygen  diffusion  in  the  three 
melts  was  investigated  isothermally  at  1350°C  for  the 
nephelinite  (fig.  9)  and  the  tholeiite  (fig.  11)  and  at 
1400°C  for  the  alkali  basalt  (fig.  10).  Figures  9,  10,  and 
11  show  that  the  pressure  dependence  of  oxygen  diffusion  in 
basaltic  liquids  is  not  simple.  The  olivine  nephelinite 
(fig.  9)  exhibits  an  increase  in  oxygen  diffusivity  with 
increasing  pressure  from  4  to  8  kilobars  followed  by  an 
apparent  sharp  decrease  between  8  and  12  kilobars.  Above  8 
kilobars  the  oxygen  diffusivity  continues  to  increase,  but 
the  pressure  dependence  is  smaller.  The  tholeiite  and  the 
alkali  basalt  show  sharp  decreases  in  oxygen  diffusivity 
with  increasing  pressure  at  6-8  kilobars  and  12-14  kilobars, 
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Figure  9.  Isothermal  (1350°C)  polybaric  diffusion  data  for 
the  olivine  nephelinite.  The  two  lines  yield  activation 
volumes  of  -39±3  cm3  (4-8  kbar)  and  -4±2  cm3  (12-16 
kbar).  The  preliminary  phase  liquidus  phase  relations 
for  C— 1 1  nephelinite  (Shimada  and  Scarfe,  unpublished) 
are  also  shown. 
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respectively  (figures  11  and  10). 

The  pressures  at  which  the  decreases  in  diffusivity 
occur  for  the  three  melts  correspond  well  to  the  pressures 
of  the  first  occurrence  of  pyroxene  as  a  liquidus  phase. 
Except  for  the  nephelinite,  phase  equilibrium  data  are  not 
presently  available  for  the  actual  melt  compositions 
studied,  but  there  are  data  available  for  composit ionally 
similar  melts.  Figure  10  shows  the  liquidus  relationships  of 
an  alkali-olivine  basalt  determined  by  Takahashi  (1980)  in 
addition  to  the  oxygen  diffusion  data.  Pyroxene  first 
appears  as  a  liquidus  phase  at  approximately  13.5  kilobars 
in  that  basalt.  Similarly,  fig.  11  shows  the  liquidus  phase 
relationships  for  1959  Kilauea  tholeiite  (Fujii, 
unpublished).  The  correlation  between  the  decrease  in  oxygen 
diffusivity  and  the  first  occurrence  of  pyroxene  as  a 
liquidus  phase  is  even  better  than  for  the  alkali  basalt. 

Except  for  the  transition  zones  noted  above  the 
pressure  dependence  of  the  diffusivity  of  oxygen  can  be 
described  by  an  Arrhenius  relationship  (eqn.  13): 

D  =  D  0EXP( -PUa/RT)  (13) 

where  D  is  the  diffusion  coefficient  at  the  pressure  of 
interest,  D  0  is  the  zero  pressure  diffusion  coefficient,  P 
is  the  pressure,  Va  is  the  activation  volume,  and  R  and  T 
are  the  same  as  in  eqn.  12.  The  solid  lines  in  figures  9, 

10,  and  11  are  least  squares  fits  of  eqn.  13  to  the 
diffusion  data.  The  activation  volume  for  the  high  pressure 
portion  of  the  nephelinite  data  does  not  include  either  the 
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Figure  10.  Isothermal  (1400°C)  polybaric  diffusion  data  for 
the  alkali  basalt.  The  two  solid  lines  give  activation 
volumes  of  -6±1  cm3  (4-12  kbar)  and  —  1 3 ± 4  cm3  (14  to  20 
kbar).  Also  shown  are  the  phase  relationships  for  an 
alkali  basalt  reproduced  from  Takahashi  (1980).  The 
coincidence  of  the  abrupt  decrease  in  oxygen  diffusivity 
with  the  first  appearance  of  pyroxene  as  a  liquidus 
phase  is  discussed  in  the  text. 
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Figure  11.  Isothermal  (!350°C)  polybaric  diffusion  data  for 
the  tholeiite.  The  two  solid  lines  give  activation 
volumes  of  -1.5±2  cm3  (4-6  kbar)  and  4±1  cm3  (8-20 
kbar).  The  phase  relationships  of  1959  Kilauea  tholeiite 
(Fujii,  unpublished)  are  also  shown.  The  coincidence  of 
the  sharp  decrease  in  oxygen  diffusivity  and  the  first 
appearance  of  clinopyroxene  as  a  liquidus  phase  is 
discussed  in  the  text. 
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data  point  at  21  kilobars  or  the  apparently  erroneous  point 
at  16  kilobars  (fig.  9).  The  point  at  21  kilobars  was 
excluded  because  the  preliminary  liquidus  curve  for  the 
nephelinite  (fig.  9)  shows  1350°C  to  be  below  the  liquidus 
at  21  kilobars.  The  coefficients  for  those  equations  are 
given  in  table  9.  It  is  of  interest  to  note,  except  for  the 
high  pressure  portion  of  the  tholeiite  data,  that  the 
activation  volumes  are  negative. 

Activation  volumes  for  diffusion  in  melts  have  been 
interpreted  as  being  related  to  the  volume  of  the  diffusing 
species  (Hamann,  1965;  Watson,  1979;  and  Watson  et  al . , 

1982)  and  also  to  the  sum  of  the  volume  of  the  diffusing 
species  and  the  volume  change  during  the  jump  from  "site”  to 
"site  (jump  volume)  (e.g.  Tyburczy  and  Waff,  1983).  Both 
interpretations  were  developed  principally  to  explain  the 
positive  activation  energies  observed  for  cation  diffusion 
and  neither  of  them  provides  an  explanation  for  the  negative 
activation  volumes  observed  in  this  study. 

Oxygen,  because  of  its  unique  structural  role  in 
silicate  melts,  is  unlikely  to  diffuse  in  strictly  the  same 
manner  as  do  cations.  Dunn  (1982)  speculates  that  oxygen 
diffuses  by  a  combination  of  the  movement  of  free  oxygen 
anions  (02~)  and  larger  silicate  anions.  Furthermore,  Dunn 
(1982)  also  suggests  that  with  changing  conditions  the 
relative  contributions  from  the  different  diffusing  species 
may  change  due  to  changes  in  the  relative  proportions  of 
those  species.  Koros  and  King  (1962)  and  Muehlenbachs  and 
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Kushiro  (1974)  propose  that  oxygen  diffuses  by  transfer 
between  different  anionic  units  during  reactions  involving 
the  breaking  of  Si-0  bonds.  The  activation  volumes  to  be 
expected  from  those  mechanisms  are  quite  different.  The 
activation  volume  for  Dunn's  suggested  mechanism  is  related 
to  the  average  volume  of  the  diffusing  species,  or  when  the 
proportions  of  the  diffusing  species  change,  to  the  sum  of 
the  average  volume  of  the  diffusing  species  and  the  change 
in  the  average  volume  due  to  changes  in  the  proportions  of 
species  of  different  sizes.  If  oxygen  diffuses  by  transfer 
during  reactions  among  various  anionic  species,  the 
activation  volume  may  be  related  to  the  volume  changes  of 
the  reactions. 

The  diffusion  coefficients  for  oxygen  determined  in 
this  study  are  approximately  the  same  as  the  diffusion 
coefficients  of  divalent  cations  in  basaltic  liquids  (see 
fig.  3,  chapter  2).  Therefore  it  seems  reasonable  to  suggest 
that  the  principal  diffusing  species  is  the  free  oxygen 
anion  (02').  With  the  exception  of  the  value  at  12  kilobars 
for  the  alkali  basalt  (Table  9),  the  activation  energies 
determined  in  this  study  are  too  small  to  be  associated  with 
the  rupture  of  either  Si-0  or  Al-0  bonds.  Therefore,  (with 
one  exception)  it  is  unlikely  that  there  is  a  large 
contribution  to  the  diffusivity  of  oxygen  in  the  melts 
studied  from  a  mechanism  involving  bond  rupture.  If,  as  Dunn 
(1982)  suggests,  the  diffusion  of  oxygen  is  by  the  movement 
of  02'  and  other  silicate  anionic  units,  then  the  observed 
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Table  9.  Parameters  for  Arrhenius  equations  for  the 
isothermal  oxygen  diffusion  data. 


Composition 

•-3 

o 

O 

P(  kb) 

1 og i o  D  o 

(cm2/sec ) 

l/a  (cm3 

r 2 

Nephelinite 

1350 

4-8 

6 . 75±0 . 05 

-3913 

0.98 

Nephelinite 

1350 

8-21 

5.5110.15 

-412 

0.92 

Alkali  basalt 

1400 

4-12 

5.9410.02 

-611 

0.91 

Alkali  basalt 

1400 

14-20 

6.9710. 18 

-1314 

0.61 

Tholei i te 

1350 

4-6 

5.6310. 10 

-1.512 

0.33 

Tholei i te 

1350 

8-20 

5.70+0.06 

411 

0.82 

The  parameters  above  were  determined  in  each  case  by  least 
squares  analyses  of  the  entire  data  set.  Errors  were 
calculated  by  the  method  of  Birge  (1932)  and  are 
reported  at  the  one  standard  deviation  level,  r2  is  the 
correlation  coefficient  for  the  regression. 


activation  volumes  can  be  related  to  the  changes  in  the 
average  volume  of  the  diffusing  species  as  the  proportions 
of  the  diffusing  species  change  with  increasing  pressure. 
Negative  activation  volumes  then  reflect  a  decrease  in  the 
average  size  of  the  diffusing  species  (an  increase  in  02'). 
As  smaller  anions  are  expected  to  be  more  mobile  than  larger 
anions,  such  a  change  in  anionic  speciation  is  compatible 
with  the  observed  increases  of  oxygen  diffusivity. 

Similarly,  the  observed  decreases  in  oxygen  diffusivity 
(figs.  9,  10,  and  11)  can  be  explained  by  a  decrease  in  the 
proportion  of  02‘,  which  imparts  a  positive  activation 
volume  as  is  required  for  those  changes.  This  model,  while 
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it  is  compatible  with  the  results,  remains  speculative 
because  the  identities,  relative  proportions,  and  partial 
molar  volumes  of  the  anionic  species  present  in  the  melts 
studied  are  unknown. 

The  observed  decreases  in  oxygen  diffusivity  (figs.  9, 
10,  and  11)  imply  that  the  melt  undergoes  a  change  of  some 
sort  at  8-12  kilobars  for  the  nephelinite,  at  12-14  kilobars 
for  the  alkali  basalt  and  at  6-8  kilobars  for  the  tholeiite. 
The  observation  that  those  pressures  are  close  to  the 
pressures  where  olivine  is  replaced  by  pyroxene  as  the 
liquidus  phase  (figs.  9,  10,  and  11),  implies  that  the 
change  in  the  melt  involves  an  increase  in  the  activity  of 
pyroxene  building  units  (short  chain  anions)  relative  to  the 
activity  of  olivine  building  units  (monomers)  for  the  two 
melts.  Tyburczy  and  Waff  (1983)  have  observed  that  the  slope 
of  the  log  of  electrical  conductivity  versus  pressure 
decreases  sharply  between  5  and  8  kilobars  for  a  Hawaiian 
tholeiite  melt.  They  postulate,  on  the  basis  of  their  data, 
that  Hawaiian  tholeiite  melt  undergoes  a  pressure  induced 
depolymerization  reaction  in  the  pressure  range  5-8 
kilobars.  The  viscosities  (Kushiro,  1976;  Scarfe  et  a/., 

1979  and  Scarfe,  1981)  and  densities  (Fujii  and  Kushiro, 
1977;  Scarfe  et  al . ,  1979;  and  Kushiro,  1980)  of  basaltic 
melts  do  not  show  sharp  changes  at  the  same  pressures  as  the 
observed  changes  in  oxygen  diffusivity.  Therefore,  the 
change  that  occurs  at  8-12  kilobars  for  the  nephelinite,  6-8 
kilobars  for  Kilauea  tholeiite  and  12-14  kilobars  for  the 


;  .  r,u 


. 


- 


60 


alkali  basalt,  does  not  involve  either  a  large  volume 
change,  or  a  major  depolymerization  of  the  melt. 

The  pressure  dependence  of  oxygen  diffusion  for  the 
melts  studied  is  sufficiently  varied  that  no  one  model  can 
adequately  explain  the  entire  data  set.  The  diffusion 
coefficients  determined  for  oxygen  in  this  study  are  of  the 
same  order  of  magnitude  as  cation  dif fusivities  in  basaltic 
liquids  (see  fig.  3,  chapter  2).  Therefore,  it  seems 
reasonable  to  speculate  that  the  principal  diffusing  species 
is  the  oxygen  anion  (O2-).  The  abrupt  decreases  in  oxygen 
diffusivity  at  approximately  the  pressures  where  olivine  is 
replaced  by  pyroxene  as  the  liquidus  phase  in  the  three 
melts  are  compatible  with  pressure  induced  depolymerization 
of  the  melt.  Because  of  the  lack  of  partial  molar  volume 
data  for  the  anionic  species  it  is  not  possible  to  write  a 
specific  depolymerization  reaction.  However,  the  changes  in 
the  phase  equilibria  suggest  that,  in  addition  to  decreasing 
the  oxygen  diffusivity,  the  reaction  that  occurs  also 
consumes  small  depolymer ized  anions  (olivine  building  units) 
and  produces  larger,  but  still  relatively  depolymer ized 
anions  (pyroxene  building  units).  Therefore,  it  is  suggested 
that  the  reaction  is  not  a  simple  depolymerization  reaction. 
Rather,  it  may  be  a  complex  reaction  in  which  large 
polymerized  anions  react  with  small  depolymer ized  anions  and 
oxygen  anions  to  produce  anions  of  intermediate  degrees  of 
polymerization.  A  hypothetical  reaction,  representative  of 
that  class  of  reactions  is  given  in  eqn .  14: 
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S i  2 0 5  2  ‘  +  S iO 4  4 '  +  0 2  ~  S i  3 0 ,  0  8  "  (14) 

where  Si2052"  represents  a  sheet  like  structure  and  Si30108- 
represents  a  short  chain  structure.  If  reactions  such  as 
reaction  14  occur,  they  result  in  a  net  depolymerization  of 
the  melt  while  simultaneously  decreasing  the  proportions  of 
small  relatively  mobile  oxygen  bearing  species.  Thus,  a 
reaction  like  reaction  14  can  account  for  both  the  observed 
changes  in  phase  equilibria  and  the  changes  in  oxygen 
di  f  f usivi ty . 

As  noted  above,  the  melts  studied  are  expected  to  show 
decreases  in  viscosity  with  increasing  pressure.  If  oxygen 
diffusion  is  inversely  related  to  viscosity  as  suggested  by 
Oishi  et  al .  (1975)  and  Dunn  (1982)  then  the  diffusivity 
should  increase  continuously  with  pressure  for  the  three 
melts  studied.  The  results  of  this  study  show  that  this  is 
not  the  case  and  lead  to  the  conclusion  that  the  mechanisms 
of  oxygen  diffusion  and  viscous  flow  are  not  related  for 
basaltic  liquids. 

The  parameters  in  Table  8  can  be  used  to  define  a 
compensation  law  for  oxygen  diffusion  in  basaltic  liquids. 
The  resultant  line  is  given  in  eqn .  15: 

log , 0D o  =  1 . 27x 1 0 “ 4E  -  5.64  r2  =  0.99  (15) 

Dunn  (1982)  has  determined  a  compensation  law  for  oxygen 
diffusion  in  a  variety  of  simple  silicate  melts  and  three 
natural  rock  melts.  Eqn.  15  is  essentially  the  same  as 
Dunn's  compensation  law  and  is  also  indistinguishable  from 
the  compensation  law  for  cationic  diffusion  in  silicate 
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melts  given  by  Hofmann  (1980).  The  compensation  law  result 
further  supports  the  contention  that  oxygen  diffuses 
principally  as  O2-. 

The  compensation  law  results  show  that  the  activation 
energies  determined  in  this  study  are  reasonable  in 
comparison  with  other  oxygen  diffusion  results.  The 
activation  energy  for  the  tholeiite  (51 ±4  kcal)  is 
indistinguishable  from  that  reported  by  Wendlandt  (1980)  for 
a  tholeiitic  basalt  (54  kcal).  While  the  physical  meaning  of 
activation  energies  in  liquids  in  not  clear  (Dunn,  1982),  it 
is  interesting  to  consider  the  variation  in  activation 
energy  of  the  alkali  basalt  as  a  function  of  pressure.  At  12 
kilobars  the  activation  energy  for  oxygen  diffusion  is 
considerably  larger  than  at  either  4  or  20  kilobars  (Table 
8).  Twelve  kilobars  is  very  near  the  transition  pressure  for 
the  alkali  basalt  (fig.  10).  The  increased  activation  energy 
at  12  kilobars  (relative  to  4  and  20  kilobars)  suggests  that 
the  activation  energy  is  related  to  structural 
transformations  in  the  melt.  Muehlenbachs  and  Kushiro  (1974) 
argue,  on  the  basis  of  their  90  kcal  activation  energy  for 
oxygen  diffusion  in  a  basalt  melt,  that  oxygen  diffusion 
involves  Si-0  bond  rupture.  Similarly,  the  86  kcal 
activation  energy  for  the  alkali  basalt  may  reflect  the 
occurrence  of  a  reaction  among  the  various  anionic  species 
at  12  kilobars  pressure. 
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D.  Conclusions 

This  study  has  shown  that  the  pressure  dependence  of 
oxygen  diffusivity  in  basaltic  liquids  is  quite  different 
from  that  shown  by  cations  in  simple  silicate  liquids 
(Watson,  1979;  Watson  and  Bender,  1980).  Variations  in 
oxygen  diffusion  coefficients  with  pressure  can  be  related 
to  the  liquidus  phase  mineralogy  and  by  extension  to  the 
anionic  constitution  of  the  melt.  As  such,  oxygen  diffusion 
measurements  provide  a  very  sensitive  probe  into  the 
variation  of  melt  structure  with  pressure. 

The  data  suggest  that  the  principal  diffusing  species 
is  the  O2 "  anion.  The  observed  changes  in  the  diffusivity 
with  pressure  can  be  attributed  to  changes  in  the  relative 
proportions  of  O2'  and  other  small  anionic  species  (Si04  4  ~ 
etc.)  during  pressure  induced  reactions  in  the  melt.  While 
it  is  not  possible  to  specifically  identify  those  reactions, 
the  correlation  of  the  changes  in  diffusivity  and  liquidus 
phase  equilibria  suggest  that  with  increasing  pressure  02' 
anions  and  olivine  building  units  are  consumed  in  reactions 
which  produce  pyroxene  building  units. 

This  study  has  shown  that  the  diffusivity  of  oxygen  in 
basaltic  liquids  is  very  sensitive  to  structural  changes  in 
the  melt.  Oxygen  diffusion  determinations  indicated  the 
presence  of  melt  structure  changes  that  do  not  show  up  in 
either  the  viscosity  or  density  data.  Therefore,  oxygen 
diffusion  measurements  constitute  one  of  the  most  sensitive 
probes  into  the  variation  of  melt  structure  with  pressure 
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IV.  The  Effect  of  Thermal  History  on  the  Structure  of 

Pb2Si04  Melts 


A.  Introduction 

Recent  work  on  the  effect  of  cooling  rate  on  the 
crystallization  of  Pb2Si04  melts  (  Gotz  et  al .  ,  1980)  and 
CaMgSi  20 6  melts  (Kirkpatrick  et  al .  ,  1981)  has  shown  that 
both  the  order  and  temperature  of  crystalline  phase 
occurrence  are  dependent  on  the  thermal  history  of  the  melt. 
Gotz  et  al .  (1980)  observed  that  different  lead  silicate 
phases  crystallized  from  Pb2Si04  melt  depending  on  the 
thermal  history  of  the  melt.  Similarly,  Kirkpatrick  et  al . 
(1981)  found  that,  depending  on  the  cooling  rate,  either 
forsterite  or  diopside  crystallized  first  from  melts  of 
diopside  composition.  In  both  cases,  the  phases  which 
crystallized  from  melts  with  different  thermal  histories  had 
different  structures.  The  crystallization  of  structurally 
different  phases  from  isochemical  melts  with  different 
thermal  histories  implies  that  the  structures  of  those  melts 
depend  on  thermal  history.  Such  differences  may  be  due  to 
changes  in  the  relative  proportions  of  the  silicate  anionic 
species  in  the  melt  without  altering  the  anionic  speciation 
(Mysen  et  al . ,  1980d),  or  to  changes  in  the  anionic 
speciation  of  the  melt  (Gotz  et  al . ,  1980). 

This  paper  presents  the  results  of  a  study  of  the 
effect  of  thermal  history  on  lead  orthosilicate  (Pb2Si04) 
melts.  There  are  several  reasons  for  choosing  to  study 
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Pb2Si04  melts  rather  than  more  geologically  significant 
compositions.  Orthosilicate  melts  have  a  non-bridging  oxygen 
to  tetrahedral  cation  ratio  (NBO/T)  of  4  (Mysen  et  al .  , 
1980d).  According  to  the  anionic  speciation  equations  of 
Mysen  et  al .  (1982b),  a  melt  with  NBO/T  =  4  should  have 
Si044~  monomers,  Si2076-  dimers  and  02"  anions  as  its  only 
anionic  species.  If  NBO/T  is  fixed  by  the  bulk  composition 
of  the  melt,  an  orthosilicate  melt  should  be  i sost ructural 
under  all  conditions.  The  results  of  Gotz  et  al .  (1980) 
suggest  that  this  is  not  the  case  for  Pb2Si04  melts.  If  the 
structure  of  Pb2Si04  melt  varies  as  a  function  of  thermal 
history,  then  NBO/T  for  that  melt  must  also  be  a  function  of 
thermal  history.  Consequently,  it  is  important  to  determine 
if  and  how  the  structure  of  Pb2Si04  melt  varies  with  thermal 
history . 

Pb2Si04  is  virtually  the  only  orthosilicate  whose  melt 
may  be  easily  quenched  to  glass.  More  geologically 
meaningful  orthosilicate  melts  (Mg2Si04 ,  Fe2Si04 ,  and 
Ca2Si04)  can  only  be  quenched  to  glass  with  the  aid  of 
exotic  rapid  quenching  techniques.  For  example,  Kusabiraki 
and  Shiraishi  (1981)  employed  a  splat  quenching  device  to 
quench  Fe2Si04  to  glass.  It  may  be  argued,  because  lead 
orhtosilicate  melts  are  easily  quenched  to  glass,  that  they 
are  not  representative  of  other  less  easily  quenched 
orthosilicate  melts.  However,  comparison  of  the  infrared 
(IR)  spectra  of  Pb2Si04  glasses,  prepared  by  rapid  quenching 
after  short  duration  fusions,  with  the  IR  spectrum  of 
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Fe2Si04  glass  (Kusabiraki  and  Shiraishi,  1981)  (fig.  12) 
shows  that  lead  orthosilicate  is  representative  of  other 
orthosilicate  melts.  However,  the  principal  Si-0  stretching 
envelope  is  broader  and  shifted  to  slightly  higher  frequency 
for  the  Fe2Si04  spectrum  as  compared  to  the  Pb2Si04 
spectrum.  As  is  discussed  later,  those  differences  imply 
that  vitreous  Fe2Si04  may  be  more  polymerized  than  vitreous 
Pb2Si04.  Thus,  conclusions  drawn  from  the  study  of  lead 
orthosilicate  melts  may  be  extended  to  other  orthosilicate 
compositions. 

Infrared  spectroscopy  can  be  a  very  powerful  probe  into 
the  structure  of  silicates  (e.g.  Billhardt,  1969;  Furukawa 
et  al . ,  1978,  1979;  Tarte,  1963).  Billhardt  (1969)  has  shown 
that  IR  spectroscopy  is  particularly  sensitive  to  the 
vibrations  of  Si-O-Si  bridging  bonds.  All  silicate  anions 
larger  than  the  monomer  (Si044~)  have  Si-O-Si  bridges. 
Therefore,  since  the  principal  aim  of  this  study  was  to 
determine  if  anionic  units  larger  than  either  the  monomer  or 
the  dimer  are  present  in  lead  orthosilicate  melts,  IR 
spectroscopy  was  chosen  as  the  structural  probe  for  use  in 
this  study. 

Due  to  the  difficulty  of  taking  IR  spectra  of  high 
temperature  silicate  melts,  glasses  formed  by  rapid 
quenching  of  melts  were  investigated.  Sweet  and  White 
(1969),  Sharma  et  a/.  (1978),  and  Seifert  et  al .  (1981)  have 
all  shown  that  glasses  formed  by  rapid  quenching  of  silicate 
melts  are  structurally  representative  of  the  melts  from 
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Figure  12.  Infrared  spectra  of  a.  vitreous  Pb2Si04  ( Pb— 1 ) 
and  b.  vitreous  Fe2Si04  (Kusubiraki  and  Shiraishi, 
1981  )  . 
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which  they  were  quenched.  However,  their  determinations  were 
not  for  orthosilicate  glasses  and  this  study  is  dependent  on 
the  validity  of  the  assumption  that  their  results  can  be 
extended  to  orthosilicates.  The  need  to  work  on  glass  was 
the  final  factor  in  the  choice  of  Pb2Si04  as  the  subject  of 
this  study. 

B.  Experimental 

The  Pb2Si04  glasses  were  prepared  by  cooling  at  various 
rates  though  the  temperature  interval  850  to  700°C.  Infrared 
spectra  were  taken  of  the  glasses  as  described  below.  All 
experiments  were  done  on  aliquots  of  two  5  gram  batches  of 
Pb2Si04  glass  prepared  from  reagent  grade  PbO  and  Si02.nH20. 
Experiments  Pb- 1  through  Pb- 1 0  were  done  on  aliquots  of  one 
batch  and  experiments  Pb- 1 1  through  Pb- 1 4  were  done  on 
splits  of  the  second  batch  (Table  10).  The  reagents  were 
dried  to  constant  weight  (PbO  -  200°C  and  Si02.nH20  - 
1400°C).  After  weighing,  the  reagents  were  thoroughly  mixed 
by  grinding  under  alcohol  in  an  agate  mortar.  The  resultant 
mixture  was  then  put  through  three  cycles  of  fusion  followed 
by  grinding.  All  fusions  were  done  at  850°C  and  had  a 
duration  of  two  hours.  Chips  of  the  glass  were  analyzed  via 
electron  microprobe  and  found  to  be  homogenous  and  on 
composition  (Table  10). 

The  electron  microprobe  analyses  were  done  on  the 
University  of  Alberta  ARL-EMX  microprobe  by  both  energy 
dispersive  (EDA)  and  wavelength  dispersive  (WDA)  methods. 
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Table  10.  Analyses  of  starting  materials  and  experimental 
products . 


No. 

Si02 

PbO 

Batch 

11.81 

88.19 

Pb-  1 

11.90 

88.10 

Pb-5 

11.92 

88.08 

Pb-6 

11.90 

88.10 

Pb-7 

11.81 

88.19 

Pb-8 

1  1.88 

88.12 

Pb-9 

11.79 

88.21 

Pb-  1  0 

11.92 

88.08 

Pb-  1  1 

10.66 

89.34 

Pb-  1  2 

10.99 

89.01 

Pb-  1  3 

11.01 

88.99 

Pb-  1  4 

10.84 

89. 16 

Analyses  were  done  on  chips  of 
as  described  in  the  text. 
100. 


glass  by  electron  microprobe 
Analyses  are  normalized  to 


Galena  and  quartz  were  used  as  standards  for  the  EDA 
analyses  and  the  data  were  reduced  with  the  computer  program 
EDATA2  (Smith  and  Gold,  1979).  The  glasses  analyzed  by  EDA 
and  a  lead  silicate  glass  standard  (EPS-22)  were  used  as  the 
standards  for  the  WDA  analyses.  The  WDA  data  were  reduced  by 
hand  using  atomic  number,  absorption,  and  fluorescence 
corrections. 

Approximately  200  milligram  aliquots  of  the  starting 
glass  were  placed  in  Pt  envelopes  which  were  crimped  shut, 
but  not  sealed.  The  capsules  were  left  open  to  air  to  retain 
the  lead  in  the  divalent  state  (Smart  and  Glasser,  1978). 
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All  of  the  experiments  were  done  in  a  Deltech  vertical  tube 
furnace.  Temperatures  were  monitored  by  a  Pt-Pt87Rh13 
thermocouple  placed  within  1-2  mm  of  the  sample. 
Thermocouples  were  calibrated  against  the  melting  point  of 
Au.  Cooling  rates  were  controlled  with  a  Love  cooling  rate 
controller.  In  all  cases,  samples  were  quenched  by  dropping 
the  Pt  envelopes  directly  from  the  furnace  into  water.  The 
resultant  quench  rate  is  estimated  to  be  in  excess  of 
500  °C/ sec . 

The  experiments  were  done  in  four  series.  The  first 
series  of  four  experiments  consisted  of  thirty  minute 
fusions  at  four  super 1 iquidus  temperatures  (760,  850,  964, 
and  1100°C)  followed  by  rapid  quenching.  The  second  series 
of  three  experiments  involved  a  thirty  minute  fusion  at 
850°C,  followed  by  controlled  cooling  at  one  of  three  rates 
(5.0,  0.7,  or  0.1°C/min.)  to  700°C  and  then  rapid  quenching. 
The  third  series  of  three  experiments  utilized  the  same 
fusion  durations,  temperatures,  and  cooling  rates  as  the 
second  series  of  experiments,  but  upon  reaching  700°C  the 
samples  were  held  for  48  hours  prior  to  quenching.  The 
fourth  series  of  four  experiments  involved  a  thirty  minute 
fusion  at  850°C,  followed  by  cooling  to  700°C  at  0.1°C/min. 
The  samples  were  then  held  for  different  times  (6,  12,  24, 
and  36  hours)  at  700°C  prior  to  quenching.  Experimental 
conditions  are  summarized  in  Table  1 1 .  There  was  no  evidence 
of  crystallization  at  700°C.  According  to  Gotz  et  al . 

(1980),  larger  degrees  of  undercooling  are  required  for 
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nucleation  and  crystal  growth.  A  chip  of  glass  from  each  of 
the  longer  duration  experiments  was  analyzed  by  electron 
microprobe  to  ascertain  that  the  samples  had  not  changed 
composition  during  the  course  of  the  experiments.  The 
analyses  of  the  samples  are  given  in  Table  10.  The  analyses 
of  the  experimental  products  are  all  within  analytical 
precision  of  the  theoretical  composition.  Experiments  Pb- 1 1 
through  Pb- 1 4  are  slightly  lower  in  Si02  than  Pb- 1  through 
Pb- 1 0 ,  but  the  difference  is  not  sufficient  to  acount  for 
the  structural  changes  noted  below. 

The  IR  spectra  were  taken  at  a  resolution  of  4  cm"1  on 
KBr  discs,  containing  1.0  wt .  %  glass,  using  a  computer 
controlled  Nicolet  7199  Fourier  transform  IR  spectrometer  in 
the  transmission  mode.  Fifty  scans  of  the  4000-400  cm-1 
spectral  region  were  collected  and  computer  averaged  for 
each  sample.  Such  averaging  enhances  the  signal  to  noise 
ratio  by  a  factor  of  approximately  7  relative  to  a  single 
scan  of  the  spectrum.  Replicate  spectra  were  taken  on 
several  samples  after  rotating  the  KBr  disc  to  assess  the 
effect  of  finite  sample  particle  size.  No  differences  were 
observed  between  replicate  spectra.  The  spectra  were 
transformed  from  transmittance  to  absorbance  prior  to 
printing.  As  there  are  no  features  attributable  to  silicate 
anionic  vibrations  above  approximately  1200  cm-1  ,  only  the 
400-1200  cm-1  portions  of  the  spectra  are  reported  here. 
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Table  11.  Summary  of  experimental  conditions  and  procedures. 


# 

T(i)°C  T(f)°C 

t  ( i  ) 

Rate 

t(f  ) 

t  ( t ) 

(min)  (deg/min) 

(min ) 

(min ) 

Pb-1 

760 

760 

30 

_ _ 

______ 

30 

Pb-2 

850 

850 

30 

— 

— 

30 

Pb-3 

964 

964 

30 

— 

— 

30 

Pb-4 

1  100 

1100 

30 

— 

— 

30 

Pb-5 

850 

700 

30 

5.0 

0 

60 

Pb-6 

850 

700 

30 

0.7 

0 

245 

Pb-7 

850 

700 

30 

0.1 

0 

1530 

Pb-8 

850 

700 

30 

5.0 

2880 

3060 

Pb-9 

850 

700 

30 

0.7 

2880 

3125 

Pb-10 

850 

700 

30 

0.1 

2880 

44  10 

Pb-1  1 

850 

700 

30 

0.  1 

360 

1890 

Pb-1  2 

850 

700 

30 

0.  1 

720 

2250 

Pb-1  3 

850 

700 

30 

0.1 

1440 

2970 

Pb-1  4 

850 

700 

30 

0.1 

2160 

3690 

T  (  i  )  - 

Initial  fusion 

temperature . 

T  (  f  )  - 

Final  experimental 

temperature . 

t  (  i  )  - 

Initial  fusion 

duration . 

Rate  - 

Cooling  rate. 

t(f)  - 

Duration  at  final 

temperature . 

t  ( t )  - 

Overall  length 

of  ' 

experiment . 

C.  Results 

The  results  of  this  study  are  in  the  form  of  IR 
spectra,  as  shown  in  figures  13-16  which  represent  the 
series  1-4  experiments,  respectively.  In  all  cases  the 
spectra  are  characterized  by  two  major  bands.  The  first 
band,  centered  at  480  cm'1,  is  of  roughly  gaussian  shape, 
but  shows  a  slight  high  frequency  shoulder  at  500-520  cm'1. 
The  other  band,  centered  around  880  cm'1,  varies 
considerably  in  form.  In  all  cases  there  is  a  definite  low 
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frequency  shoulder  in  the  630-680  cm-1  region.  Details  of 
the  spectra  for  each  series  are  given  below.  Vibrational 
mode  assignments  are  discussed  in  a  following  section. 

Series  one 

The  spectra  of  the  series  one  experiments  are  shown  in 
fig.  13.  It  is  apparent  from  fig.  13  that  the  spectra  of 
glasses  rapid  quenched  after  short  duration  fusions  at 
various  super 1 iquidus  temperatures  are  essentially 
identical.  This  is  in  agreement  with  results  obtained  from 
similar  experiments  for  glasses  of  diopside  composition 
(Mysen  and  Virgo,  1980;  Dunn  -  unpublished  data).  These 
spectra  are  also  qualitatively  indistinguishable  from 
spectra  of  Pb2Si04  glass  reported  by  Furukawa  et  al .  (1978) 
and  Gotz  et  al .  (1980). 

As  mentioned  above,  the  spectra  are  dominated  by  two 
bands  at  480  and  880  cm-1.  While  the  480  cm-1  band  is  very 
nearly  gaussian  in  shape  (there  is  a  slight  high  frequency 
asymmetry  which  implies  a  shoulder  around  500  cm-1),  the  880 
cm-1  band  is  definitely  a  composite  band.  The  880  cm-1  band 
has  a  high  frequency  shoulder  (approximately  900  cm-1  )  and 
a  low  frequency  shoulder  (approximately  800  cm-1  ),  both  of 
which  are  indicated  by  slope  changes.  In  addition,  there  is 
a  more  prominent  shoulder  on  the  low  frequency  side  centered 
around  650  cm- 1  . 

Series  two 

Fig.  14  shows  the  results  of  the  series  two  experiments 
(30  minute  fusions  at  850°C  followed  by  controlled  cooling 
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Figure  13.  Infrared  spectra  of  the  series  one  experiments. 
The  spectra  are  stacked  in  order  of  increasing  fusion 
temperature  (760°C  -  Pb- 1 ;  850°C  -  Pb-2;  964°C  -  Pb-3; 
and  1100°C  -  Pb-4).  Salient  features  discussed  in  the 
text  are  marked  by  arrows. 
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Figure  14.  Infrared  spectra  of  the  series  two  experiments. 
Samples  were  cooled  from  850  to  700°C  at  different 
cooling  rates  (5.0°C/min  -  Pb-5;  0.7°C/min  -  Pb-6;  and 
0.1°C/min  -  Pb-7).  Salient  features  discussed  in  the 
text  are  marked  by  arrows. 
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at  various  rates  to  700°C  and  then  rapid  quenching). 
Qualitatively,  the  spectra  are  very  similar  to  the  series 
one  results.  However,  there  are  notable  differences.  The  480 
cm-1  band  is  very  similar  to  that  in  the  series  one  spectra, 
except  that  the  high  frequency  shoulder  is  more  pronounced 
in  the  spectrum  of  the  sample  cooled  at  0.1 °C  per  minute. 

The  880  cm'1  band  is  significantly  different  from  that 
observed  in  the  series  one  experiments  and  also  varies 
systematically  in  form  as  a  function  of  cooling  rate.  In 
general,  the  high  frequency  band  is  characterized  by  a 
maximum  at  860-870  cm'1  ,  with  definite  low  frequency 
shoulders  centered  at  800  and  650  cm-1  ,  as  well  as  an 
asymmetry  on  the  high  frequency  side.  With  decreasing 
cooling  rate,  the  800  cm"1  shoulder  increases  in  intensity 
as  does  the  high  frequency  asymmetry.  The  650  cm-1  band 
shows  no  discernable  dependence  on  cooling  rate. 

Series  three 

The  results  of  the  series  three  experiments  (30  minute 
fusion  at  850°C  followed  by  controlled  rate  cooling  at 
various  rates  to  700°C  and  then  holding  at  700°C  for  48 
hours  prior  to  quenching)  are  shown  in  fig.  15.  The  spectra 
show  three  distinct  morphologies.  The  spectrum  of  the  sample 
cooled  at  5°C/min  is  indistinguishable  from  both  the  series 
two  experiment  that  has  the  same  cooling  rate  and  the  series 
one  experiments.  Similarly,  the  0.7°C/min  spectrum  is 
essentially  identical  to  its  series  two  counterpart. 

However,  the  0.1°C/min  experiment  shows  significant  changes 
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Figure  15.  Infrared  spectra  of  the  series  three  experiments. 
Samples  were  cooled  from  850  to  700°C  at  different 
cooling  rates  (5.0°C/min  -  Pb~8;  0.7°C/min  -  Pb-9;  and 
0.1°C/min  -  Pb- 1 0 )  and  then  held  at  700°C  for  48  hours. 
Salient  features  discussed  in  the  text  are  marked  with 
arrows . 
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as  a  function  of  thermal  soaking  at  700°C  relative  to  the 
corresponding  series  two  experiment.  The  480  cm-1  band  in 
that  spectrum  shows  two  maxima  (465  cm- 1  and  485  cm-1  )  and 
a  definite  shoulder  at  about  525  cm-1  .  The  largest 
differences  are  observed  in  the  880  cm'1  band.  That  band 
displays  three  distinct  maxima  at  920  cm'1  ,  880  cm-1  ,  and 
820  cm' 1  .  The  shoulder  at  650  cm" 1  appears  unchanged 
relative  to  other  series  1,  2,  and  3  spectra,  except  for 
experiment  PB- 1 0  where  it  is  of  increased  intensity.  In 
addition,  the  center  of  the  band  (determined  at  half 
maximum)  is  shifted  to  856  cm'1  ,  approximately  20  cm"1 
lower  in  frequency  than  observed  for  any  other  series  1,  2, 
or  3  experiments.  It  appears  that  thermal  soaking  at  700°C 
has  no  significant  effect  on  melts  cooled  at  0.7°C/min  or 
faster . 

Series  four 

The  results  of  the  series  four  experiments  (30  minute 
fusion  at  850°C  followed  by  cooling  at  0.1°C/min  to  700°C 
and  then  thermal  soaking  at  700°C  for  various  times  prior  to 
quenching)  and  the  series  three  experiment  that  was  cooled 
at  0.1°C/min  are  shown  in  fig.  16.  The  spectra  in  fig.  16 
are  stacked  in  order  of  increasing  thermal  soaking  time  with 
the  shortest  soaking  time  (Pb-11:  6  hours)  at  the  bottom  and 
the  longest  soaking  time  experiment  (Pb-10:  48  hours)  at  the 
top.  With  increasing  soaking  time  the  480  cm-1  band  shows  an 
increase  in  the  intensity  of  its  high  frequency  shoulder 
(500-520  cm"1).  Similarly,  the  high  frequency  shoulder  and 
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Figure  16.  Infrared  spectra  of  the  series  four  experiments. 
Samples  were  cooled  from  850  to  700°C  at  0.1°C/min  and 
then  held  at  700°C  for  varying  times  ( 6hr .  -  Pb— 11; 
12hr.  -  Pb-12;  24  hr.  -  Pb- 1 3 ;  36hr.  -  Pb- 1 4 ;  and  48hr. 
-  Pb- 1 0 ) .  Salient  features  discussed  in  the  text  are 
marked  by  arrows. 
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the  820  cm' 1  shoulder  in  the  880  cm" 1  band  also  increase  in 
intensity  with  increasing  soaking  time.  After  48  hours  of 
thermal  soaking  both  of  those  shoulders  can  be  resolved  as 
separate  bands.  In  all  cases  the  center  of  the  880  cm'1  band 
(determined  at  half  maximum)  is  shifted  approximately  20 
cm-1  to  855-860  cm"1  .  The  650  cm'1  shoulder  increases 
slightly  in  intensity  with  increasing  isothermal  soaking 
time.  In  addition,  the  band  width  (full  width  at  half 
maximum  -  FV7HM)  shows  a  continuous  increase  (from  270  cm"1 
to  400  cm"1  )  with  increasing  soaking  time. 

All  of  the  series  2,  3,  and  4  experiments  had  700°C  as 
their  final  temperature.  The  melting  temperature  of  lead 
orthosilicate  is  approximately  740°C.  Even  though  the 
samples  were  cooled  to  below  the  melting  temperature  and 
held  at  that  temperature  for  considerable  lengths  of  time, 
there  was  no  evidence,  either  in  the  form  of  visible 
crystals,  or  as  sharp  crystalline  bands  in  the  IR  spectra, 
that  any  crystallization  occured. 

In  summary,  the  effects  on  the  melt  structure  of  slow 
cooling  and  isothermal  soaking  are  similar.  Both  treatments 
result  in  an  increase  in  the  number  of  bands  in  the  high 
frequency  envelope,  but  the  effect  of  isothermal  soaking 
appears  to  be  larger  than  that  of  slow  cooling. 
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D.  Discussion 

Workers  such  as  Gotz  et  al .  (1980  and  1976),  Furukawa 
et  al .  (1978),  Smart  and  Glasser  (1978),  and  Hess  (1975) 
have  argued  on  the  basis  of  both  spectroscopic  and 
chromatographic  evidence  that  Pb2Si04  melts  contain  a 
diverse  assemblage  of  silicate  anions  as  well  as  Pb2+  and 
02~  ions.  Another  possibility  is  that  with  increasing 
duration  of  heating,  Pb2Si04  melts  begin  to  phase  separate 
into  PbO  and  Si02-rich  regions.  However,  if  such  phase 
separation  occurs,  it  is  on  a  very  small  scale.  Furukawa  et 
al  (1978)  report  no  evidence  of  phase  separation  at  a  scale 
greater  than  50  angstroms.  The  glasses  prepared  in  this 
study  are  both  optically  and  composit ionally  (microprobe) 
homogenous.  Therefore,  it  is  concluded  that  phase  separation 
does  not  occur  in  Pb2Si04  melts. 

The  silicate  anions  presumed  to  exist  in  Pb2Si04  melts 
include  Si044'  (Furukawa,  et  al . ,  1978;  Smart  and  Glasser, 
1978;  and  Gotz  et  al . ,  1976),  Si2076-,  Si3O108-,  Si40128“, 
and  larger  polysilicate  anions  (Smart  and  Glasser,  1978  and 
Gotz  et  al . ,  1976).  In  addition,  Worrell  and  Henshall  (1978) 
and  Furukawa  et  al .  (1978)  speculate  as  to  the  existence  of 
Pb-0  units.  In  as  much  as  the  Raman  active  bands  of 
crystalline  lead  oxides  occur  at  frequencies  below  400  cm-1 
(Worrell  and  Henshall,  1978),  it  seems  unlikely  that  any  IR 
active  bands  for  Pb-0  vibrations  will  fall  within  the 
spectral  range  (400-1200  cm-1)  examined  in  this  study. 
Consequently,  the  vibrational  modes  of  hypothetical  lead 
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oxide  will  not  be  considered  further.  The  contributions  made 
by  the  various  silicate  anions  to  the  IR  spectrum  of  the 
melt  may  be  estimated  by  considering  the  spectra  of  suitable 
crystalline  analogs.  Ideally,  the  crystalline  phases  used  as 
analogs  should  have  both  the  same  metal  cations  as  the  melt 
and  the  appropriate  anionic  units.  Therefore,  lead  silicate 
crystalline  phases  have  been  used  wherever  possible.  The 
only  exception  is  for  the  Si04  4  ~  anion,  for  which  there  is 
no  lead  silicate  crystalline  analog. 

The  highest  possible  symmetry  for  an  Si04  4  '  anion  is 
tetrahedral.  From  the  character  table  for  tetrahedral 
molecules  (e.g.  Harris  and  Bertolucci,  1978)  it  can  be  seen 
that  a  tetrahedral  Si04  4  '  ion  has  four  normal  modes  of 
vibration  (A,,  E,  and  2F  2 )  .  Those  normal  modes  are  v, ,  a 
totally  symmetric  Si-0  stretching  mode  (A,),  v2r  a  doubly 
degenerate  O-Si-O  bending  mode  (E) ,  v3f  a  triply  degenerate 
antisymmetric  Si-0  stretching  mode  (F 2 )  ,  and  \/4,  a  triply 
degenerate  O-Si-O  bending  mode  (F2)  (Nakamoto,  1978).  All 
four  of  the  normal  modes  are  Raman  active,  but  only  \/3  and 
\/4  are  IR  active.  The  relative  frequencies  of  the  normal 
modes  generally  follow  the  rules  v3  greater  than  v ,  and  v4 
greater  than  v2  (Nakamoto,  1978).  On  the  basis  of  Raman 
depolarization  ratios,  Furukawa  et  al .  (1978)  conclude  that 
the  \/,  vibration  of  Si044'  occurs  at  approximately  830  cm'1. 
The  IR  spectra  of  olivine  structure  minerals  can  be  used  to 
determine  the  positions  of  the  v3  and  v4  bands.  The  Si04 
tetrahedra  in  Co2Si04  (Morimoto  et  al . ,  1974)  are  less 
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distorted  than  for  other  olivines  (Birle  et  al .  ,  1968),  and 
this  is  illustrated  by  the  smaller  site  splitting  observed 
in  the  IR  bands  of  Co2Si04  than  for  forsterite,  fayalite, 
and  Ni 2Si04-f orster ite  solid  solutions  (Tarte,  1963). 
Additionally,  the  IR  bands  in  olivines  have  been  observed  to 
shift  to  lower  frequencies  with  increasing  cation  mass  and 
radius  (Tarte,  1963;  Burns  and  Huggins,  1972).  For  those 
reasons,  the  IR  spectrum  of  Co2Si04  (Tarte,  1963)  was  used 
to  determine  the  positions  of  the  v3  and  v 4  bands  of  the 
Si044"  anion.  In  Co2Si04  those  bands  fall  at  880  and  480 
cm-1  ,  respectively.  The  \/2  vibration  is  expected  to  fall 
below  480  cm-1  ,  based  on  the  relative  frequency  rules  noted 
above.  Thus,  in  the  400-1200  cm-1  portion  of  the  IR 
spectrum,  a  tetrahedral  Si044*  anion  will  only  contribute 
two  bands  at  approximately  880  and  480  cm-1  . 

If  the  Si044"  anion  is  distorted  from  tetrahedral 
symmetry,  the  degeneracies  noted  above  will  be  split  and  all 
of  the  normal  modes  of  vibration  will  become  IR  active.  The 
IR  spectrum  of  Si044-  anion  will  then  resemble  those  of  the 
olivine  minerals. 

In  the  case  of  the  dimer  (Si2076-)  there  are  two 
crystalline  phases  which  have  been  shown  to  contain  Si20,6' 
units.  Billhardt  (1969),  Petter  et  al .  (1971),  and  Smart  and 
Glasser  (1978)  have  shown  that  Pb3Si207  contains 
pyrosilicate  structural  units.  Similarly,  Gotz  et  al .  (1975b 
and  1976)  and  Smart  and  Glasser  (1978)  have  shown  that  the 
phase  L-Pb2Si04  (using  the  nomenclature  of  Smart  and 
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Glasser,  1978)  also  has  pyrosilicate  structural  units. 
Because  crystallographic  data  (Petter  et  al .  ,  1971)  are 
available  for  Pb3Si207,  but  not  for  L-Pb2Si04,  the  former 
was  used  as  the  crystalline  analog  of  the  dimer.  Infrared 
spectral  data  were  taken  from  Billhardt  (1969). 

The  dimer  IR  spectrum  is  of  particular  interest  because 
the  dimer  possesses  a  non-linear  Si-O-Si  bridge.  The 
spectrum,  consequently  includes  the  vibrational  modes 
characteristic  of  that  bridge.  Normal  coordinate  analyses  of 
other  X2Y 7  (Y3X-Y-XY3  type)  ions  (Brown  and  Ross,  1972)  have 
shown  that  there  are  three  normal  modes  attributable  to  the 
X-Y-X  bridge.  They  are  a  symmetrical  X-Y-X  stretch,  an 
ant i symmet r ical  X-Y-X  stretch,  and  an  in-plane  X-Y-X  bending 
deformation.  Wing  and  Callahan  (1969)  have  shown  that  the 
ant isymmetr ical  stretch  lies  to  high  frequency  of  the 
symmetrical  stretch  and  that  the  separation  is  proportional 
to  the  X-Y-X  bridging  angle.  Billhardt  (1969)  identifies  a 
"characteristic"  pyrosilicate  band  at  679  cm"1  in  the  IR 
spectrum  of  Pb3Si207,  which  by  analogy  to  Na4(P207)  (Bues  et 
al . ,  1963)  and  C1  207  (Roziere  et  al . ,  1973),  is  most  likely 
due  to  the  Si-O-Si  symmetrical  stretch.  Billhardt  (1969) 
also  noted  a  strong  negative  linear  correlation  between 
metal  cation  size  and  the  frequency  of  that  vibrational 
mode.  Similarly,  the  band  at  822  cm'1  may  be  assigned  to  the 
Si-O-Si  antisymmetr ical  stretching  mode.  The  normal 
coordinate  analyses  and  assignments  of  Brown  and  Ross  (1972) 
show  that  the  bending  vibration  lies  to  low  frequency  of 
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both  the  stretching  modes  (generally  between  100  cm'1  and 
300  cm"1  ).  Thus,  it  is  not  expected  to  appear  in  the 
spectra  reported  here. 

The  trimer  provides  a  problem  because  there  is  no  lead 
silicate  crystalline  phase  containing  Si30108'  structural 
units  for  which  crystallographic  data  are  available. 

However,  Smart  and  Glasser  (1978)  found  that  hexagonal 
PbSi03  contained  Si30108‘  units.  In  addition  Furukawa  et  al . 
(1978  and  1979)  found  from  vibrational  spectroscopy  that 
hexagonal  PbSi03  is  not  a  chain  silicate  such  as  alamosite, 
which  supports  the  chromatographic  determination  that  the 
structural  unit  in  hexagonal  PbSi03  is  the  trimer.  Thus, 
hexagonal  PbSi03  was  used  as  the  crystalline  analog  of  the 
trimer.  Infrared  spectral  data  were  taken  from  Furukawa  et 
al .  (1978  and  1979). 

The  tetramer  presents  an  additional  problem  because 
there  are  four  structural  isomers.  Those  are:  chain, 
branched  chain,  three-membered  ring  with  one  branch,  and 
f our-membered  ring.  Smart  and  Glasser  (1978)  found  that 
silicate  glasses  along  the  join  Pb0-Si02  contain  tetramers 
in  the  form  of  f our-membered  rings  (Si40128').  As  far  as  a 
crystalline  analog  is  concerned,  Gotz  et  al .  (1975a  and 
1976),  Smart  and  Glasser  (1978),  and  Furukawa  et  al .  (1979) 
all  confirm  that  H-Pb2Si04  (using  the  nomenclature  of  Smart 
and  Glasser,  1978)  contains  Si40128'  rings.  Therefore, 
H-Pb2Si04  was  used  as  the  crystalline  analog  of  the 
tetramer.  IR  spectral  data  were  taken  from  Furukawa  et  al . 
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(  1979)  . 

In  a  melt,  the  anionic  units  are  not  held  rigidly  as  in 
a  crystal.  The  anions  will  presumably  undergo  continuous 
deformation  as  a  consequence  of  the  dynamic  melt 
environment.  Thus,  upon  quenching  to  the  glass  state,  the 
individual  anions  of  a  given  species  will  exhibit  a  range  of 
shapes  clustered  around  some  average  structure.  Such  small 
differences  will  lead  to  broader  IR  bands  than  is  the  case 
for  the  crystalline  analogs.  Consequently,  the  crystal 
analog  spectra  described  above  were  broadened  and  smoothed 
by  fitting  Gaussian  curves  to  the  bands.  The  method  is 
described  below  in  some  detail.  The  "quasi-melt"  spectra  for 
the  anionic  units  up  to  Si40128'  are  shown  in  fig.  17. 

Larger  anionic  units  were  not  considered  due  to  the  large 
number  of  structural  isomers  possible  for  each  case  and 
because  there  currently  are  no  data  which  allow  a  particular 
structure  (or  structures)  to  be  chosen  for  a  given  polymer. 
However,  qualitatively  speaking,  the  IR  spectra  of  such 
anions  will  tend  to  approach  that  of  vitreous  Si02  (fig.  18) 
with  increasing  anionic  size. 

The  smoothed  "quasi-melt"  spectra  in  fig.  17,  when 
taken  in  conjunction  with  the  discussion  above,  show  that 
the  IR  spectra  of  silicate  anions  may  be  divided  into  three 
principal  regions.  Those  are:  the  O-Si-O  bending  region 
below  about  600  cm-1  ,  the  Si-O-Si  stretching  region  between 
600  and  830  cm"1  ,  and  the  Si-0  stretching  region  above  830 
cm-1  .  Further,  several  trends  are  apparent.  With  increasing 
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Figure  17.  Smoothed  infrared  spectra  of  the  crystalline 

analogs  of  the  monomer  (M) ,  dimer  (D) ,  trimer  (T) ,  and 
tetramer  (Te ) . 
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Figure  1 
and 


8.  Infrared  spectrum  of  vitreous  Si02  from  Ishikawa 
Akagi  ( 1 978 ) . 
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anionic  size,  the  median  of  the  Si-0  stretching  band  shifts 
to  higher  frequency  and  the  Si-O-Si  stretching  mode  shifts 
to  lower  frequency.  These  features  allow  the  qualitative 
interpretation  of  the  IR  spectra  of  silicate  melts  and 
glasses . 

It  is  also  possible  to  use  the  smoothed  crystalline 
analog  spectra  (fig.  17)  to  qualitatively  deconvolute  the 
glass  spectra.  However,  by  simply  examining  the  spectrum,  it 
is  possible  to  tell  with  reasonable  certainty  if  there  are 
contributions  attributable  to  a  particular  anionic  species. 
In  order  to  do  more  than  that  by  deconvolution  methods, 
there  must  either  be  data  about  the  relative  intensities  of 
the  vibrational  modes  for  different  anionic  species,  or 
assumptions  must  be  made  about  those  intensities.  This 
presupposes  that  the  possible  anionic  species  are  known. 
Without  such  data,  deconvolution  relies  on  the  assumption 
that  the  relative  intensities  of  vibrational  modes  in  the 
same  class  (Si-0  stretches  etc.)  are  the  same  for  different 
anionic  species.  Deconvolution  schemes  currently  in  use  for 
silicate  glass  spectra  (e.g.  Mysen  et  a/.,  1982a)  depend  on 
a  determination  of  the  number  of  bands  in  a  particular 
composite  band.  The  bandwidths,  intensities,  and  positions 
are  then  varied  in  order  to  minimize  the  difference  between 
the  real  and  synthetic  spectra.  While  such  methods  have  met 
with  some  success  for  Raman  spectra  of  silicate  glasses 
(Mysen  et  al . ,  1980d  and  1982a),  the  large  number  of  IR 
active  bands  for  even  simple  silicate  anions  (fig.  17) 
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significantly  complicates  the  problem.  Additionally,  the 
results  of  spectral  deconvolution  by  iteratively  varying 
parameters  (band  width,  intensity,  and  band  position)  and 
evaluating  the  goodness  of  fit  by  minimizing  residuals  are 
strongly  dependent  on  the  order  in  which  the  variables  are 
changed.  Because  of  these  difficulties,  the  spectra 
presented  in  this  paper  were  not  deconvoluted .  The 
interpretations  presented  here  are  based  on  inspection  and 
comparison  of  the  spectra  with  the  crystalline  analog 
spectra . 

Series  one 

As  was  previously  noted,  the  series  one  spectra  are 
very  similar  and  are  treated  as  one.  The  480  cm-1  band 
indicates  that  O-Si-O  bending  modes  are  present.  There  is  no 
way  of  telling  whether  the  oxygens  involved  in  those  modes 
are  non-bridging  oxygens  (NBO)  or  are  bridging  oxygens  (BO) 
between  two  silicons.  The  880  cm-1  band,  on  the  other  hand, 
provides  more  information.  The  principal  maximum  in  that 
band  is  most  likely  due  to  the  ant isymmetr ical  Si-0 
stretching  mode  of  Si044-  (\/3).  However,  the  presence  of 
both  low  and  high  frequency  shoulders  in  that  band  argues 
either  for  contributions  from  larger  silicate  anions,  or  for 
a  deformed  Si044'  anion.  The  pronounced  shoulder  at 
approximately  650  cm-1  on  the  low  side  of  the  Si-0 
stretching  envelope  is  in  the  Si-O-Si  stretching  region  and 
thus  indicates  the  presence  of  silicate  anionic  units  larger 
than  the  monomer.  The  most  likely  candidate  is  the  dimer 
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(Si2076')  which  has  bands  at  680,  820  and  910  cm-1  (fig. 

17).  The  spectra  do  not  show  other  features  which  may  be 
confidently  assigned  to  particular  anionic  units,  but  that 
does  not  preclude  the  presence  of  small  proportions  of  such 
units . 

Smart  and  Glasser  (1978)  report,  based  on 
chromatographic  resolution  of  tr imethylsilyl  derivatives  of 
glass,  that  Pb2Si04  glasses,  prepared  by  short  duration 
fusion  and  rapid  quenching,  contain  21%  monomer,  30%  dimer, 
17%  trimer,  7%  ring  tetramer,  2-3%  ring  hexamer,  and  22-23% 
polysilicate  (larger  than  hexamer)  anions.  In  order  to 
estimate  the  appearance  of  the  IR  spectrum  of  such  a  melt,  a 
synthetic  spectrum  was  constructed. 

The  synthetic  spectrum  was  constructed  by  adding  the 
digitized  (%  transmission)  spectra  of  the  crystalline 
analogs  together  in  the  proportions  reported  by  Smart  and 
Glasser  (1978).  Because  the  spectra  of  Pb2Si04  glasses  do 
not  have  any  features  above  1000  cm-1  which  are  attributable 
to  the  vibrations  of  silicate  anions,  only  contributions 
from  the  monomer,  dimer,  trimer,  and  tetramer  were  included. 
The  sum  of  the  contributions  was  then  transformed  to 
absorbance  and  plotted.  This  procedure  assumes  that  the 
maximum  infrared  activity  of  the  strongest  absorption  for 
each  of  the  anionic  species  is  the  same.  The  synthetic 
spectrum  was  smoothed  and  broadened  by  fitting  Gaussian 
bands  to  the  observed  maxima.  The  band  widths  (FWHM)  of  the 
Gaussian  curves  were  assigned  according  to  the  position  of 
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the  band.  Bands  below  600  cm-1  were  assigned  a  band  width  of 
100  cm-1.  Band  widths  of  200  and  160/200  cm'1  were  used  for 
bands  between  600  and  800  cm' 1  and  above  800  cm' 1 , 
respectively.  The  result,  shown  in  fig.  19  along  with  one  of 
the  series  one  spectra  (Pb-1;  Table  11),  is  very  similar  to 
the  series  one  spectrum,  except  that  the  series  one  spectrum 
has  a  more  intense  shoulder  at  approximately  650  cm' 1 .  That 
result,  while  not  proving  that  higher  order  polymers  are 
present  in  Pb2Si04  melts,  lends  support  to  the  results  of 
Smart  and  Glasser  (1978).  The  conclusion  for  the  series  one 
Pb2Si04  melts  is  that  both  monomeric  and  dimeric  anionic 
units  are  present  and  that  small  proportions  of  larger 
anionic  units  may  also  be  present.  That  conclusion,  in 
conjunction  with  mass  balance  considerations,  requires  that 
a  sizeable  fraction  of  the  oxygen  in  Pb2Si04  melts  be 
present  as  the  free  anion  (O2 " ) . 

Series  two,  three,  and  four 

These  results  show  progressive  changes  in  the  IR 
spectra  either  as  a  consequence  of  decreasing  cooling  rate, 
or  as  a  function  of  increasing  isothermal  soaking  time.  The 
sense  of  those  changes  is  the  same  for  all  three  series  of 
experiments.  The  most  polymerized  melt  is  that  from 
experiment  PB-10  (Table  11  and  figs.  15  and  16),  which 
underwent  both  the  slowest  cooling  and  the  longest 
isothermal  soaking.  The  spectra  of  the  other  experiments  are 
intermediate  between  those  of  series  one  and  that  of  PB-10. 
Thus,  a  discussion  of  the  PB-10  results  serves  to  delineate 
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Figure  19.  Synthetic  infrared  spectrum  of  Pb2Si04 
constructed  by  adding  the  infrared  spectra  of 
crystalline  analogs  of  various  silicate  anions  in  the 
proportions  determined  by  Smart  and  Glasser  (1978).  Both 
raw  (a)  and  smoothed  (b)  spectra  are  shown.  The 
smoothing  parameters  are  given  in  the  text.  The  infrared 
spectrum  of  glass  from  experiment  Pb- 1  is  shown  for 
compar i son  (c ) . 


95 


the  entire  series  of  progressive  changes. 

The  IR  spectrum  of  PB-10  includes  identifiable  bands  at 
465,  485,  820,  880,  and  920  cm'1  .  There  are  also  definite 
shoulders  at  525,  650,  and  at  greater  than  950  cm-1  . 
Furthermore,  the  band  width  ( FWHM )  of  the  Si-0  stretching 
envelope  is  much  greater  than  that  of  the  series  one  spectra 
(400  cm'1  versus  270  cm'1  ),  and  the  mid-point  of  that 
envelope  is  shifted  by  20-30  cm-1  to  lower  frequency.  Those 
band  positions  and  changes  (relative  to  series  one)  are 
consistent  with  the  presence  of  both  monomer  and  dimer  as 
well  as  some  larger  anionic  units  (the  shoulder  above  950 
cm'1  implies  at  least  the  tetramer).  The  650  cm-1  shoulder 
was  also  noted  to  have  increased  intensity  relative  to 
series  one,  which  argues  for  an  increasing  contribution  by 
larger  anionic  units. 

Relative  to  the  series  one  experiments,  the  spectrum  of 
PB-10  shows  decreasing  contributions  from  the  monomer 
coupled  with  increasing  contributions  from  dimeric  and 
larger  silicate  anions.  These  results  show  that  both  types 
of  thermal  treatment  (slow  cooling  and  isothermal  soaking) 
lead  to  increasing  polymerization  of  Pb2Si04  melts  (i.e. 
NBO/T  decreases).  Mass  balance  considerations  require,  as  a 
corollary  to  that  conclusion,  that  such  thermal  treatments 
also  increase  the  proportion  of  free  oxygen  anions  (02‘)  in 
Pb2Si04  melts.  It  is  possible  to  calculate  the  proportion  of 
O2"  anion  in  the  glass  from  the  mass  balance  using  the 
anionic  proportions  reported  by  Smart  and  Glasser  (1978).  If 
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the  polysilicate  reported  by  Smart  and  Glasser  (1978)  is 
assumed  to  be  completely  polymerized  (Si02),  the  calculation 
yields  15%  O2  ~ . 

At  this  point  it  is  necessary  to  consider  the  IR 
spectrum  of  vitreous  fayalite  and  the  interpretations 
reported  by  Kusabiraki  and  Shiraishi  (1981).  The  IR  spectrum 
of  vitreous  fayalite  has  strong  bands  centered  at  508  and 
933  cm"1  ,  as  well  as  a  weak  band  at  695  cm-1  . 

Additionally,  the  933  cm-1  band  has  a  pronounced  high 
frequency  shoulder  above  1000  cm-1  .  Kusabiraki  and 
Shiraishi  (1981)  interpret  that  spectrum  as  being  totally 
due  to  the  vibrations  of  Si044'  monomers.  They  assign  the 
695  cm-1  band  to  the  symmetric  Si-0  stretching  mode  of 
tetrahedral  Si04  4  “  (v, )  and  completely  neglect  the  high 
frequency  shoulder  in  the  933  cm'1  band.  Their  assignment  of 
the  695  cm'1  band  is  incorrect.  From  the  discussion  of  the 
crystal  spectra  above,  it  is  clear  that  the  695  cm"1  band  is 
within  the  Si-O-Si  stretching  region,  not  the  Si-0 
stretching  region.  Furukawa  et  al  (1978)  report  the  v/,  mode 
of  tetrahedral  Si044'  at  approximately  830  cm'1  .  The 
position  of  the  695  cm'1  band  is  in  very  close  accord  with 
that  predicted  from  the  ionic  radius  of  Fe2+  and  the 
correlation  noted  by  Billhardt  (1969)  between  metal  cation 
radius  and  the  position  of  the  Si-O-Si  stretching  band  for 
various  pyros i 1 icates .  Similarly,  the  high  frequency 
shoulder  on  the  933  cm-1  band  is  indicative  of  silicate 
anions  at  least  as  large  as  the  tetramer.  The  IR  spectrum  of 
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vitreous  fayalite  is  consistent  with  the  interpretation  that 
vitreous  fayalite  contains  silicate  anions  larger  than 
Si044".  However,  the  low  intensity  of  the  Si-O-Si  stretching 
band  (695  cm' '  )  implies  that  the  dominant  anionic  species 
is  the  monomer. 

E.  Conclusions 

This  study  has  shown  that  Pb2Si04  melts  contain 
silicate  anions  larger  than  Si044-  monomers.  The  anionic 
species  that  are  present  and  their  relative  proportions  have 
been  shown  to  depend  on  the  thermal  history  of  the  melt. 
These  observations  require  that  some  of  the  oxygen  in  the 
melt  not  be  bound  to  Si  (i.e.  oxygen  anions  exist  in  the 
melt).  That  ionic  oxygen  exists  in  Pb2Si04  melts  and  that 
those  melts  are  structurally  similar  to  Fe2Si04  melts 
(Kusabiraki  and  Shiraishi,  1981)  can  be  taken  as  implying 
that  ionic  oxygen  may  exist  in  all  orthosilicate  melts.  The 
results  also  show  that  NBO/T  for  Pb2Si04  melts  is  less  than 
4  and  that  NBO/T  depends  on  the  thermal  history  of  the  melt. 


- 1 


- 


V.  The  Oxidation  State  of  Iron  in  Silicate  Melts:  New  Data 


for  Basaltic  Melts 


A.  Introduction 

The  relationship  between  the  oxidation  state  of  iron  in 
natural  rock  melts  and  the  bulk  composition,  temperature  and 
oxygen  fugacity  has  been  investigated  by  a  number  of  workers 
in  recent  years  (Kennedy,  1948;  Fudali  1965;  Thornber  et 
al . ,  1980;  and  Sack  et  al .  ,  1980).  Other  workers  including 
Paul  and  Douglas  (1965),  Mysen  and  Virgo  (1978),  and 
Dickensen  and  Hess  (1981)  have  investigated  the  same 
relationship  for  iron  in  simple  synthetic  silicate  melts. 

The  volume  of  work  on  synthetic  melts  far  surpasses  that  for 
natural  melts.  The  result  is  that,  except  for  the  equation 
presented  by  Sack  et  a/.  (1980),  there  is  no  way  to  predict 
the  iron  redox  equilibria  in  natural  melts  for  given 
conditions.  Furthermore,  much  of  the  work  on  natural  melts 
was  done  at  low  oxygen  fugacities  (10'6-10"s  bars).  Only 
Kennedy  reports  data  determined  at  high  oxygen  fugacity 
(0.21  bar).  In  addition,  only  Kennedy  (1948)  has  done 
experiments  on  a  single  composition  at  more  than  two 
temperatures.  Consequently,  while  the  relationship  presented 
by  Sack  et  al .  (1980)  has  a  diverse  compositional  base  (143 
data  sets  extending  from  nephelinitic  to  granitic 
compositions),  the  range  of  oxygen  fugacities  is  quite 
limited  and  the  data  were  determined  at  only  3  or  4 
different  temperatures.  Therefore,  this  study  was  undertaken 
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to  expand  upon  the  data  base  used  by  Sack  et  al .  (1980)  by 
providing  new  data  for  basaltic  melts  at  high  oxygen 
fugacity  and  over  a  range  of  temperatures. 

The  results  of  27  experiments  on  three  basaltic 
compositions  are  presented  in  this  paper.  The  bulk  of  the 
experiments  were  done  in  air  and  span  the  temperature  range 
1 250- 1 437°C .  The  new  data  were  combined  with  those  of 
Kennedy  (1948),  Fudali  (1965),  Thornber  et  al .  (1980),  and 
Sack  et  al .  (1980)  and  an  equation,  relating  the 
ferric/ferrous  ratio  of  the  melt  to  temperature,  oxygen 
fugacity,  and  bulk  composition,  was  determined  by  multiple 
regression  analysis  of  the  data.  The  results  presented  in 
this  paper  expand  the  previous  data  base  in  terms  of 
temperature  and  oxygen  fugacity.  This  increases  the  validity 
of  extrapolating  the  predictive  equation  outside  the 
temperature  and  oxygen  fugacity  range  of  the  data  set.  In 
addition,  including  27  new  data  sets  for  basalts  improves 
the  accuracy  of  the  method  for  basaltic  liquids. 

B.  Experimental  and  Analytical 

The  experiments  were  done  on  three  different  "basaltic" 
melts.  Those  compositions  are  C-11  (a  nephelini te ) ,  KR- 1 3 
(an  alkali  basalt),  and  1921  Kilauea  tholeiite.  The  sample 
of  C-11  was  provided  by  Dr.  T.  Fujii,  the  sample  of  KR- 1 3 
was  provided  by  Dr.  C.M.  Scarfe,  and  the  sample  of  1921 
Kilauea  tholeiite  was  provided  by  Dr.  I  Kushiro.  The 
starting  materials  were  initially  ground  to  powder  in  a 
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tungsten  carbide  swing  mill.  The  powders  were  then  fused  in 
covered  Pt  crucibles  in  air  at  1350°C  for  at  least  2  hours. 
The  experiments  were  done  at  both  low  ( logfO 2 =-8 . 0 )  and  high 
( logf02 =-0 . 68 )  oxygen  fugacities  in  an  electrically  heated 
vertical  tube  furnace  (Deltech  DT-31VT)  equipped  with  a 
gas-tight  Al  203  muffle  tube.  Oxygen  fugacities  were  set 
either  by  leaving  the  furnace  open  to  air,  or  by  passing 
mixtures  of  CO  and  C02  through  the  furnace.  The  ratio  of  CO 
to  CO 2  used  in  the  low  oxygen  fugacity  experiments  was  set 
with  an  oxygen  fugacity  probe  of  the  type  used  by  Arculus 
and  Delano  (1981).  The  low  oxygen  fugacities  are  estimated 
to  be  accurate  to  within  ±0.25  log  units  of  f 02 .  The 
temperatures  during  the  experiments  were  monitored  with  a 
Pt-Pt 8 ,Rh, 3  thermocouple  which  was  calibrated  against  the 
melting  point  of  Au  (1064. 5°C). 

All  but  three  of  the  experiments  were  done  by  the  Pt 
wire  loop  method  (Donaldson,  1975).  The  other  three 
experiments  were  done  as  batch  fusions  in  Pt  crucibles.  The 
samples  were  placed  on  the  Pt  loops  as  slurries  with 
acetone.  After  the  acetone  evaporated  from  the  slurry,  the 
sample  was  lightly  fused  onto  the  loop  with  an  oxy-acetylene 
torch.  The  low  oxygen  fugacity  experiments  were  done  in  the 
same  manner.  Iron  loss  to  the  Pt  loops  was  not  thought  to  be 
a  problem  as  the  samples  were  analyzed  after  the 
experiments.  Each  experiment  included  several  samples  of  the 
same  basalt.  Each  sample  was  split  into  several  pieces.  One 
piece  was  used  for  bulk  chemical  analysis  and  the  other 
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Figure  20.  Calculated  time  required  for  equilibration  in 

iron  oxidation  state  experiments  for  a  spherical  50  mg. 
sample  with  a  density  of  2.75  g/cm3 .  The  method  is 
explained  in  the  text. 
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pieces  were  analyzed  for  FeO.  In  most  cases  the  samples 
analyzed  for  FeO  were  also  analyzed  for  total  iron  by  atomic 
absorption  using  BCR- 1  (Flanagan,  1973)  as  a  standard. 

The  bulk  chemical  analyses  were  done  on  the  University 
of  Alberta  ARL-SEMQ  electron  microprobe.  Suitable  standards 
were  used  to  allow  data  reduction  by  the  method  of  Bence  and 
Albee  (1968)  using  the  alpha  factors  of  Albee  and  Ray 
(1970).  The  FeO  analyses  were  done  by  the  method  of  Wilson 
(1960).  Microsyringes  were  used  for  the  titrations  in  the 
FeO  analyses,  thus  allowing  a  precision  of  ±0.01  wt .  %  to  be 
achieved  for  most  of  the  analyses.  The  solutions  used  were 
frequently  standardized  against  one  another  to  assure  that 
the  ferrous  ammonium  sulphate  solution  had  not  become 
oxidized.  The  accuracy  of  the  FeO  analyses  was  checked  by 
analyzing  25  aliquots  of  BCR- 1  (Flanagan,  1973).  The  result 
of  those  analyses  was  8.66±0.21  wt  %  as  compared  to  the 
reported  value  of  8.80  wt .  %  (Flanagan,  1973).  All  of  the 
experimental  products  produced  in  this  study  were  totally 
glassy  as  determined  by  optical  examination  of  the  crushed 
samples . 

If  experiments  of  the  type  done  in  this  study  are  to 
have  any  meaning,  they  must  be  equilibrium  experiments. 
Thornber  et  al .  (1980)  argue,  on  the  basis  of  a  time  study, 
that  only  5-7  hours  are  required  to  reach  equilibrium  in 
wire  loop  experiments  with  samples  of  up  to  100  mg.  However, 
their  replicate  experiments  all  fail  to  produce  the  same 
Fe3+/Fe2+  value  for  given  conditions,  but  different 
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experimental  durations.  The  rate  of  approach  to  equilibrium 
in  iron  oxidation  experiments  is  controlled  by  the  diffusion 
of  oxygen  either  into  or  out  of  the  sample  (see  chapter  6). 
Thus,  measured  oxygen  dif fusivities  in  silicate  melts  may  be 
used  to  calculate  the  time  required  for  equilibration  in 
such  experiments.  Most  of  the  experiments  reported  here  were 
done  on  KR- 1 3  alkali  basalt.  The  zero  pressure  oxygen 
diffusivity  for  that  composition  was  determined  in  chapter  3 
(Table  9)  and  is  approximately  10‘6cm2/sec  at  1400°C. 
Assuming  an  activation  energy  of  approximately  70  kcal 
(Table  8)  the  diffusion  coefficient  is  approximately 
3x 1 0 “ 7cm2/sec  at  1300°C.  That  diffusion  coefficient  can  be 
used  in  conjunction  with  eqn .  16: 

r  =  (Dt)  1/2  (16) 

to  calculate  an  equilibration  time  for  wire  loop  experiments 
by  letting  r  equal  the  radius  of  the  sample,  D  equals  the 
oxygen  diffusion  coefficient,  and  t  equal  the  duration  of 
the  experiment.  That  was  done,  assuming  a  50  mg.  sample  mass 
and  a  melt  density  of  2.75  g/cm3.  The  results  are  presented 
in  fig.  20  in  which  equilibration  time  is  plotted  as  a 
function  of  experimental  temperature.  Fig.  20  was  used  to 
determine  the  experimental  durations  used  in  this  study. 

C.  Results  and  discussion 

The  experimental  results  are  given  in  Table  12.  Sack  et 
a/.  (1980)  used  the  technique  of  multiple  regression  to 
derive  an  equation  of  the  form: 
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Table  12.  Analyses  of  experimental  products.  Analyses  were 
done  by  electron  microprobe.  FeO  was  determined  by  the 
method  of  Wilson  (1960).  Fe203  was  determined  by 
difference  from  total  iron  determined  either  by  electron 
microprobe  or  by  atomic  absorption. 


Exp.  # 

15 

16 

17 

18 

19 

20 

f  o2 

air 

air 

air 

air 

air 

air 

T°C 

1262 

1278 

1300 

1300 

1300 

1300 

t  (hr ) 

91 

70 

23 

82 

48 

144 

Si02 

47.20 

47.90 

49.  14 

46.20 

46.90 

47.60 

TiO  2 

2.50 

2.40 

2.14 

2.80 

2.70 

2.50 

A1203 

15.20 

15.80 

16.14 

14.90 

15.30 

15.60 

MgO 

5.70 

6.00 

5.99 

6.00 

5.80 

5.90 

CaO 

9.60 

10.00 

9.60 

9.40 

9.50 

9.50 

FeO 

1  .49 

1 .76 

2.05 

1.84 

1  .92 

2.11 

Fe  203 

10.24 

10.16 

9.22 

12.07 

10.54 

10.70 

MnO 

n  .a . 

n  .a . 

n .  d. 

n .  a . 

n  .a . 

n .  a . 

Na  20 

3.60 

3.40 

4.03 

3.10 

3.40 

3.10 

K20 

1.11 

1  .02 

1  .23 

0.87 

0.95 

0.61 

p2o5 

n  .a . 

n  .a . 

0.41 

n  .a . 

n  .a . 

n  .a . 

Total 

97.63 

99.44 

99.94 

98.37 

98.12 

98.57 

Exp.  # 

21 

22 

25 

26 

27 

28 

f  02 

air 

air 

air 

air 

air 

air 

T°C 

1305 

1320 

132  1 

1340 

1350 

1360 

t  (hr ) 

81 

17 

64 

48 

47 

24 

Si02 

47.00 

47.80 

47.20 

50.40 

47.70 

51.42 

TiO  2 

2.50 

2.28 

2.60 

2.11 

2.40 

2.32 

A1203 

14.80 

16.00 

15.60 

15.64 

15.50 

16.05 

MgO 

5.80 

5.82 

5.80 

6.05 

5.90 

5.89 

CaO 

9.30 

9.54 

9.60 

9.82 

9.40 

9.57 

FeO 

2.01 

1  .91 

2.08 

2.37 

2.40 

2.62 

Fe  203 

10.88 

11.91 

10.45 

8.68 

9.94 

9.06 

MnO 

n .  a . 

0.14 

n  .a . 

0.19 

n  .a . 

n  ,d. 

Na  20 

3.30 

3.68 

3.30 

3.73 

3.40 

3.53 

K20 

0.78 

1.18 

0.84 

0.96 

0.86 

0.97 

p2o5 

n .  a . 

0.36 

n .  a . 

0.  14 

n .  a . 

0.21 

Total 

97.47 

100.63 

98.59 

100.09 

97.50 

101.64 

.  -  not  detected. 
.  -  not  analyzed. 


n .  d 
n .  a 
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Table  12.  Analyses  of  experimental  products  t(cont.). 


Exp.  # 

29 

30 

31 

32 

33 

34 

fo2 

air 

air 

air 

air 

air 

air 

T°C 

1380 

1380 

1400 

1400 

1420 

1420 

t  (hr ) 

18 

27 

20 

23 

23 

26 

Si02 

46.32 

47.10 

49.24 

47.70 

47.30 

47.20 

TiO  2 

3.57 

2.50 

2.26 

2.50 

3.20 

2.70 

Al  2  0  3 

14.78 

15.20 

16.27 

15.30 

15.46 

15.50 

MgO 

5.95 

5.90 

6. 1  1 

5.80 

6.01 

5.90 

CaO 

8.82 

10.70 

9.67 

9.50 

9.17 

9.30 

FeO 

3.28 

2.93 

3.51 

2.66 

2.88 

3.16 

Fe  203 

13.18 

10.34 

7.63 

9.69 

13.80 

10.57 

MnO 

n.d. 

n .  a . 

0.15 

n .  a . 

0.13 

n .  a . 

Na  20 

3.32 

3.20 

3.65 

3.50 

3.10 

3.00 

K20 

0.90 

0.95 

1.10 

1  .05 

0.75 

0.85 

P2Os 

0.23 

n .  a . 

0.17 

n .  a . 

0.04 

n  .a . 

Total 

100.35 

98.82 

99.67 

97.70 

101.84 

98.  18 

Exp.  # 

35 

39 

41 

43 

93 

94 

f  02 

air 

-8.  18 

-8.30 

-8.15 

air 

air 

T°C 

1438 

1250 

1300 

1350 

1350 

1400 

t  (hr ) 

21 

91 

61 

45 

45 

24 

Si02 

48.70 

47.60 

49.40 

49.60 

41.36 

39.63 

Ti02 

2.57 

2.80 

2.50 

2.60 

2.93 

2.97 

Al  203 

15.81 

15.60 

16.70 

16.50 

12.37 

12.59 

MgO 

5.93 

6.20 

6.10 

6.40 

9.38 

9.48 

CaO 

9.  13 

9.90 

10.50 

10.40 

13.04 

12.80 

FeO 

2.41 

10.63 

9.00 

8.24 

2.19 

2.65 

Fe  2O3 

10.65 

0.96 

0.79 

0.66 

14.57 

13.75 

MnO 

n.d. 

n .  a . 

n .  a . 

n  .a . 

0.25 

0.27 

Na  20 

4.75 

3.00 

2.90 

2.30 

3.34 

3.63 

K20 

1.10 

0.95 

0.94 

0.79 

0.76 

0.91 

p2o5 

0.02 

n  .a . 

n  .a . 

n  .a . 

0.72 

0.82 

Total 

100.14 

97.69 

98.83 

97.49 

100.91 

99.50 

n .  a 
n .  d 


not  analyzed 
not  detected 
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Table  12.  Analyses  of  experimental  results  (cont.) 


Exp.  # 

f  02 

T°C 

t  (hr ) 

Si02 

Ti02 

A1203 

MgO 

CaO 

FeO 

F  e  2  0  3 

MnO 

Na  20 

K20 

P20s 

Total 


C-1  1 
air 
1350 
24 
37.77 
2.75 
12.00 

9.29 
12.59 

2.08 

13.46 

0.22 

4.29 
2.17 
2.37 

98.99 


KR-  1  3 
air 
1400 
43 

49.63 
2.43 

14.63 
5.91 
9.38 
2.26 

10.10 
0.21 
4.06 
1  .25 
0.50 
100.36 


1921 

air 

1350 

23 

48.08 

2.97 

12.11 

11.30 

9.44 

2.44 
10.05 

0.18 

2.25 

0.51 

0.28 

99.68 


n.a.  -  not  analyzed, 
n.d.  -  not  detected. 

Experiments  15-43  were  done  on  KR- 1 3 ,  experiments  93  and  94 
were  done  on  C— 11. 


7r?{X(Fe  203  )/X(FeO)  }  =  a  Inf  (02  )  +b/T+c+  £  d  (  i  )X  (  i  )  (17) 

i 

where  X ( i )  is  the  mole  fraction  of  oxide  component  i  in  the 
melt  and  a,  b,  c,  and  d  are  coefficients  determined  from  the 
multiple  regression  analysis.  In  their  analysis  Sack  et  al . 
(1980)  excluded  the  oxides  Ti02,  MnO,  Cr203,  and  P205  on  the 
basis  that  those  oxides  can  have  very  low  concentrations  in 
silicate  melts  and  are  frequently  not  reported  in  analyses. 
The  oxides  Ti02 ,  MnO,  and  P205  were  included  in  the  analysis 
reported  here  because  examination  of  the  data  set  showed, 
that  while  the  proportions  of  the  oxides  are  generally  low, 
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there  was  considerable  variation  in  their  concentrations. 
Cr203  was  excluded  from  the  analysis  for  the  same  reasons 
given  by  Sack  et  al .  (1980). 

The  multiple  regression  analyses  were  done  using  the 
program  MULTR  (Davis,  1973)  (translated  into  basic).  The 
program  was  modified  to  allow  the  multiple  regression 
coefficients  to  be  standardized.  The  regression  coefficients 
were  standardized  using  equation  18 

b  =  d{(std.  dev.  I)/(std.  dev  B)}  (18) 

where  B  is  the  dependent  variable,  I  is  an  independent 
variable,  d  is  the  multiple  regression  coefficient  of 
variable  I,  and  b  is  the  standardized  multiple  regression 
coefficient  of  variable  I.  Standardizing  the  regression 
coefficients  allows  direct  comparison  of  the  regression 
coefficients  in  order  to  determine  which  of  them  accounts 
for  more  of  the  variance  of  the  dependent  variable  (Davis, 
1973)  . 

The  multiple  regressions  were  done  incrementally  as 
different  groups  of  analyses  were  added  to  the  data  set.  The 
grouping  of  analyses  for  the  multiple  regressions  was:  (A) 
the  27  analyses  presented  in  this  study  (Table  12);  (B)  the 
57  analyses  presented  by  Sack  et  al .  (1980);  (C)  the 
combination  of  A  and  B  above;  (D)  the  results  of  Thornber  et 
a/.  (1980)  plus  C  above;  and  (E)  the  results  of  Kennedy 
(1948)  and  Fudali  (1965)  plus  D  above  (150  analyses).  The 
results  of  the  multiple  regression  analyses  are  presented  in 
Table  13.  The  largest  data  set  (E)  is  not  significantly 
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Table  13.  Results  of  multiple  regression  analyses.  The 

columns  are  labelled  according  to  the  cases  defined  in 
the  text. 


Variable 

A 

B 

C 

D 

E 

Info  2 

0.207 

0.297 

0.231 

0.230 

0.228 

1/T 

6344 

25124 

12415 

1  1535 

11676 

Const . 

-103.24 

-34.93 

-112.46 

-56.35 

0.16 

Si02 

100.21 

22.44 

106.74 

50.88 

-5.78 

TiO  2 

77.94 

20.37 

103.30 

49.00 

-5.65 

A1203 

97.94 

10.34 

94.25 

42.75 

-13.35 

MgO 

99.61 

18.09 

102.26 

46.43 

-10.29 

CaO 

93.08 

25.93 

109.11 

53.  18 

-3.59 

FeO 

112.22 

18.45 

107.92 

52.30 

-4.44 

MnO 

148.34 

60.32 

69.29 

46.67 

-13.00 

Na  20 

119.54 

31.90 

119.13 

57.46 

0.48 

K20 

18.00 

30.53 

115.92 

59.31 

1.87 

p2o5 

199.50 

13.72 

112.58 

56.56 

2.61 

r  2 

0.99 

0.54 

0.97 

0.96 

0.96 

r2  is  the  goodness  < 

of 

fit  from  the  regression  analyses. 

larger  than 

that  analyzed  by  Sack  et  al 

.  ( 1980)  due 

to  the 

rather  more 

stringent 

equilibrium  criteria  applied 

in  this 

study.  Table  14  lists 

the  variables  in 

order  of  the 

relative 

proportion 

of  the  variance  attributable 

to  the  variable  for 

each  of  the 

cases  above. 

The  regression 

coefficients 

for  the  complete  data  set 

(case  E)  compare  favorably  with 

those  determined  by 

Sack  et 

al .  (1980) 

except  that 

the  sign 

of  the 

coef  f ic ient 

for  CaO 

is  negative 

in  this 

study,  while 

it  was 

positive  in 

Sack ' s 

study.  The  coefficients  for  lnf(02),  X(Na20),  X(K20) ,  and 
X(P205)  are  all  positive.  Thus,  an  increase  in  any  of  those 
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components  can  be  expected  to  coincide  with  an  increase  of 
Fe3+  in  the  melt.  However,  it  is  of  interest  to  note  that 
all  of  the  regression  coefficients  are  positive  for  cases  A, 
B,  C,  and  D.  Only  in  case  E  are  there  any  negative 
regression  coefficients.  Standardizing  the  regression 
coefficients  (Table  14)  shows  that  the  fraction  of  the 
variance  of  the  dependent  variable  attributable  to  any  given 
independent  variable  is  dependent  on  the  variance  of  all  the 
variables  in  the  data  set.  In  other  words,  if  an  independent 
variable  varies  greatly  in  the  data  set,  then  it  will  have  a 
proportionately  large  fraction  of  the  variance  of  the 
dependent  variable  attributed  to  it.  Therefore,  multiple 
regression  analysis,  while  providing  a  useful  way  of 
determining  a  predictive  model,  provides  little  or  no 
insight  into  the  real  relationships  between  the  dependent 
and  independent  variables. 

D.  Applications 

Sack  et  a/.  (1980)  have  demonstrated  that  their 
multiple  regression  equation  reproduces  measured  oxidation 
states  of  iron  in  natural  melts  quite  well.  Rather  than  go 
through  the  exercise  of  making  the  same  comparisons  as  were 
made  by  Sack  et  al .  (1980),  the  results  of  the  equation 
determined  for  case  E  were  compared  directly  with  the 
results  produced  by  Sack's  equation.  The  comparison  was  made 
for  two  diverse  melt  compositions  (BCR-1  -  Flanagan,  1973 
and  average  rhyolite  -  Nockolds,  1954)  as  a  function  of  both 
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Table  14.  Variables  ranked  according  to  the  relative 

fraction  of  the  total  variance  attributable  to  them  (#1 
=  highest).  The  columns  correspond  to  the  cases 
described  in  the  text. 


Rank 

A 

B 

C 

D 

E 

1 

S  iO  2 

Si02 

S  iO  2 

Si02 

Info 

2 

MgO 

CaO 

CaO 

CaO 

S  iO  2 

3 

FeO 

MgO 

MgO 

MgO 

MgO 

4 

Info  2 

FeO 

FeO 

1/T 

1/T 

5 

CaO 

K20 

A1203 

FeO 

A120 

6 

Al  20  3 

Na  20 

K20 

k2o 

CaO 

7 

Na  20 

1/T 

Info  2 

Na  20 

FeO 

8 

p2o5 

7  nf  o2 

Na  20 

Al  2  0  3 

TiO  2 

9 

TiO  2 

Ti02 

TiO  2 

Ti02 

K20 

10 

1/T 

a  1 2  0  3 

1/T 

Info  2 

MnO 

1  1 

MnO 

MnO 

P20s 

p2o5 

Na  20 

12 

k2o 

p2o5 

MnO 

MnO 

p2o5 

temperature  (fig.  21)  and  oxygen  fugacity  (fig.  22).  Fig.  21 
was  calculated  for  ah  oxygen  fugacity  of  10“ 8  bar.  The 
results  for  BCR-1  show  that  the  equation  determined  in  this 
study  agrees  quite  well  with  that  from  Sack  et  al  .  (1980) 
for  basaltic  liquids.  However,  the  agreement  is  not  as  good 
for  the  rhyolite  composition,  particularly  at  low 
temperatures.  Fig.  22  shows  that  the  two  equations  give 
virtually  the  same  results  as  a  function  of  oxygen  fugacity 
at  1250°C.  All  four  curves  fall  within  the  stipled  area  in 
fig.  22.  Figs.  21  and  22  also  illustrate  that  the  effects  of 
temperature  and  oxygen  fugacity  on  the  oxidation  state  of 
iron  in  the  melt  are  large  when  compared  to  the  effects  of 
bulk  composition. 
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Figure  21.  The  relationship  between  temperature  and  the 
oxidation  state  of  iron  at  an  oxygen  fugacity  of  10‘ 8 
bars  in  silicate  melts  calculated  from  the  case  E 
regression  equation  (Table  13)  for  BCR- 1  (a)  and  average 
rhyolite  (b).  The  results  calculated  from  the  equation 
of  Sack  et  al .  (1980)  for  BCR- 1  (dashed  line)  and  the 
average  rhyolite  (dotted  line)  are  shown  for  comparison. 
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Figure  22.  The  relationship  between  the  oxidation  state  of 
iron  in  silicate  melts  as  calculated  at  1250°C  from  the 
case  E  regression  equation  and  the  equation  of  Sack  et 
a/.  (1980)  for  BCR- 1  and  average  rhyolite.  The  results 
for  both  compositions  and  both  equations  are  included 
within  the  stipled  area. 
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The  multiple  regression  coefficient  modifying  the 
variable  1/T  gives  the  enthalpy  change  of  the  reaction 
(Fudal i ,  1 965 )  : 

Fe01t5^=^Fe0  +  1/40  2  (  19) 

which  describes  the  redox  equilibrium  of  iron  in  silicate 
melts.  This  study  yields  an  enthalpy  change  of  23200  cal  for 
reaction  19.  The  value  determined  from  the  equation  of  Sack 
et  al .  (1980)  is  26200  cal.  These  results  can  be  compared 
with  the  enthalpy  change  calculated  by  Kennedy  (1948)  of 
18500  cal  for  a  basaltic  liquid.  The  differences  between 
Kennedy's  result  and  those  calculated  from  the  results  of 
Sack  et  al .  (1980)  and  this  study  most  likely  reflect  the 
failure  of  the  low  temperature  experiments  done  by  Sack  et 
a/.  (1980)  to  reach  equilibrium,  particularly  for  the  more 
silicic  compositions.  This  may  be  seen  by  comparing  the 
experimental  durations  employed  by  Sack  et  al .  (1980)  with 
those  calculated  from  oxygen  diffusion  data  (fig.  20)  for  a 
temperature  of  1300°C  and  a  basaltic  composition.  The 
diffusivity  of  oxygen  decreases  as  the  silica  content  of  the 
melt  increases  and  as  the  temperature  decreases  (Dunn, 

1982).  Thus,  fig.  20  shows  that  most  of  the  low  temperature 
experiments  done  by  Sack  et  al .  (1980)  were  unlikely  to  have 
reached  equilibrium.  Because  the  experiments  of  Sack  et  al . 
(1980)  were  done  under  reducing  conditions,  the  failure  to 
achieve  equilibrium  will  result  in  an  erroneously  high  ratio 
of  Fe3+  to  Fe2+  being  incorporated  into  the  data  set. 
Furthermore,  since  it  is  most  likely  that  equilibrium  was 


not  achieved  at  the  lower  temperatures,  such  erroneous  data 
will  tend  to  increase  the  apparent  magnitude  of  the 
temperature  dependence  of  the  oxidation  state  of  iron  when 
it  is  expressed  as  the  ratio  of  ferric  to  ferrous  iron. 
Correspondingly,  as  the  enthalpy  change  of  reaction  19  is 
given  by  (bR)  where  b  is  the  regression  coefficient 
modifying  the  variable  1/T  and  R  is  the  gas  constant,  the 
calculated  enthalpy  change  will  be  too  large. 

The  observation  that  the  regression  coefficient 
modifying  1/T  may  be  too  large  limits  the  range  of 
temperatures  over  which  the  regression  equations  determined 
both  in  this  study  and  by  Sack  et  al .  (1980)  may  be  applied. 
Use  of  the  equation  at  temperatures  much  below  1200°C  will 
result  in  erroneously  high  fractions  of  Fe203  being 
calculated.  The  temperature  coefficient  determined  in  this 
study  is  smaller  than  that  determined  by  Sack  et  al .  (1980). 
Therefore,  the  equation  determined  here  seems  to  be  a  better 
choice  if  an  extrapolation  below  1200°C  is  required. 

E.  Conclusion 

Three  different  basaltic  liquids  have  been  equilibrated 
in  terms  of  iron  oxidation  state  at  various  temperatures 
from  1262  to  1437°C,  either  in  air  (24  experiments),  or  at 
an  oxygen  fugacity  of  approximately  10" 8  bars  (3 
experiments).  The  27  experimental  data  points  have  been 
added  to  the  data  base  used  by  Sack  et  al .  (1980)  and  the 
entire  data  set  analyzed  by  multiple  regression  to  yield  a 
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predictive  equation  for  the  oxidation  state  of  iron  in 
silicate  melts. 

The  predictive  equation  is  quite  similar  to  that 
determined  by  Sack  et  a/.  (1980).  However,  the  temperature 
dependencies  of  both  equations  are  probably  too  large. 
Consequently,  extrapolation  to  low  temperatures  will  yield 
erroneously  high  Fe3+/Fe2+  ratios.  In  that  respect,  the 
equation  determined  in  this  study  appears  to  be  preferable 
as  its  temperature  dependence  is  the  smaller  of  the  two. 

While  little  can  be  said  about  the  relative  magnitudes 
of  the  effects  of  changes  in  the  various  oxide  components  on 
the  iron  oxidation  state  of  the  melt,  it  was  observed  that 
the  effect  of  bulk  composition  is  small  relative  to  the 
effects  of  either  temperature  or  oxygen  fugacity. 


VI.  Oxygen  Activity  in  an  Alkali  Basalt  as  a  Function  of 

Temperature 


A.  Introduction 

The  chemical  and  physical  properties  of  iron-bearing 
silicate  melts  are  strongly  affected  by  the  prevailing 
oxygen  fugacity.  Recently,  Tran  and  Brungs  (1980b)  and  Loh 
et  al .  (1981)  have  used  oxygen  activity  probes  similar  to 
those  described  by  Sato  (1971)  to  directly  measure  the 
activity  of  oxygen  in  sodium  disilicate  and  soda-lime  silica 
melts,  respectively.  Their  results  show  that  the  oxygen 
activity  in  simple  silicate  melts  containing  polyvalent 
cations  can  change  very  sharply  with  changes  in  temperature. 
Their  results  prompted  this  study  of  the  feasibility  of 
using  oxygen  activity  probes  to  directly  measure  the 
activity  of  oxygen  in  natural  rock  melts. 

This  study  reports  the  results  of  oxygen  activity 
measurements  in  an  alkali  basalt  melt  over  the  temperature 
range  1257-1373°C.  The  oxygen  potential  measurements  were 
made  with  an  Y203  doped  Zr02  electrolyte  cell  (Arculus  and 
Delano,  1981).  The  half  cell  reactions  upon  which  cells  of 
this  type  are  based  are: 

20  2  "  ;=^  02  (ref  )  +  4e"  (20) 

at  the  reference  electrode  and: 

02  (melt )  +  4e-^=^  202  '  (21) 

at  the  measuring  electrode.  The  overall  reaction  is: 

02 (melt ) ~  s 0  2 ( ref ) 
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(22) 
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The  free  energy  change  for  reaction  22  is: 

AG(rxn)  =  AG°(rxn)  +  RT7n(a02 ( ref )/a02 (melt ) )  (23) 

The  half  reaction  which  occurs  in  the  melt  is  given  by 
reaction  21.  As  four  electrons  are  transferred  in  the 
reaction  the  relationship  between  the  half  cell  potential 
and  the  free  energy  change  of  the  reaction  is: 

AG(rxn)  =  -4 FE  (24) 

where  E  is  the  cell  emf  and  F  is  the  Faraday  constant. 
Substituting  eqn.  24  into  eqn .  23,  and  noting  that  AG 0  = 

-4 FE°  (where  E°  is  the  standard  half  cell  potential  at  unit 
activity  which  goes  to  zero  for  a  pure  oxygen  reference), 
gives : 

£  =  RT  7  n ( aO 2 ( r  e  f ) /aO 2 (melt) )  (25) 

Eqn.  25  holds  in  the  absence  of  ionic  conduction  in  the 
electrolyte  cell. 

B.  Experimental  procedure 

The  experiment  was  done  on  a  20  gram  batch  of  alkali 
basalt  from  Kettle  River  B.C.  (KR-13;  Table  15).  A  sample  of 
the  rock  was  crushed  in  a  tungsten  carbide  swing  mill  and  a 
split  of  the  rock  powder  was  placed  in  a  Pt  crucible  and 
fused  for  48  hours  at  1300°C  in  air.  The  resultant  glass  was 
crushed  (several  chips  were  mounted  for  microprobe 
analysis),  returned  to  the  Pt  crucible,  mounted  on  a 
pedestal,  and  raised  into  the  hot  zone  of  an  electrically 
heated  vertical  tube  furnace  (Deltech  DT-31VT).  The  sample 
was  held  at  1250°C  for  20  hours  prior  to  starting  the 
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Figure  23.  Schematic  drawing  of  the  extper  imental 

configuration,  showing  the  oxygen  electrolyte  cell  and 
the  position  of  the  sample  melt. 


Table  15.  Composition  of  the  alkali  basalt.  Wt .  %  oxide 
determined  by  electron  microprobe  analysis  of  glass 
chips . 


S  i02 

48.39 

Ti02 

2.41 

A1203 

15.91 

MgO 

5.96 

CaO 

9.62 

FeO1 

11.35 

MnO 

0.25 

Na  20 

3.98 

K20 

1 .22 

p2o5 

0.51 

Total 

99.60 

-  Total  iron  as  FeO 


experiment  to  assure  oxidative  equilibrium  in  the  melt. 

The  microprobe  analysis  was  done  on  the  University  of 
Alberta  ARL-SEMQ  electron  microprobe  using  wavelength 
dispersive  methods.  Suitable  standards  were  used  to  allow 
reduction  of  the  raw  data  by  the  method  of  Bence  and  Albee 
(1968)  using  the  alpha  factors  of  Albee  and  Ray  (1970).  The 
analysis  of  the  starting  material  is  given  in  Table  15. 

The  oxygen  probe  used  in  this  study  was  the  same  as 
those  used  by  Arculus  and  Delano  (1981).  However,  several 
modifications  were  made  to  the  experimental  configuration. 
Only  one  probe  was  used  and  oxygen  gas  was  used  to  provide 
the  reference.  The  measuring  electrode  consisted  of  a  Pt 
wire  with  a  6  cm  length  of  high  purity  A1203  tube  sintered 
onto  it  approximately  3  cm  from  one  end.  The  exposed  end  of 
the  Pt  wire  was  attached  directly  to  the  z i rconia/yttr ia 
electrolyte  (fig.  23).  The  Al203  tube  extends  into  the  melt 
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when  the  probe  is  in  use,  preventing  the  establishment  of 
the  triple  junction  Pt-air-melt .  This  assures  that  the 
oxygen  activity  in  the  melt  rather  than  at  the  melt/gas 
interface  is  measured.  Tran  and  Brungs  (1980b)  noted  that, 
for  conventional  Sato  type  oxygen  probes  (Sato,  1971), 
diffusion  of  oxygen  through  the  electrolyte  resulted  in 
anomalously  high  oxygen  potentials  near  the  oxygen  probe.  To 
alleviate  that  problem  they  separated  the  measuring 
electrode  from  the  the  oxygen  probe.  That  configuration  was 
not  used  in  this  study  because  it  was  thought  highly 
unlikely  that  any  oxygen  would  diffuse  through  the  one 
centimeter  length  of  electrolyte  in  the  probe  used. 
Temperatures  were  measured  with  a  Pt-Pt87Rh13  thermocouple 
placed  inside  the  oxygen  probe  (fig.  23). 

The  emf  output  of  the  electrolyte  cell  was  measured 
with  both  a  digital  voltmeter  (10  megaohm  input  impedance) 
and  a  chart  recorder  (8  megaohm  input  impedance).  All  of  the 
emf' s  reported  here  were  measured  with  the  digital 
voltmeter.  The  digital  voltmeter  was  also  used  in  an  attempt 
to  measure  the  cell  current.  No  measurement  was  obtained  on 
the  most  sensitive  scale  (200  microamps  full  scale),  thus 
satisfying  the  requirement  that  the  cell  current  be  zero. 

The  experiment  was  initiated  by  placing  the  electrolyte 
cell  into  the  melt.  The  melt  was  then  stepped  to  1257°C  and 
held  at  that  temperature  until  the  cell  emf  appeared  to 
stabilize.  The  melt  temperature  was  then  increased  by  steps 
to  1373°C.  The  intermediate  temperatures  were  1286,  1316, 
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and  1345°C.  At  each  step,  the  melt  was  held  until  the  cell 
emf  appeared  to  stabilize.  After  reaching  1373°C  the  system 
was  cooled  back  to  1257°C  using  the  same  temperature  steps. 
After  the  system  appeared  to  stabilize  at  1257°C  another 
heating  cycle  was  begun,  but  was  only  carried  to  1316°C 
where  the  experiment  was  terminated. 

At  the  end  of  the  experiment  the  electrolyte  cell  was 
withdrawn  from  the  melt.  Examination  of  the  cell  showed 
that,  except  for  a  thin  coating  of  basalt  melt,  the  cell  was 
apparently  unaffected  by  its  76  hours  immersion  in  the 
basalt  melt.  However,  the  lack  of  a  visibly  obvious  reaction 
does  not  mean  that  the  cell  did  not  react  with  the  melt  to  a 
limited  extent.  The  evidence  for  limited  reaction  between 
the  electrolyte  cell  and  the  melt  is  presented  below. 

C.  Results  and  discussion 

The  results  of  the  experiment  are  summarized  in  Table 
16  and  shown  graphically  in  figs.  24  and  25.  At  the 
initiation  of  the  experiment  the  emf  of  the  electrolyte  cell 
stabilized  within  10  minutes.  The  temperature  was  then 
abruptly  increased  to  1257°C.  As  soon  as  the  temperature  of 
the  system  began  to  increase,  the  cell  emf  began  to  decrease 
reaching  a  minimum  after  8  to  10  minutes,  when  the  system 
temperature  reached  a  maximum.  After  the  temperature 
stabilized,  the  cell  emf  began  a  slow  asymptotic  approach  to 
its  final  (equilibrium)  value  (fig.  24).  Similar  behavior 
was  observed  for  each  temperature  increase.  Apparent 
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Figure  24.  Plot  of  the  emf  output  of  the  oxygen  electrolyte 
cell  as  a  function  of  time.  Two  cooling  (a  and  b)  steps 
and  one  heating  step  (c)  are  shown.  The  temperature 
stabilized  at  the  new  setting  within  10  minutes  after 
each  step. 
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Figure  25.  Plot  of  lna(02)  versus  1000/T  (K)  showing  the 

experimental  brackets  and  the  best  fit  line  to  the  data 
(eqn.  27).  The  solid  triangles  are  values  used  to 
determine  the  best  fit  line.  Open  triangles  are  data 
points  which  were  not  used.  The  reasons  for  excluding 
those  points  are  discussed  in  the  text. 
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Table  16.  Oxygen  probe  readings  after  apparent  stabilization 
and  calculated  oxygen  activities. 


T°C 

Mode 

t  (hr ) 

emf ( v ) 

ao2 

1257 

H 

1.6 

0.0501 

0.216 

1257 

C 

16.0 

0.0540 

0.194 

1286 

H 

2.2 

0.0465 

0.251 

1286 

C 

7.5 

0.0541 

0.200 

1286 

H 

10.0 

0.0480 

0.240 

1316 

H 

3.5 

0.0458 

0.263 

1316 

C 

9.5 

0.0516 

0.220 

1316 

H 

16.3 

0.0483 

0.244 

1345 

H 

2.5 

0.0392 

0.325 

1345 

C 

4.5 

0.0491 

0.245 

1373 

H 

2.3 

0.0371 

0.351 

Mode  -  Indicates  whether  the  temperature  was  approached  by 
heating  (H)  or  cooling  (C)  the  melt. 


stabilization  times  ranged  from  as  little  as  6  to  longer 
than  16  hours  and  were  not  systematically  related  to  the 
system  temperature.  The  cell  response  was  the  same,  but  in 
the  opposite  sense,  when  the  temperature  was  decreased.  When 
the  melt  temperature  was  decreased  the  cell  emf  rapidly 
increased,  passed  through  a  maximum  at  the  minimum 
temperature  and  then  asymptotically  approached  equilibrium 
(fig.  24).  The  maxima  and  minima  observed  in  the  cell  emf 
correspond  to  oxygen  activities  as  low  as  0.1  and  as  high  as 
0.6.  Those  are  significant  variations  relative  to  the 
equilibrium  value  of  0.21.  While  the  magnitude  of  the 
changes  to  be  expected  at  low  oxygen  fugacity  is  unknown,  if 
those  changes  are  of  the  same  relative  magnitude  as  they  are 
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at  high  oxygen  fugacities,  the  redox  equilibria  in  the  melt 
could  be  greatly  affected. 

Eqn .  25  shows  that  the  oxygen  activity  in  the  melt  is 
inversely  related  to  the  electrolyte  cell  emf.  Thus  the 
decrease  in  the  cell  emf  upon  heating  corresponds  to  an 
increase  in  oxygen  activity,  and  the  emf  increase  upon 
cooling  indicates  a  decrease  in  oxygen  activity.  Such  sharp 
changes  in  oxygen  activity  upon  heating  and  cooling  (fig. 

24)  can  be  attributed  to  rapid  oxidation/reduction  reactions 
in  the  melt  (Loh  et  al . ,  1981). 

The  only  cation  present  in  significant  proportions  in 
basaltic  melts  which  undergoes  a  valence  change  under  the 
conditions  of  the  experiment,  is  iron.  The 

oxidation/reduction  behavior  of  iron  in  silicate  melts  can 
be  expressed  by  equation  26  (Fudali,  1965  and  Sack  et  al . , 
1980)  : 


FeO  ,.5  - - ±FeO  +  1/4  02  (26) 

where  FeO,.s  is  ferric  iron  and  FeO  is  ferrous  iron.  Equation 
26  shows  that  the  reduction  of  ferric  iron  produces  oxygen 
and  thus  causes  an  increase  in  the  activity  of  oxygen  in  the 
melt.  Conversely,  the  oxidation  of  ferrous  iron  consumes 
oxygen  and  decreases  the  activity  of  oxygen.  Many  workers 
(e.g.  Sack  et  a/.,  1980;  Thornber  et  al . ,  1980;  Fudali, 

1965,  and  Kennedy,  1948)  have  observed  that  the  proportion 
of  ferrous  iron  in  basaltic  liquids  increases  with 
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increasing  temperature.  Thus,  it  is  reasonable  to  equate  the 
abrupt  increase  in  oxygen  activity  observed  with  increasing 
temperature  with  the  production  of  oxygen  through  the 
reduction  of  Fe3+  to  Fe2+  and  vice  versa  for  decreasing 
temperature.  The  abrupt  nature  of  the  change  in  oxygen 
activity  with  changes  in  temperature,  indicates  that  the 
oxidation/reduction  of  iron  in  silicate  melts  is  fast. 
Similarly,  Semkow  et  a/.  (1982)  have  observed  that 
oxidation/reduction  reactions  for  Ni,  Co,  and  Zn  in  diopside 
composition  melts  are  rapid. 

The  slow  decrease  in  the  activity  of  oxygen  (increase 
in  cell  emf )  following  a  temperature  increase  is  due  to  the 
diffusion  of  oxygen  (produced  by  the  reduction  of  ferric 
iron)  out  of  the  system.  The  reverse  is  the  case  following  a 
temperature  decrease.  The  time  required  for  the  system  to 
reach  the  equilibrium  oxygen  activity  is  dependent  on  the 
geometry  and  dynamics  (whether  or  not  the  system  is 
convecting)  of  the  system.  Thus,  it  seems  unlikely  that 
equilibrium  was  achieved  at  any  time  during  the  course  of 
the  experiment. 

Fig.  25  is  a  plot  of  the  log  of  the  oxygen  activity, 
determined  with  the  probe  after  stabilization,  versus  1/T. 
The  brackets  determined  during  the  heating  and  cooling 
cycles  are  shown.  The  line  in  fig.  25  is  a  least  squares  fit 
to  the  data  points  shown  as  solid  triangles.  Data  points 
shown  as  open  triangles  were  not  included  in  the  regression 
analysis  because  their  inclusion  caused  the  line  to  fall 
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outside  one  or  more  of  the  experimental  brackets.  The 
principle  criterion  used  to  exclude  those  points  is  that  the 
process  responsible  for  the  non-equilibrium  behavior  has  a 
simple  exponential  temperature  dependence  over  the 
temperature  interval  studied.  In  that  case  the  line  shown  in 
fig.  25  must  be  a  straight  line.  The  equation  of  the  line  in 
fig.  25  is: 

lna(02 )  =  ( 1 . 7 1 ± 1 .46)-(5030±2300)/T  r2=0.49  (27) 

The  non-zero  slope  of  the  line  given  in  equation  27  shows 
that  the  oxygen  probe  measured  anomalously  high  oxygen 
activities  relative  to  the  imposed  oxygen  fugacity  at  higher 
temperatures.  That  observation  can  be  explained  by 
permeation  of  oxygen  through  either  the  electrolyte  or  the 
alumina  tube  portions  of  the  probe.  Tran  and  Brungs  (1980b) 
report  similar  behavior,  but  only  when  highly  reduced  melts 
were  examined.  Due  to  the  design  differences  between  the 
oxygen  probes  used  by  Tran  and  Brungs  (1980b)  and  those  used 
in  this  study  (fig.  23)  and  the  fact  that  an  oxidized  melt 
was  studied,  it  is  likely  that  the  anomolous  activities 
observed  are  not  due  to  the  permeation  of  oxygen  through  the 
probe  walls.  It  seems  likely,  that  at  the  higher 
temperatures  there  was  some  solution  of  the  Al203  from  the 
alumina  tube  into  the  melt.  Solution  of  Al203  would  promote 
polymerization  of  the  melt  (Mysen  et  al . ,  1980d),  thus 
increasing  the  activity  of  oxygen  in  the  vicinity  of  the 
probe.  That  problem  does  not  have  a  significant  effect  on 
the  conclusions  of  this  study,  but  must  be  taken  into 
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consideration  if  oxygen  probes  of  the  type  used  in  this 
study  are  to  be  used  for  quantitative  studies  in  natural 
rock  melts. 

D.  Petrologic  implications 

This  study  has  shown  that,  upon  heating  and  cooling, 
the  oxidation  state  of  iron  in  the  melt  changes  very  rapidly 
either  releasing  (upon  the  reduction  of  iron)  or  consuming 
(upon  the  oxidation  of  iron)  oxygen  in  the  process.  The 
release  or  consumption  of  oxygen  rapidly  changes  the 
oxidation  state  of  the  melt.  When  a  melt  cools,  the  iron  in 
the  melt  becomes  increasingly  oxidized  in  response  to  the 
temperature  change  (Kennedy,  1948;  Fudali,  1965;  and  Sack  et 
al . ,  1980).  Simultaneously,  the  activity  of  oxygen  in  the 
melt  decreases,  creating  disequilibrium  conditions.  In  an 
open  system,  oxygen  can  diffuse  into  the  melt  from  the 
surroundings  and  increase  the  oxygen  activity  back  to  the 
equilibrium  value.  Open  system  conditions  of  this  sort  exist 
in  gas  buffered  experiments.  If  sufficient  time  is  allowed 
for  oxygen  diffusion  to  restore  the  equilibrium  oxygen 
activity,  then  the  equilibrium  iron  oxidation  state  in  the 
melt  is  attained.  In  crystallization  experiments,  the 
samples  are  generally  cooled  from  super 1 iquidus  temperatures 
to  the  run  temperature.  That  will  cause  the  oxygen  activity 
to  be  depressed  below  the  equilibrium  value  for  the 
conditions.  Such  k inet ically- imposed  reducing  conditions  may 
induce  the  crystallization  of  phases  which  will  not  be  in 
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equilibrium  when  the  oxygen  activity  in  the  melt 
equilibrates.  If  such  crystalline  phases  do  not  dissolve 
back  into  the  melt,  they  can  only  reach  equilibrium  with  the 
melt  through  solid  state  diffusion. 

Heating  or  cooling  a  melt  under  closed  system 
conditions  (relative  to  oxygen  transfer)  also  increases  or 
decreases  the  oxygen  activity  in  the  melt.  However,  in  the 
closed  system  case,  oxygen  cannot  diffuse  into  or  out  of  the 
system.  The  result,  upon  cooling,  is  that  the  oxygen 
activity  becomes  increasingly  reducing.  If  the  melt  was 
originally  quite  reduced,  it  is  conceivable  that  the  system 
could  autoreduce  itself  into  the  Fe2+  or  even  the  Fe°  field. 
Such  a  situation  could  lead  to  the  crystallization  of  phases 
enriched  in  Fe 2  + ,  relative  to  the  composition  expected  from 
the  bulk  composition  and  estimated  equilibrium  oxygen 
fugacity  of  the  melt.  In  extreme  conditions,  this  process 
could  even  lead  to  the  precipitation  of  metallic  iron  from 
the  melt  as  is  observed  for  lunar  basalts  (Levinson  and 
Taylor,  1971).  Closed  system  conditions  of  this  sort  are 
easy  to  conceive  of  for  natural  melts.  The  case  of  a  basalt 
erupted  onto  the  surface  of  the  earth  is  one  instance  where 
the  melt  is  rapidly  cooled  and  closed  system  conditions  may 
be  provided  by  the  chilled  margin  of  the  flow. 

Autoreduction  of  the  melt  by  the  process  described 
above  may  also  affect  the  volatile  phase  speciation  in  the 
melt.  If  a  melt  is  cooled  quickly  over  a  large  temperature 
range,  as  is  the  case  for  sea  floor  basalts,  it  is 
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conceivable  that  some  of  the  C02  dissolved  in  the  melt  could 
be  reduced  to  CO  or  even  to  graphite  (which  could  then 
precipitate).  A  large  decrease  in  the  activity  of  oxygen  in 
the  melt  would  also  promote  the  dissociation  of  water  to  H2 
and  02 ,  thus  lowering  the  apparent  water  content  of  the  melt 
(if  the  H 2  diffuses  out  of  the  system).  This  process  can 
account  for  the  presence  of  elemental  carbon  in  sea  floor 
basalts  (Mathez  and  Delaney,  1981)  as  well  as  for  the  low 
water  contents  observed  for  some  sea  floor  basalts  (Muenow 
et  al . ,  1979).  Furthermore,  because  of  this  effect,  the  iron 
oxidation  state  of  quenched  basalts  may  not  be  a  good 
indicator  of  the  prevailing  oxygen  fugacity  in  the  melt 
prior  to  eruption. 

E.  Conclusion 

This  study  has  shown  that  the  oxidation  state  of  a 
basaltic  melt  changes  rapidly  in  response  to  temperature 
changes.  Cooling  of  a  melt  results  in  an  immediate  decrease 
in  the  activity  of  oxygen  as  a  consequence  of  the  rapid 
oxidation  of  Fe2+  to  Fe3+.  In  a  system  closed  to  oxygen 
transfer,  the  decrease  in  oxygen  activity  may  result  in  the 
quenched  melt  being  out  of  oxidative  equilibrium  with  the 
crystalline  phase  assemblage.  Consequently,  the  oxidation 
state  of  quenched  melts  (glasses)  may  not  reflect  the 
oxidation  state  of  the  magma  from  which  the  glass  was 
quenched . 
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VI I .  Conclusion 

A.  Summary  of  conclusions 

The  five  studies  presented  in  this  thesis  illustrate 
the  value  of  studying  oxygen  in  silicate  melts.  Chapter  2 
showed  that  the  principal  species  contributing  to  the 
diffusivity  of  oxygen  is  no  larger  than  the  Si044-  monomer, 
even  in  melts  which  are  thought  of  as  highly  polymerized.  It 
was  also  observed  that  oxygen  conforms  to  the  compensation 
law  for  cationic  diffusion  in  silicate  melts.  That 
observation  suggests  that  the  diffusing  species  is 
"cation-like"  (O2").  The  magnitude  of  oxygen  diffusivity  was 
found  to  be  similar  to  the  magnitudes  of  divalent  cation 
diffusivity,  which  further  supports  the  suggestion  that 
oxygen  diffusivity  receives  a  significant  contribution  from 
02-  anions. 

Chapter  3  showed  that  the  relationship  between  oxygen 
diffusion  in  basaltic  liquids  and  pressure  is  not  at  all 
like  the  pressure  dependences  observed  for  cations  in 
silicate  melts.  Furthermore,  the  observation  of  negative 
activation  volumes  for  oxygen  diffusion  argues  strongly  that 
variations  in  the  diffusivity  of  oxygen  with  pressure  are 
not  related  to  changes  in  the  free  volume  of  the  melt.  The 
correspondence  between  abrupt  changes  in  oxygen  diffusivity 
and  changes  in  the  liquidus  mineralogy  of  the  melts  studied, 
suggests  that  the  pressure  dependence  of  oxygen  diffusion  is 
related  to  changes  in  the  anionic  speciation  of  the  melt. 
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The  pressure  dependence  of  oxygen  diffusion  is  not  related 
to  the  pressure  dependence  of  viscosity  for  the  melts 
studied.  Therefore,  the  correspondence  between  viscous  flow 
and  oxygen  diffusivity  suggested  by  Oishi  et  al .  (1975)  and 
Yinnon  and  Cooper  (1980)  may  be  fortuitous.  Oxygen 
di f f usivit ies  in  basaltic  liquids  are  of  the  same  magnitude 
as  divalent  cation  di f f us ivi t ies .  Therefore,  in  cases  of 
diffusion  controlled  or  slower  crystal  growth,  crystals  will 
not  be  out  of  oxidative  equilibrium  with  the  body  of  the 
melt . 

Chapter  4  demonstrated  that  the  structure  of  Pb2Si04 
melts  depends  on  thermal  history.  The  melts  become 
increasingly  polymerized  (NBO/T  decreases)  with  either 
decreasing  cooling  rate,  or  increasing  duration  of 
isothermal  soaking.  The  presence  of  silicate  anions  larger 
than  Si044~  in  lead  orthosilicate  melt  requires  a 
corresponding  proportion  of  02~  anions  in  the  melt. 

Chapter  5  presents  an  improved  predictive  equation  for 
calculating  either  the  oxidation  state  of  iron  in  a  melt 
(given  the  temperature  and  the  oxygen  fugacity),  or  the 
oxygen  fugacity  of  a  melt  (given  the  iron  oxidation  state 
and  the  temperature).  The  oxidation  state  of  iron  in 
silicate  melts  was  shown  to  depend  much  more  on  temperature 
and  oxygen  fugacity  than  on  bulk  composition. 

Chapter  6  presented  new  data  on  the  change  of  the 
oxidation  state  of  a  basaltic  melt  with  changing 
temperature.  The  oxygen  activity  of  a  basalt  melt  was  shown 
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to  decrease  rapidly  on  cooling  in  response  to  oxidation  of 
Fe2+.  Reequilibration  of  the  oxygen  activity  then  occurred 
by  the  diffusion  of  oxygen  into  the  melt  from  the 
surroundings.  In  the  case  of  a  system  closed  to  oxygen 
transfer,  this  process  provides  a  possible  explanation  for 
the  occurrence  of  reduced  phases  in  an  otherwise  apparently 
more  oxidized  system. 

B.  Concluding  remarks 

The  results  presented  in  chapters  2,  3,  and  4  suggest 
that  02"  is  present  not  only  in  orthosilicate  melts,  but 
also  in  much  more  polymerized  melts.  That  observation 
requires  that  models  of  silicate  melt  structure  which  do  not 
incorporate  the  02~  anion  be  reassessed.  The  presence  of  02" 
in  relatively  polymerized  silicate  melts  allows  the  presence 
of  a  greater  variety  of  anionic  species  in  the  melt  than  is 
the  case  if  02"  is  not  present.  Furthermore,  Mysen  et  al . 

( 1 9 8 Od )  argue  that  NBO/T,  as  they  calculate  it,  is 
definitive  of  the  structure  of  the  melt  and  depends  only  on 
bulk  composition.  The  results  presented  here  do  not  agree 
with  that  assertion.  The  presence  of  oxygen  anion  (02“)  in 
silicate  melts  and  the  suggestion  that  it  participates  in 
reactions  with  anionic  polymers  (chapter  3)  implies  that  the 
structure  of  the  melt  is  dependent  on  oxygen  fugacity. 
Finally,  the  results  of  this  thesis  have  shown  that  oxygen 
can  act  as  a  very  sensitive  probe  into  the  structure  of 


silicate  melts. 
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